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ABSTRACT 

Recent  advances  in  signal  processing  and  remote  sensing  have  highlighted  the 
importance  of  bistatic  radar  systems  for  the  purposes  of  environmental  monitoring, 
surveillance,  and  tracking  radar.  The  calibration  of  such  systems  has  been  problematic- 
much  more  so  than  similar  monostatic  systems,  primarily  as  a  result  of  the  lack  of 
reference  objects  suitable  for  calibrating  at  any  given  bistatic  angle.  This  research  deals 
with  the  problems  of  calibrating  full-polarimetric  laboratory-environment  bistatic  radar 
systems,  including  the  lack  of  suitable  calibration  targets  and  procedures,  and  operational 
considerations  such  as  alignment  and  mounting.  Several  popular  bistatic  calibration 
techniques  are  classified,  evaluated,  and  comparisons  are  made  between  the  relative 
merits  of  various  calibration  objects.  The  analysis  addresses  sensitivity  to  target 
alignment  error,  sensitivity  to  polarization  impurity,  and  ease  of  implementation.  Both 
theoretical  concepts  and  practical  considerations  are  discussed,  based  on  measurements 
accomplished  at  the  European  Microwave  Signature  Laboratory  (EMSL)  of  the  Joint 
Research  Center  (JRC)  in  Ispra,  Italy.  Significant  gains  in  co-polarized  channel  accuracy 
and  cross-polarization  purity  are  realized  with  calibrations  that  utilize  the  complete 
system  distortion  model,  and  these  conclusions  are  discussed  in  detail. 
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THE  CALIBRATION  OF  BISTATIC  RADAR  CROSS  SECTION  MEASUREMENTS 

1  INTRODUCTION 

Radar  calibration  is  a  procedure  used  to  compensate  for  irregularities  in  the 
hardware  or  measurement  systems  and  to  account  for  systematic  error-producing 
phenomenon  present  in  the  measurement.  A  good  radar  calibration  technique  takes  into 
consideration  the  uniqueness  of  each  measurement  system  as  well  as  its  application. 
Applications  for  radar  calibration  come  in  two  main  configurations:  monostatic  and 
bistatic. 

Figure  1  shows  the  difference  between  the  monostatic  and  bistatic  geometries. 
Monostatic  radar  is  defined  as  a  system  that  transmits  a  radio-frequency  electromagnetic 
wave  and  receives  the  scattered  wave  at  the  same  location  for  the  purposes  of  detection 
and  ranging  of  objects  within  the  antenna’s  illumination  area.  Typically,  a  monostatic 
radar  will  employ  a  single  antenna  that  is  switched  between  transmit  and  receive  modes. 
Often,  however,  two  closely  placed  antennas  are  used;  one  for  transmission  and  one  for 
reception.  This  mitigates  problems  due  to  cross  talk  between  the  Tx/Rx  channels,  and 
allows  each  antenna  to  be  designed  specifically  for  transmitting  or  receiving. 
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Monostatic 


Figure  1:  Monostatic  and  Bistatic  Configurations  (Eigel,  1999) 

The  proper  calibration  of  radars  is  highly  dependent  on  the  type.  A  monostatic 
calibration  is  one  that  is  designed  to  calibrate  based  on  the  considerations  for  monostatic 
measurements.  This  type  of  calibration  is  relatively  straightforward  and  the  procedures 
follow  well-known  heuristics.  Bistatic  calibration,  however,  is  much  more  complicated. 
Additional  measurement  concerns  such  as  the  presence  of  more  than  one  antenna,  cabling 
and  non-linearities  in  antenna  phase  response  over  the  angular  sweep,  and  calibration 
object  selection  make  this  type  of  calibration  difficult.  Much  less  is  known  about  bistatic 
calibration,  though  bistatic  radar  itself  has  been  in  development  almost  as  long  as 
monostatic  radar.  It  is  a  direct  result  of  recent  advantages  in  signal  processing 
capabilities  that  bistatic  radar,  and  thus  bistatic  calibration,  have  become  popular. 

The  use  of  bistatic  radar  has  its  advantages  over  the  monostatic  case.  The 
electromagnetic  scattering  signature  (most  often  referred  to  as  radar  cross-section,  or 
RCS)  can  be  determined  in  a  much  more  complete  way  as  opposed  to  monostatic  radar. 
Monostatic  radar  can  only  observe  the  backscatter  of  an  object,  or  the  radiation  that 
returns  from  the  object  directly  back  in  the  direction  of  the  source.  We  know  from 


electromagnetic  scattering  theory  that  objects,  in  general,  radiate  in  all  directions,  not  just 
the  backscatter  direction.  A  bistatic  measurement  can  observe  how  the  object  scatters  in 
every  direction  in  a  plane  that  intersects  the  object,  for  a  given  look  angle.  The 
disadvantage  of  bistatic  radar  is  directly  coupled  to  its  main  advantage.  A  complete 
bistatic  measurement  can  require  hundreds  of  times  the  amount  of  data  as  a  monostatic 
measurement  for  a  given  look  angle.  Analysis  of  the  data  becomes  computationally 
expensive  and  much  less  intuitive  than  in  the  monostatic  case. 

The  inherent  advantage  to  bistatic  radar  has  a  direct  impact  the  future  of  U.S.  Air 
Force  operations.  Currently,  low  observable  (LO)  aircraft  have  been  designed  to  have  a 
low  backscatter  cross-section.  It  can  be  shown  that  using  aircraft  shaping  techniques  to 
reduce  the  backscatter  return  of  LO  aircraft  tends  to  increase  the  bistatic  RCS  of  the  same 
aircraft  (Tuley,  Alexander,  1995),  thereby  increasing  its  probability  of  detection.  With 
the  increased  use  of  bistatic  radar  for  aircraft  detection,  the  threat  to  U.S.  LO  aircraft 
increases  accordingly.  It  is  for  this  reason  that  it  is  necessary  to  understand  bistatics  and 
how  best  it  can  be  used  for  the  detection  of  low  observables.  A  logical  first  step  to  this 
end  would  be  to  develop  a  way  to  maximize  the  accuracy  with  which  one  can  measure 
the  bistatic  scattering  characteristics  of  a  given  object,  ergo  the  need  for  efficient  and 
precise  calibration  techniques. 

Two  major  problems  exist  that  make  bistatic  calibration  difficult.  First  of  all,  the 
reference  objects  preferred  in  monostatic  calibration  such  as  dihedrals  and  trihedrals  are 
not  appropriate  for  bistatic  measurements.  For  monostatic  backscatter  measurements, 
dihedrals  and  trihedrals  have  a  high  radar  return,  which  ensures  data  collection  within  the 
dynamic  range  of  the  radar  system.  In  addition,  the  theoretical  backscatter  of  these 
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objects  can  be  easily  and  precisely  calculated  based  on  specular  scattering.  For  any 
measurement  other  than  backscatter,  i.e.  bistatic  measurements,  these  objects  are  no 
longer  desirable.  The  response  of  the  object  is  no  longer  dominated  by  specular 
scattering,  and  as  a  result  the  theoretical  response  is  much  more  difficult  to  calculate 
precisely,  and  the  radar  return  will  be  much  lower  at  some  bistatic  angles,  which  makes 
accurate  calibration  difficult.  Other  objects  such  as  spheres,  cylinders,  and  plates  are 
often  used  for  bistatic  reference  objects,  since  their  solution  can  be  calculated  more 
easily,  and  the  radar  return  is  sufficiently  high  for  most  bistatic  angles. 

The  second  problem  is  that  the  theoretical  RCS  of  the  reference  objects  for  a 
bistatic  calibration,  which  are  necessary  to  compare  with  the  measurement  in  order  to 
calibrate,  are  often  more  rigorous  than  in  a  monostatic  calibration.  In  a  complete,  full- 
polarimetric  calibration  for  instance,  usually  three  objects  are  used  and  each  must  have 
precisely  known  scattering  characteristics,  which  often  means  several  days  of  moment 
method  calculations  on  a  fast  computer. 

This  work  addresses  the  problem  of  finding  suitable  calibration  objects  for  bistatic 
calibration,  and  also  evaluates  the  various  techniques  that  are  currently  used  and  are 
available  in  current  literature.  Calibration  techniques  detailed  in  (Alexander  et  al.,  1995), 
(Jersak,  1993),  (Kahny,  et  al.,  1992),  (Mortensen,  1995),  (Ulaby,  Elachi,  1990),  and 
(Whitt  et  al,  1991)  were  performed  at  the  European  Microwave  Signature  Laboratory 
(EMSL)  of  the  European  Commission  Joint  Research  Center  (JRC)  in  Ispra,  Italy.  The 
EMSL  was  chosen  as  the  location  to  perform  the  measurements  because  of  its  capability 
and  availability  to  produce  far-field,  fully  polarimetric,  bistatic  measurements  on  a 
variety  of  objects.  The  precision  of  this  facility  is  of  paramount  importance  for  the 
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objects  can  be  easily  and  precisely  calculated  based  on  specular  scattering.  For  any 
measurement  other  than  backscatter,  i.e.  bistatic  measurements,  these  objects  are  no 
longer  desirable.  The  response  of  the  object  is  no  longer  dominated  by  specular 
scattering,  and  as  a  result  the  theoretical  response  is  much  more  difficult  to  calculate 
precisely,  and  the  radar  return  will  be  much  lower  at  some  bistatic  angles,  which  makes 
accurate  calibration  difficult.  Other  objects  such  as  spheres,  cylinders,  and  plates  are 
often  used  for  bistatic  reference  objects,  since  their  solution  can  be  calculated  more 
easily,  and  the  radar  return  is  sufficiently  high  for  most  bistatic  angles. 

The  second  problem  is  that  the  theoretical  RCS  of  the  reference  objects  for  a 
bistatic  calibration,  which  are  necessary  to  compare  with  the  measurement  in  order  to 
calibrate,  are  often  more  rigorous  than  in  a  monostatic  calibration.  In  a  complete,  full- 
polarimetric  calibration  for  instance,  usually  three  objects  are  used  and  each  must  have 
precisely  known  scattering  characteristics,  which  often  means  several  days  of  moment 
method  calculations  on  a  fast  computer. 

This  work  addresses  the  problem  of  finding  suitable  calibration  objects  for  bistatic 
calibration,  and  also  evaluates  the  various  techniques  that  are  currently  used  and  are 
available  in  current  literature.  Calibration  techniques  detailed  in  (Alexander  et  al.,  1995), 
(Jersak,  1993),  (Kahny,  et  al.,  1992),  (Mortensen,  1995),  (Ulaby,  Elachi,  1990),  and 
(Whitt  et  al,  1991)  were  performed  at  the  European  Microwave  Signature  Laboratory 
(EMSL)  of  the  European  Commission  Joint  Research  Center  (JRC)  in  Ispra,  Italy.  The 
EMSL  was  chosen  as  the  location  to  perform  the  measurements  because  of  its  capability 
and  availability  to  produce  far-field,  fully  polarimetric,  bistatic  measurements  on  a 
variety  of  objects.  The  precision  of  this  facility  is  of  paramount  importance  for  the 
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accurate  evaluation  of  these  sometimes-sensitive  measurements  and  calibration 
techniques. 

The  scope  of  the  project  provided  for  a  performance  analysis  of  existing  techniques 
based  on  the  following  factors: 

•  Sensitivity  of  calibration  to  reference  object  alignment  error 

•  Sensitivity  to  antenna  polarization  impurity 

•  Robustness  and  repeatability  of  the  technique 

•  Applicability  to  fixed-bistatic  or  true  bistatic  measurements 

•  Ease  and  efficiency  of  implementation 

The  following  pages  introduce  the  theoretical  and  practical  concepts  of  bistatic 
measurements  and  calibration,  and  also  document  the  process,  analysis,  results,  and 
conclusions  of  the  experimentation. 
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2  BACKGROUND 


2.1  Definition  of  Terms 

2.1.1  Scattering 

When  an  electromagnetic  wave  strikes  a  material  discontinuity  (in  the  case  of  radar, 
a  wave  traveling  in  free  space  strikes  an  object),  electric  currents  are  induced  in  the 
region  of  the  discontinuity.  These  currents  travel  within  or  on  the  surface  of  the  body  and 
re-radiate  an  electromagnetic  field.  The  radiation  can,  in  general,  be  in  any  direction. 
This  re-radiation  is  called  scattering.  The  manner  in  which  an  object  scatters  radiation  is 
highly  dependent  on  the  geometry  of  an  object  and  the  direction  of  the  incident  radiation, 
and  thus  any  given  object  has  its  own  unique  scattering  signature. 

2.1.2  Radar  cross  section 

Radar  cross-section  (RCS)  is  a  measure  of  the  magnitude  of  the  electromagnetic 
scattering  of  an  object  with  reference  to  a  metal  sphere.  If  a  metal  sphere  were 
constructed  to  produce  the  same  scattered  power  that  the  object  scatters  the  cross- 
sectional  area  of  that  sphere  would  be  the  object’s  radar  cross-section.  RCS  is  defined 
mathematically  as: 

IE  f  (1) 

<j  -  lim  4  7rR2  | — s-~ 

|E,| 

Where  R  is  the  range  from  the  transmitting  antenna  to  the  scatterer,  Es  is  the  scattered  E- 
field  (V/m),  and  £,  is  the  incident  E-field. 
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2.1.3  Bistatic  scattering  regions 


At  different  bistatic  angles,  measurement  considerations  can  change  significantly. 
For  the  purposes  of  this  project,  the  regions  of  bistatic  scattering  are  separated  into  three 
zones:  quasi-monostatic,  normal  bistatic,  and  forward  scatter.  Figure  2  shows  these 
regions  and  the  placement  of  the  receiving  antenna  at  the  boundaries. 


Rx 


Tx  — < 
Rx  ^ 


Rx 
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Figure  2:  Bistatic  Scattering  Regions 


In  the  quasi-monostatic  region,  the  transmitter  and  receiver  are  not  co-located  as 
in  the  true  monostatic  case,  but  have  a  small  angular  separation.  The  size  of  the  quasi¬ 
monostatic  region  is  dependent  on  the  type  of  object  to  be  measured.  For  many  simple 
objects,  the  quasi-monostatic  RCS  will  be  almost  identical  to  the  monostatic  RCS,  with  a 
correction  applied.  Corrections  of  this  type,  called  monostatic-to-bistatic  equivalence 
theorems  (MBETs),  have  been  studied  in  depth  in  (Kell,  1965)  and  will  not  be  elaborated 
on  any  further  in  this  paper.  More  complex  objects  may  have  an  RCS  that  varies 
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significantly  at  small  angles,  and  the  region  defined  to  be  quasi-monostatic  may  be  much 
smaller. 

In  the  normal  bistatic  region,  MBETs  in  general  do  not  apply.  RCS  in  this  region  is 
generally  less  than  in  the  quasi-monostatic  region,  but  large  variances  in  RCS  may  arise 
as  a  result  of  the  object  geometry.  As  a  consequence,  reference  object  selection  for 
bistatic  calibration  in  this  region  is  very  difficult. 

In  the  forward  scattering  region,  the  electromagnetic  signature  of  the  object  is  better 
characterized  by  it’s  shadow ,  or  absence  of  scattered  energy  in  this  area.  However,  this  is 
exceedingly  difficult  to  measure.  The  power  received  at  the  receiver  begins  to  be 
dominated  by  the  incident  field  rather  than  the  scattered  field,  as  the  receiver  approaches 
the  main  beam  of  the  transmitter.  In  this  region,  accurate  measurements  are  nearly 
impossible  due  to  the  large  ratio  of  incident  energy  to  scattered  energy. 

2.1.4  Scattering  mechanisms 

The  type  of  scattering  mechanism  that  dominates  on  a  particular  body  is  an 
important  factor  in  reference  object  selection.  It  is  desired  that  the  bodies  scatter  in  a 
predictable  way  for  many  bistatic  angles,  and  some  mechanisms  are  more  predictable 
than  others. 

Specular  reflection  is  the  dominant  scattering  mechanism  for  bodies  with  features 
much  larger  than  one  wavelength  (X).  It  occurs  when  a  wave  strikes  a  smooth  surface 
and  reflects  back  at  the  same  angle  from  the  normal  vector  to  the  surface  (Snell’s  Law). 
The  method  of  Geometrical  Optics  (GO)  predicts  this  mechanism  by  the  use  of  ray 
tracing  (Balanis,  1989).  Specular  reflection  tends  to  produce  a  larger  radar  return  from  an 
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incident  plane  wave,  because  the  wave  spreading  is  much  less  than  in  the  case  of 
diffraction,  dependent  on  the  geometry.  Spreading  is  a  term  used  to  describe  the 
attenuation  of  an  electromagnetic  wave,  due  to  the  divergence  of  rays  as  it  propagates. 
For  a  plane  wave,  the  rays  spread  at  a  rate  proportional  to  the  square  root  of  the  distance. 
As  the  wave  reflects  off  an  object,  the  spreading  in  the  return  path  will  be  a  function  of 
the  curvature  of  the  illuminated  surface.  A  convex  surface  with  a  small  radius  of 
curvature  will  cause  the  wave  to  spread  more,  and  a  flat  surface  will  introduce  no 
additional  spreading.  Thus,  for  objects  with  relatively  large,  flat  features,  specular 
reflection  will  be  the  dominant  scattering  mechanism. 

Diffraction  can  be  produced  in  a  number  of  different  ways.  For  complex 
geometries,  it  is  likely  that  diffraction  may  be  the  dominant  scattering  mechanism  on  the 
body.  For  a  large  set  of  geometries,  diffraction  can  be  modeled  as  a  local  phenomenon 
and  solved  for  numerically  with  the  Uniform  Geometrical  Theory  of  Diffraction  (UTD) 
(Kouyoumjian,  Pathak,  1974).  This  technique  predicts  diffraction  from  features  that 
create  an  instantaneous  change  of  the  radius  of  curvature  on  an  object.  It  models  an 
arbitrary  object  as  an  assembly  of  canonical  objects  that  the  code  can  calculate  precisely. 
This  set  includes  objects  such  as  plates,  cylinders,  spheroids,  ellipsoids,  cone  frustrums, 
etc.  Other  mechanisms  of  diffraction  are  surface  traveling  waves ,  which,  depending  on 
the  polarization  of  the  incident  wave,  will  propagate  along  the  surface  of  the  body  and 
produce  diffractions  at  edges.  Surface  waves  that  travel  in  the  non-illuminated  region  (or 
shadow  region)  of  the  object  are  called  creeping  waves. 

An  object  with  a  large  number  of  diffraction-producing  features  can  be  calculated 
using  asymptotic  high-frequency  codes.  Codes  utilizing  GO,  Physical  Optics  (PO),  UTD, 
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or  the  Physical  Theory  of  Diffraction  (PTD),  are  often  used,  but  it  is  sometimes  not  well- 
known  whether  the  solution  that  is  output  will  be  completely  accurate,  especially  in  the 
case  of  bistatic  predictions.  For  complex  bodies  it  is  usually  desirable  to  use  a  prediction 
code  based  on  the  Method  of  Moments  (MoM)  (Harrington,  1993),  which  can  produce  an 
exact  solution  at  the  cost  of  computational  time. 

2.1.5  Polarization 

The  term  polarization  refers  to  the  alignment  of  the  electric  field  vector  with  respect 
to  whatever  coordinate  system  is  in  use.  In  general,  the  electric  field  will  rotate  as  it 
propagates  through  the  medium.  This  is  elliptical  polarization.  In  the  case  of  most  radar 
measurements,  the  linearly  polarized  waves  are  used;  meaning  the  alignment  of  the 
electric  field  vector  does  not  change  as  a  function  of  time  or  distance. 

The  E-field,  H-field,  and  direction  of  propagation  (P)  are  all  mutually  orthogonal,  as 
seen  in  Figure  3. 

z 


Figure  3:  Wave  Polarization 

The  typical  naming  convention  for  polarization  is  horizontal  (H-pol)  and  vertical 
(V-pol).  An  H-polarized  wave  is  transmitted  with  the  E-field  vector  aligned  parallel  to 
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the  horizontal  axis  (in  the  x-y  plane),  and  a  V-polarized  wave  is  transmitted  with  the  E- 
field  vector  aligned  perpendicular  to  the  horizontal  axis  (parallel  to  the  z-axis).  A 
polarimetric  measurement  is  one  that  transmits  and  receives  every  combination  of 
polarizations.  In  other  words,  HH,  VV,  (co-polarized);  HV,  and  VH  (cross-polarized) 
data  are  collected,  where  the  first  letter  identifies  the  received  polarization,  and  the 
second  letter  identifies  the  transmitted  polarization.  Thus  the  polarization  “HV”  means 
that  a  vertically  polarized  field  was  transmitted,  and  the  receiver  collected  the 
horizontally  polarized  field.  The  collection  of  cross-polarized  data  gives  a  more  unique 
description  of  the  object  scattering  characteristics.  Some  object  geometries  may 
depolarize  one  or  more  polarizations,  meaning  the  polarization  that  is  reflected  had  a 
component  perpendicular  to  the  polarization  that  was  transmitted.  The  field  of  radar 
polarimetry  is  dedicated  to  the  analysis  of  this  type  of  measurement. 

2.1.6  Scattering  matrix 

A  polarimetric  measurement  gathers  four  complex  RCS  values  for  a  measurement 
of  a  single  bistatic  angle  and  frequency.  These  values  are  arranged  into  a  2  x  2  complex 
scattering  matrix.  Scattering  matrices  usually  take  the  form 

[HH  VH~\ 

S  = 

HV  VV 

where  the  elements  HH,  VH,  HV,  and  VV  are  the  amplitude  and  phase  of  each  respective 
measurement.  The  scattering  matrix  is  also  referred  to  as  the  Sinclair  matrix.  The  matrix 
representation  of  a  polarimetric  measurement  becomes  very  useful  in  performing 
calibrations  for  analysis  and  calibration,  as  we  will  see  in  the  chapters  to  follow. 
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2.1.7  Subsystem  Distortion  matrix 


The  errors  introduced  into  a  measurement  as  a  result  of  non-ideal  hardware  systems 
can  be  calculated  using  two  2x2  complex  matrices  called  subsystem  distortion  matrices. 
The  two  matrices  correspond  to  the  distortion  of  the  receiver  subsystem  and  the 
transmitter  subsystem,  as  given  by: 


R  = 

Rhh 

RVh 

T  = 

~Thh 

I - 

ll 

Rhv 

RVv  _ 

Jhv 

Tyy  J 

Where  R  is  the  receiver  subsystem  distortion  matrix  and  T  is  the  transmitter  subsystem 
distortion  matrix.  The  matrices  operate  on  the  actual  value  of  the  object  scattering  matrix 
(thereby  distorting  the  actual  RCS)  to  yield  the  value  of  the  measurement  of  the  same 
object.  The  object’s  measured  scattering  matrix  can  be  expressed  as 

M  =  RST  (2) 

Where  M  is  the  measured  2x2  scattering  matrix  and  S  is  the  actual  scattering  matrix. 

2.1.8  Near-field  and  far -field 

In  RCS  measurements,  the  incident  wave  is  assumed  to  be  a  plane  wave.  The 
reason  is  that  in  an  actual  radar  scenario,  especially  that  involving  ground-based  radar 
illuminating  an  aircraft,  the  object  is  located  at  a  long  range  from  the  illuminating 
antenna.  The  wavefront  of  this  antenna  is  roughly  spherical,  so  the  radius  of  curvature  of 
the  wavefront  is  the  range  from  the  antenna  to  the  object.  At  long  ranges  such  as  these, 
the  radius  of  curvature  is  large  enough  so  that  the  wavefront  looks  planar. 

The  region  in  which  a  spherical  wave  front  can  be  assumed  to  be  planar  is  called  the 
far-field.  At  ranges  smaller  than  this,  it  is  called  the  near-field. 


12 


2.2  Coordinate  Convention 


With  the  wide  variety  of  measurement  facilities  and  conventions,  the  H-  and  V- 
polarization  conventions  can  sometimes  be  ambiguous.  For  this  reason  it  is  desirable  to 
link  the  convenience  of  the  H  and  V  notation  with  the  spherical  coordinate  system  in 
order  to  eliminate  confusion.  For  the  measurement  geometry  of  the  EMSL,  the  plane  in 
which  the  receiver  tracks  is  taken  to  be  the  y-z  plane,  and  the  z-axis  is  perpendicular  to 
the  ground.  The  coordinate  convention  most  convenient  for  this  configuration  is 
described  in  (Ulaby,  Elachi,  1990),  and  is  depicted  in  Figure  4  below,  although  the  same 
system  is  used  for  measurements  in  other  configurations  (for  example,  in  (Sarabandi,  et 
al,  1990)). 


Figure  4:  Forward  Scattering  Coordinate  Convention  (Ulaby  and  Elachi,  1990) 

Here  the  h  and  v  unit  vectors  corresponding  to  their  respective  polarizations 

correspond  directly  to  the  <|>  and  0  unit  vectors,  respectively.  The  v  and  h  directions  for 
the  incident  field  (vj  and  hi)  and  the  scattered  field  (vs  and  hs)  are  defined  with  respect  to 
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the  direction  of  propagation  (kj  and  ks),  making  this  a  wave  oriented  coordinate  system. 
The  set  of  equations  governing  the  directions  of  propagation  and  polarization  in  relation 
to  the  spherical  coordinate  system  are: 


kj  =  cos  fa  sin  6  ix  +  sin$  sin  6 fa  +  cos  Op 


zxkj 

h;  = - —  -  -sin^.x  +  cos^.y 

zxkj 


Vj  =  hj  x  kj  =  cos  fa  cos  Gfa  +  sin  fa  cos  6 fa  -  sin  Gfa 


ks  =  cos  fa  sin^x  +  sin^  sin^y  +  cosGfa 


zxks 

hs  =  =  -sin^x  +  cos^y 

zxk$ 


vs  =  hs  x  ks  =  cos^  cos  6sx  +  sin  fas  cos Gs y  -  sin  Gs 


(3) 

(4) 

(5) 

(6) 

(7) 

(8) 


For  bistatic  calibration,  the  primary  purpose  of  such  a  coordinate  system  is  to 
eliminate  the  ambiguity  in  calling  a  wave  H-  or  V-polarized. 


2.3  Sources  of  Measurement  Inaccuracy 

The  motivation  for  radar  calibration  is  to  minimize  the  effect  of  systemic  errors  on 
the  measurement.  When  raw,  uncalibrated  data  is  taken,  it  includes  undesirable 
information  that  must  be  filtered  out  before  the  data  can  have  any  value.  There  are 
several  sources  of  this  measurement  inaccuracy. 
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2.3. 1  Reflections  from  chamber  walls  and  mounting  apparatus 

An  indoor  anechoic  chamber  is  covered  with  radar-absorbing  material  (RAM)  to 
minimize  reflections  from  the  transmitting  antenna  that  would  interfere  with  the  scattered 
field  from  the  object.  RAM  typically  covers  the  walls  of  the  chamber  as  well  as  the 
object-mounting  pedestal.  This  usually  reduces  the  magnitude  of  the  reflection  down  to  a 
reasonable  level,  but  it  is  often  desirable  to  reduce  it  further  by  means  of  background 
subtraction  and  range  gating,  but  these  methods  introduce  errors  of  their  own.  These 
methods  are  discussed  in  the  next  section. 


2.3.2  Antenna  coupling 

At  large  bistatic  angles,  as  the  receiving  antenna  moves  into  the  main  beam  of  the 
transmitter,  the  incident  field  of  the  transmitter  becomes  much  larger  than  the  scattered 
field  from  the  object,  and  it  is  difficult  to  subtract  the  two  (Figure  5).  There  are  two 
techniques  for  mitigating  this  error,  discussed  in  the  next  section. 
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Incident  Field,  Ej 

Figure  5:  Antenna  Coupling  in  Forward  Scattering  Region 
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2.3.3  Subsystem  distortion 


Subsystem  distortion,  or  cross-polarization  impurities  occurring  in  the  transmitter  or 
receiver  subsystems,  can  have  a  subtle  yet  significant  affect  on  the  accuracy  of  a 
polarimetric  measurement,  especially  when  considering  depolarizing  objects,  or  objects 
which  scatter  both  polarizations  with  a  single,  linearly  polarized  incident  wave.  Due  to 
hardware  imperfections  and/or  the  positioning  of  either  antenna,  antenna  cross¬ 
polarization  leakage  is  inevitable.  This  leakage  can  occur  in  either  the  transmitter  or 
receiver,  and  often  both.  When  a  wave  of  a  particular  polarization  is  transmitted,  the 
transmitted  wave  will  have  a  small  component  of  the  opposite  polarization.  In  the 
receiver,  a  particular  antenna  oriented  to  receive  a  single  polarization  may  have  a  rotation 
that  lets  it  receive  some  of  the  opposite  polarization  simultaneously. 

2. 3. 4  Object  alignment  error 

It  can  sometimes  be  very  difficult  to  precisely  align  an  object  for  measurement  or 
calibration.  A  misalignment  of  an  object  will  always  produce  a  deviation  from  the 
theoretical  behavior  of  the  object,  but  some  objects  may  be  more  sensitive  to  alignment 
error,  and  produce  deviations  that  are  much  worse  than  another  object.  If  the  object  is 
used  for  calibration,  much  care  needs  to  be  taken  to  choose  an  object  that  is  not  sensitive 
to  alignment  error,  or  a  precise  and  sometimes  complex  alignment  system  is  necessary. 

A  sphere,  for  instance,  is  completely  insensitive  to  alignment  error,  because  it  has 
rotational  symmetry  about  any  axis.  This  makes  it  a  popular  calibration  object.  The 
radar  return  from  a  sphere,  however,  is  relatively  low,  and  the  location  of  the  noise  floor 
of  the  measurements  may  preclude  its  use.  Louder  objects  such  as  dihedrals  and 
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trihedrals  are  often  used  in  monostatic  measurements,  but  these  objects  have  a  degree  of 
alignment  sensitivity  that  may  introduce  error  into  the  measurement. 

2.3.5  Near-field  effects 

Scattering  behavior  for  an  object  is  not  the  same  in  the  near-field  as  it  is  in  the  far- 
field.  A  spherical  wavefront  does  not  illuminate  a  distributed  object  with  a  single  phase. 
Scattering  points  along  the  object  are  illuminated  at  a  phase  dictated  by  their  location 
relative  to  the  axis  defined  by  the  direction  of  propagation.  The  object  will  scatter  these 
various  phases,  and  they  will  interfere,  as  in  the  far-field  case;  but  this  interference  will 
obviously  produce  different  patterns  because  the  relative  phase  of  each  scattering  point 
from  a  near-field  excitation  are  different  than  that  of  a  far-field  excitation. 

Simulating  a  planar  wave  front  in  an  anechoic  chamber  is  a  major  problem.  Indoor 
ranges  are  limited  in  size,  so  parabolic  reflectors  are  often  used  to  distort  the  wave  front 
to  make  it  planar.  Measurement  chambers  that  use  parabolic  reflectors  to  create  plane 
waves  are  called  compact  ranges. 

2.4  Range  Gating  and  Background  Subtraction 

Background  subtraction  is  done  by  performing  and  empty  room  measurement:  that 
is,  a  measurement  with  the  object  removed  and  the  pedestal  and  antennas  in  exactly  the 
same  position.  The  empty  room  return  is  coherently  subtracted  from  the  measurement  in 
order  to  isolate  the  object  from  its  environment. 

Two  phenomena  are  often  present  in  the  background  of  RCS  measurements  that 
cannot  be  subtracted:  first,  interaction  between  the  object  and  its  pedestal,  most  often 
involving  multipath  reflections  between  the  portion  of  the  object  facing  the  pedestal  and 


17 


the  more  specular  (flat)  portions  of  the  pedestal.  Since  this  interaction  is  not  possible 
with  the  object  absent,  it  cannot  be  subtracted.  Second,  shadowing  as  a  result  of  the 
object  blocking  the  antenna’s  view  of  the  pedestal  can  also  be  a  problem.  For  objects 
with  a  large  footprint  on  the  pedestal,  in  the  empty  room  measurement,  the  antenna  sees  a 
portion  of  the  pedestal  not  seen  in  the  measurement,  which  turns  up  in  the  time  domain  as 
a  spurious  return  overlaid  upon  the  object.  Often  this  shadowing  cannot  be  removed  by 
range  gating  since  the  shadowed  surface  is  many  times  very  close  to  the  object. 

The  next  set  of  time-domain  plots  help  to  visualize  the  effects  of  the 
aforementioned  error  sources.  Figure  6  is  the  calibrated  RCS  of  a  metallic  sphere  with 
and  without  background  subtraction. 

The  vertical  axis  is  the  bistatic  angle  (4  -  184  degrees),  the  horizontal  is  the  range 
(0  ns  is  the  center  of  the  chamber),  and  the  color  dimension  is  the  RCS  in  dBsm.  The 
feature  labeled  (1.)  is  the  coupling  between  the  Tx  and  Rx  antennas.  One  can  gain 
intuition  into  the  reason  for  the  arc-like  shape  of  the  coupling  by  visualizing  the  chamber 
with  the  Rx  antenna  in  the  forward  scattering  position.  At  the  edge  of  the  forward  scatter 
region  (around  150  degrees),  the  angular  position  of  the  receiver  is  offset  from  the 
equiphase  plane  of  the  incident  wave.  This  means  that  the  coupling  that  strikes  the 
antenna  at  this  point  is  from  a  more  recent  phase  front,  and  the  apparent  location  of  the 
coupling  at  this  angle  would  be  closer  to  the  Tx  antenna.  As  the  bistatic  angle  increases, 
the  apparent  location  of  the  coupling  grows  closer  to  the  pedestal,  as  the  plane  of  the 
aperture  of  the  Rx  antenna  moves  perpendicular  to  the  incident  phase  front.  After 
background  subtraction,  feature  1  is  greatly  suppressed,  but  is  still  visible  in  the  time- 
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domain  plot.  However,  the  level  is  around  20-30  dBsm  lower  than  the  return  from  the 
sphere. 


184  deg 
2f».0  ~ 


■6.7 


-33.3 


-60.* 


4  deg 

Figure  6:  Time-Domain  Sphere  RCS  (HH-pol)  With  Background  Subtraction 
(right)  and  Without  Subtraction  (left) 

The  red  vertical  line  labeled  feature  2  is  the  RCS  of  the  sphere  itself.  The  line  shifts 
somewhat  as  a  function  of  bistatic  angle,  which  is  to  be  expected.  In  the  backscatter 
direction  (4  degrees)  the  specular  point  occurs  at  about  -1  ns — closer  to  the  Tx  antenna. 
This  is  because  the  sphere  is  centered  at  0  ns,  and  the  specular  point  occurs  at  one  radius- 


length  behind  the  center.  As  the  bistatic  angle  increases,  the  specular  return  moves  as  the 
specular  point  moves  to  the  top  of  the  sphere,  and  the  creeping  wave  grows  larger  and 
interferes  with  the  specular  return.  As  expected,  the  forward  scatter  RCS  of  the  sphere  is 
quite  large.  This  is  not  the  energy  radiated  by  the  sphere  in  the  forward  direction,  but  the 
absence  of  radiation,  or  shadow  of  the  sphere. 


The  phenomenon  of  shadowing  can  be  seen  in  feature  3  of  Figure  6.  It  is  clear  that 
this  return  is  not  part  of  the  object  return,  since  the  return  of  the  object  is  largely  localized 
to  the  region  between  —1  and  0  ns.  At  near-monostatic  angles,  the  receiver  sees  the  same 
pedestal  in  the  empty  room  measurement  as  it  does  with  the  object  in  place.  As  the 
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receiver  moves  past  90  degrees,  the  absence  of  the  sphere  reveals  a  portion  of  the 
pedestal  that  was  shadowed;  in  this  case  it  is  the  flat  horizontal  surface  and  the  concave 
“cup”  that  the  sphere  rests  on. 

Feature  4  is  a  reflection  within  the  transmission  lines  of  the  system.  The  return 
from  the  sphere  enters  the  receiver,  propagates  along  the  coaxial  cables  to  the  network 
analyzer,  and  due  to  imperfect  impedance  matching  in  the  system,  it  is  attenuated, 
bounces  back  to  the  antenna,  and  bounces  again  back  to  the  network  analyzer  and  appears 
as  a  delayed  and  attenuated  replica  of  the  object. 

Feature  5  is  a  by-product  of  the  time-domain  transform.  To  get  a  time-domain  plot, 
an  FFT  is  performed  over  the  bandwidth  of  the  measurement.  Since  the  bandwidth  is 
finite,  a  rippling  occurs  as  a  result  of  the  transform,  which  gives  rise  to  feature  5.  This 
rippling  occurs  everywhere  in  the  plot,  but  is  most  apparent  at  184  degrees  because  the 
return  is  much  larger  at  this  angle. 

Figure  7  shows  the  same  calibrated  and  subtracted  data  after  a  range  gate  is  applied 
with  a  width  of  10  ns.  In  the  vertical  regions  around  the  main  object  return  the  rippling 
due  to  the  range  gate  is  visible.  Comparing  this  plot  with  Figure  6,  one  can  see  that  the 
object/pedestal  interaction  is  still  present.  It  is  possible  and  quite  easy  to  make  the  range 
gate  tight  enough  to  eliminate  this  interaction,  and  still  leave  all  of  the  object  return,  but 
the  tighter  range  gate  will  make  the  rippling  effects  more  prevalent.  As  it  is,  the  data 
taken  at  the  edge  of  the  band  (here,  the  data  around  6  GHz  and  14  GHz)  is  already 
unusable  due  to  the  range  gate.  If  it  is  necessary  to  use  these  frequencies  for  analysis  or 
for  calibration,  measurements  must  be  done  on  a  wider  bandwidth. 
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Figure  7:  Time-Domain  Sphere  RCS  (HH-pol)  with  10  ns  Range  Gate  Applied 

2.5  Types  and  Applicability  of  Calibration 

For  a  given  measurement,  there  is  more  than  one  type  of  calibration  that  can  be 
used.  The  determination  of  the  calibration  type  will  be  dependent  on  the  quality  of  the 
measurement  facility,  the  measurement  conditions,  and  the  degree  of  accuracy  that  is 
required.  Radar  calibration  can  be  of  three  basic  types: 

•  Amplitude  and  phase  (Type- 1 ) 

•  Simple  polarimetric  (Type-2) 

•  Full  Polarimetric  (Type-3) 

2.5.1  Amplitude  and  phase  (Type-1)  calibration 

A  Type-1  calibration  applies  a  complex  constant  scale  factor  to  each  scattering 
matrix.  The  scale  factor  is  usually  calculated  by  measuring  an  object  with  a  known  RCS 
and  phase  center,  and  comparing  the  ratio  of  the  measured  and  exact  responses.  Type-1 
calibrations  all  have  a  common  form.  For  each  polarization  channel  (HH,  VH,  HV,  and 
VV),  the  calibrated  RCS  is: 
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cal  _  tar  theoretical 


(9) 


where  (7cai  is  the  calibrated  response  of  the  object.  The  calibrated  scattering  matrix  takes 
the  form  of: 


G cat,. m 

G  cal, VH 

° cal,. HV 

G cal, VV 

(10) 


where  Scai  is  the  calibrated  scattering  matrix  of  the  unknown  object  and  the  cfs 
correspond  to  the  amplitudes  an  phases  of  their  corresponding  polarizations. 

The  advantage  of  using  a  Type-1  calibration  technique  is  in  its  simplicity.  Only  one 
calibration  reference  object  is  necessary  and  it  is  computationally  very  simple  to 
implement.  The  disadvantages  of  using  only  a  Type-1  technique,  however,  are 
numerous.  An  object  must  be  selected  that  has  high  RCS  for  every  polarization,  in  order 
for  the  calibration  to  be  accurate  for  the  complete  scattering  matrix.  Often,  however,  a 
Type-1  technique  is  used  only  to  calibrate  co-polar  channels.  It  is  practical  to  perform  a 
Type-1  calibration  with  more  than  one  reference  object,  but  at  this  point  the  technique  is 
not  much  more  efficient  than  a  Type-2  or  3  technique.  In  addition,  a  Type-1  technique 
does  not  account  for  distortion  occurring  in  the  transmitter  and  receiver  subsystems. 


2.5.2  Simple  polarimetric  (Type-2)  calibration 

In  order  to  account  for  system  distortion,  a  Type-2  or  3  calibration  technique  is 
usually  employed.  The  Type-2  calibration  technique,  or  simple  polarimetric  calibration, 
accounts  for  zero-order  antennae  polarization  distortion  effects.  Each  of  the  co-polarized 
subsystem  distortion  terms,  Rhh,  Thh,  Rvv,  Tvv,  cannot  be  determined  distinctly,  but  can 
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be  found  as  products  of  the  receiver  subsystem  distortion  (R)  and  transmitter  subsystem 
distortion  (T)  for  each  polarization  channel  combination.  The  distortion  terms  are 
expressed  as  a  4x4  complete  system  distortion  matrix,  C  is  introduced,  which  carries  the 


same  information  as  the  2x2  matrices  R  and  T.  The  C  matrix  has  the  form: 


RhvThh 

RhhTvh 

& 

a* 

RyH^HH 

Rw^hh 

RvhTvh 

Rw'Rvh 

(Mortensen,  1995) 

Rhv^hv 

RhhTw 

RHVTyv 

RyV^HH 

Rw^hv 

RvhTw 

i 

d 

where  the  elements  R**  and  T**  correspond  directly  to  the  elements  of  the  R  and  T 
matrices  mentioned  earlier.  For  a  simple  polarimetric  calibration,  only  the  elements  lying 
on  the  diagonal  of  C  (analogous  to  the  products  of  the  elements  lying  on  the  diagonals  of 
R  and  T)  are  calculated.  Thus  the  form  of  the  C  matrix  for  a  simple  polarimetric 
calibration  is: 


c  = 


RhH^HH 
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0 
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RyV^HH 

0 
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0 
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^HH^VV 


0 

0 

0 

RvyTw 


(11) 


All  other  terms  correspond  to  cross-polarization  impurity  of  the  subsystems,  the 
calculation  of  which  requires  a  full-polarimetric  calibration  (Type-3).  It  is  important  to 
note  that,  though  a  Type-2  calibration  solves  for  the  products  of  the  R  and  T  coefficients, 
this  calibration  can  give  no  information  about  R  and  T  separately,  so  the  combined 
matrix  (C)  must  be  used  instead. 

A  Type-2  calibration  will  requires  a  minimum  of  two  calibration  objects:  a  single 
reference  calibration  object,  for  which  the  theoretical  solution  is  precisely  known,  and  an 
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additional  depolarizing  object  for  cross-polarization  correction,  also  with  a  known 
scattering  matrix. 

2.5.3  Fully  polarimetric  (Type-3)  calibration 

A  fully  polarimetric  calibration  is  able  to  solve  for  all  eight  coefficients  of  the  R 
and  T  matrices  (and  thereby  all  16  terms  of  the  C  matrix).  The  first  and  second  order 
subsystem  distortion  effects  are  all  accounted  for,  including  the  cross-polarization 
impurity  of  the  individual  transmit  and  receive  antennas. 

The  Type-3  technique  is  the  best  available  for  accuracy  in  calibration,  at  the  cost  of 
complexity.  A  fully  polarimetric  technique  requires  three  calibration  objects  whose 
theoretical  scattering  matrices  are  precisely  known  (Whitt,  et  al,  1991).  By  definition,  a 
fully  polarimetric  technique  must  solve  for  the  coefficients  of  the  distortion  matrices 
separately,  as  opposed  to  a  Type-2  technique.  This  neglects  the  dependent  nature  of  each 
distortion  term  on  another.  A  rigorous  mathematical  process  is  necessary  in  order  to 
solve  for  the  subsystem  distortion  matrices  simultaneously.  Once  the  distortion  matrices 
are  computed,  the  corrected  scattering  matrix  for  an  unknown  object  is  given  by: 

S  =  R-IMT-1  (12) 

Where  S  is  the  corrected  scattering  matrix  and  M  is  the  measured  matrix,  with 
background  subtraction  performed.  This  relationship  can  be  derived  from  equation  (2)  by 
re-arranging  terms  under  the  assumption  that  the  subsystem  distortion  matrices  are 
invertible  (Ulaby  and  Elachi,  1990).  A  technique  of  this  type  requires  a  minimum  of 
three  calibration  reference  objects  with  known  scattering  matrices. 
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3  BISTATIC  CALIBRATION  METHODS 


Several  existing  calibration  methods  were  sampled,  for  the  reasons  listed  above,  and 
they  are  listed  in  this  section. 

3.1  Basic  Type-1  Calibration 

This  technique  is  not  documented  in  the  literature,  but  it  is  essentially  the  only  way 
one  would  perform  a  Type-1  calibration.  An  object  is  selected  which  has  zero  or 
negligible  cross-polarization  return  (a  sphere  or  circular  disk,  for  example).  The 
calibrated  measurement  values  for  each  polarization  (HH  and  W)  are  computed  by  Eqn 

(9). 


3.2  EMSL  Simple  Polarimetric  Calibration 

This  is  the  bistatic  calibration  method  (Type-2)  currently  used  in  the  EMSL 
(Mortensen,  1995).  A  circular  disk  is  used  for  the  reference  object.  This  object  provides 
a  high  specular  RCS  and  well-defined  phase  center,  as  well  as  an  accurate  theoretical 
prediction  by  Physical  Optics,  although  the  disk  is  somewhat  difficult  to  align  in  the 
quasi-monostatic  configuration.  The  reference  object  is  used  to  solve  for  the  (1,1)  and 
(4,4)  elements  of  the  C  matrix,  or  the  coefficients  RhhThh  and  RwTw.  The  (2,2)  and 
(3,3)  elements  of  the  matrix  are  computed  first  by  correcting  the  amplitudes  of  the 
depolarizing  object  with  the  (1,1)  and  (4,4)  coefficients  already  found,  and  computing  the 
ratio  of  the  measured  cross-polar  amplitude  to  the  absolute  corrected  amplitude.  The  PO 
approximation  for  the  backscatter  RCS  of  the  disk,  assuming  a  perfectly  conducting 
surface  is  (in  dB), 
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(13) 
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where  r  is  the  radius  of  the  sphere  in  meters,  and  X  is  the  wavelength  in  meters.  For 
matrix  multiplication  by  a  4x4  matrix,  the  normal  2x2  scattering  matrix  format  is 
modified  to  a  4x1  vector  of  the  form: 


^  disk 


s 


disk 


0 

0 


rr  W 
r disk 


(14) 


The  elements  Ci,i  and  C4,4  are  solved  by  relating  the  measured  scattering  matrix  of 
the  disk  (Mdisk)  to  the  theoretical  prediction: 


m”,"  =  cus™k  +  c 


,2  ctVH  ,  ^1,3  nHV  ,  iol,4  nW 
° disk  +  C  +  U  °disk 


(15) 
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(16) 


For  the  circular  disk,  the  cross-polarized  elements  are  approximately  zero,  and  we 
assume  that  the  elements  due  to  the  cross-talk  in  the  system,  C1’4  and  C4’1  ,  are  negligible 
so  the  equation  simplifies  to: 
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(18) 


We  can  now  rearrange  this  formula  to  solve  for  C1,1  and  C4’4  and  use  these 
coefficients  to  find  the  corrected  amplitude  of  the  depolarizing  object  This  allows  us  to 
use  the  wire  mesh  as  a  calibration  object  without  knowing  its  exact  RCS.  The  EMSL 
uses  a  mesh  of  parallel  wires  as  a  depolarizer.  Two  measurements  are  made  on  the  mesh: 
one  in  the  vertical  orientation,  and  one  with  the  mesh  tilted  45  degrees  from  vertical  in 
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the  counter-clockwise  direction.  In  the  vertical  orientation,  it  is  reasonable  to  assume  that 
only  the  VV-polarized  wave  will  be  reflected,  and  all  others  will  be  transmitted  through 
the  mesh,  giving  the  4x1  scattering  vector  the  form 
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In  the  tilted  configuration,  the  vector  takes  the  form 
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where  0  is  the  bistatic  angle.  This  approximation  is  valid  for  small  bistatic  angles  (Bicci, 

1992).  In  the  tilted  orientation,  the  amplitude  of  the  monostatic  (0  =  0)  RCS  of  the  mesh 

should  be  the  same  for  all  polarizations.  The  cosine  dependence  of  the  VH,  HV,  and  VV 
polarizations  can  be  seen  intuitively  by  visualizing  the  measurement  geometry.  In  the 
EMSL,  the  receiver  tracks  in  the  vertical  direction  and  therefore  over  the  object  (the  tilted 
wire  mesh,  for  instance)  rather  than  around  it  in  azimuth.  As  the  receiver  moves,  the  V- 
polarized  component  of  the  wave  that  is  transmitted  is  not  completely  V-polarized  from 
the  receiver’s  perspective.  This  polarization  distortion  at  the  receiver  is  expressed 
mathematically  as  a  loss  in  the  V-polarized  RCS  as  a  function  of  the  bistatic  angle. 

To  find  the  corrected  values  of  the  mesh  in  the  vertical  orientation,  we  simply 
divide  the  VV-component  (the  only  component  that  is  non-zero)  of  the  measurement  by 
the  element  of  C  corresponding  to  VV-pol,  C4’4: 
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From  the  corrected  magnitude  of  Smesh,vert  ,  we  can  find  the  corrected  monostatic 
RCS  of  the  vertical  mesh,  Smesh: 


^mesh  ~  ' 


cos2(-) 


Substituting  equation  (22)  into  equation  (20),  and  selecting  the  HV  and  VH 
polarizations,  we  get 
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And  this  is  all  the  information  we  need  to  solve  for  the  remaining  coefficients  of  the 
diagonal  C  matrix: 
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After  the  four  calibration  coefficients  have  been  computed,  the  calibration  can  be 
performed  for  a  single  frequency  by  a  simple  matrix  multiplication: 


S  =  C_1  (M-B) 


Where  S  is  the  calibrated  scattering  matrix,  M  is  the  measured  scattering  matrix. 


and  B  is  the  scattering  matrix  of  the  background. 
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3.3  Generalized  Dual-Antenna  Calibration 


This  technique,  proposed  in  (Whitt,  et  al,  1991)  and  elaborated  in  (Ulaby  and 
Elachi,  1990),  is  the  most  general  technique  available  in  current  literature.  Its  advantage 
is  that  it  makes  few  assumptions  about  the  form  of  the  scattering  matrix  of  each  object.  It 
requires  that  one  of  the  scattering  matrices  must  be  invertible,  and  assumes  that  the 
scattering  matrices  are  all  symmetric  (cross-polarized  components  are  equal).  This  is  a 
good  assumption  for  any  passive  object  in  the  monostatic  or  quasi-monostatic 
configurations,  but  is  not  always  true  for  bistatic  measurements,  in  general.  Therefore 
this  assumption  practically  limits  the  use  of  the  technique  to  dual-antenna  systems 
calibrated  in  the  quasi-monostatic  configuration. 

The  main  limitation  of  the  technique  is  the  requirement  that  the  theoretical  RCS  of 
each  of  the  three  reference  objects  be  exactly  known.  In  addition,  the  technique  is  much 
more  mathematically  complex  than  other  calibration  techniques,  but  this  obviously  has  no 
effect  on  the  ease  of  implementation  of  the  technique  once  the  mathematical  routine  is 
successfully  implemented. 

The  technique  is  able  to  account  for  all  eight  error  terms  included  in  the  matrix 
formula  of  equation  (2)  but  the  notation  includes  a  complex  scale  factor  (ce'^)  so  that  the 
relationship  between  the  measured  scattering  matrix  (M)  and  the  theoretical  matrix  (S) 
expands  to: 

M  =  c  e'*  RST  (28) 

The  technique  requires  three  reference  objects,  and  exact  theoretical  predictions  for 
each,  producing  three  similar  equations: 
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M.  =  ce"'RS,TJ  =  1,2,3 


(29) 


The  real  constant,  c,  will  be  the  same  for  each  object,  as  it  is  a  characteristic  of  the 
radar.  The  phase  constant  {(/>j)  is  subscripted  to  allow  for  the  phase  centers  of  each  object 
to  differ,  thus  eliminating  the  need  to  place  the  objects  such  that  their  phase  centers  all 
coincide.  We  must  solve  for  the  unknowns  R,  T,  c,  and  (fo .  Following  the  procedure  of 
Ulaby  and  Elachi,  the  steps  are  described  in  the  following  pages. 

We  first  assume  that  one  of  the  calibration  reference  objects  has  an  invertible 
scattering  matrix.  For  the  purposes  of  this  derivation  we  assume  that  matrix  to  be  the 
object  corresponding  toy  =  1.  We  first  define  the  similarity  transformations: 


Nx  =  e«A-*»>r,PTT 

(30) 

N  =  e,(*-wT-'pTT 

(31) 

Where  NT  =  M ;'M, ,  NT  =  M^M3 ,  PT  =  P,_,P2 ,  and  PT  =  P^'P, . 

It  is  now  necessary  to  express  the  above  matrices  Pt  and  Nj 

in  terms  of  their 

eigenvectors  and  eigenvalues. 
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the  eigenvalues  and  eigenvectors  are  related  by 

A;  =  Arem~*2) 

(34) 
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(35) 

30 


The  matrices  Xy  and  Yy  are  defined  as  2x2  matrices,  the  columns  of  which  are  the 
eigenvectors  of  Pt  and  Nj,  respectively.  and  AT ,  respectively,  are  diagonal  matrices 

of  the  eigenvalues.  Using  the  equation  (35)  and  manipulating  the  matrices  to  solve  for 
the  transmitter  subsystem  distortion  matrix,  T,  and  introducing  an  multiplicative  scale 
factor  to  account  for  the  arbitrary  constant  that  cannot  be  solved  for  in  the  above 
equations,  we  arrive  at  the  expression  of  equation  (36): 

T  =  XtCYt-'  (36) 


The  matrix  C  is  a  diagonal  matrix  of  real  constants  (ci  and  C2).  We  can  now  define  a  new 
set  of  equations, 


PTXT  =  Xy  At 

(37) 

NTYT  =  YyAT 

(38) 

and  then  say  that 


t=xtcy;’ 


(39) 


Again,  C  is  a  diagonal  matrix  of  constants  (c,  and  c2).  Understanding  that  the  T  of 
equation  (36)  and  of  equation  (39)  are  identical,  we  now  have  enough  information  to 
solve  for  the  arbitrary  constants  cj  and  C2  (as  well  as  c,  and  c2 ,  however  we  only  need 

one  set  of  constants).  Through  a  process  of  solving  the  above  matrix  equations  in  terms 
of  sets  of  scalar  equations,  we  find  the  constants  to  be: 


c\  =  A(YT)(*ny22 -—x12y2l)~' 
C\ 


(40) 
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(41) 


Q 

c2  =  A(Yt)(  -^11^22  -^12^21) 
c2 

where  the  scalars  x*»  and  y**  are  the  corresponding  elements  of  the  matrices  Xy  and  Yy. 
After  substituting  the  right  side  of  equations  (40)  and  (41)  into  the  matrix  C  in  equation 
(36),  we  now  have  a  unique  value  for  the  transmitter  subsystem  distortion  matrix. 

The  receiver  subsystem  distortion  matrix  is  solved  for  in  an  identical  matter.  First 
the  similarity  transformations  of  equations  (42)  and  (43)  are  created, 


Nr  =e'^-^T',PRT 

(42) 

Nr  =e'^“A)T“1PRT 

(43) 

Where  NR  =  M2M“' ,  NR  =  M3M^ ,  PR  =  P^,1 ,  and  PR  =  P3P, 1 . 

By  subsequently  expressing  the  N  and  P  matrices  of  (42)  and  (43)  in  terms  of 
eigenvectors  and  eigenvalues,  as  in  (37)  and  (38),  we  can  specify  the  matrix  R.  The 
constant  amplitude  and  phase  factors  included  in  equation  (29),  however,  remain 
unknown.  The  amplitude  factor  c  can  be  solved  for  by  substituting  the  known  R  and  T 
matrices  back  into  equation  (28),  and  solving  for  c : 

|(RP,T)„|  (44) 

The  matrix  Pj  should  be  selected  such  that  it  is  the  most  accurate  of  the  three, 
thereby  giving  a  more  accurate  estimation  of  c.  The  mnth  elements  of  the  matrices  should 
probably  be  selected  to  correspond  to  one  of  the  co-polar  channels  (HH  or  VV),  since 
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they  are  in  general  larger  in  magnitude,  and  therefore  are  less  affected  by  noise.  The 
absolute  phase  factor  (<f>r)  is  given  by: 


<!>'=  ^  +  tan"'(- 


ImW 

Re{c) 


(45) 


Where  the  unknown  constant  <f>  represents  the  phase  represented  by  the  propagation 
distance  between  the  antenna  and  object. 

Finally,  the  calibrated  scattering  matrix  of  an  unknown  object  (Pu)  can  be  found 
from  the  measured  scattering  matrix  (Nu)  by  the  expression: 


P»  =e-^'i7R-1NuT-1 


(46) 


3.4  Bistatic  via  Monostatic  Calibration 


AataWad  Uonoalntlc  Ppw*r  M  P^i  ^ 

*•»*•(»  BtafirflC  PRyi 

Figure  8:  BICOMS  Calibration  Geometry  (Alexander,  et  al,  1995) 

Used  by  the  Bistatic  Coherent  Measurement  System  (BICOMS)  at  Holloman  AFB, 

NM,  this  Type-3  calibration  technique  utilized  the  calibration  procedure  described  in 
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(Whitt,  et  al,  1991).  This  procedure  is  designed  around  the  unique  measurement 
conditions  of  this  facility.  Initially,  each  monostatic  leg  of  the  bistatic  path  is  calibrated 
separately,  as  illustrated  in  Figure  8.  The  normalized  full-polarimetric  distortion  matrices 
(PDM’s),  R  and  T,  are  solved  for  in  each  case,  exactly  as  outlined  in  Section  3.3, 
yielding  the  relationship  stated  in  equation  (28)  for  each  monostatic  leg.  This  process 
gives  four  matrices;  a  receiver  PDM  and  a  transmitter  PDM  for  each  of  the  two  antennas 
as  follows: 


Ri  = 


•vh 


V  Y 
.  hv  w. 


(47) 


R2  = 


*  vh 


Y  Y 

hv  w 


(48) 


where  rpq=r’pq/r’hh,  and  tpq^’pq/t’hh  (Tuley,  1995).  The  above  equations  need  only  solve 
for  the  PDMs  within  a  complex  constant,  the  exact  scattering  solutions  for  the  three 
monostatic  reference  objects  offers  no  useful  information,  since  the  amplitude  and  phase 
constant  for  the  bistatic  path  cannot  be  found  from  the  monostatic  path.  The  advantage  of 
this  technique  is  realized  when  objects  are  selected  such  that  the  structure  of  the 
scattering  matrix  is  well  known,  but  the  precise  solution  is  not.  For  instance,  the 
theoretical  scattering  matrix  that  would  be  input  to  the  procedure  for  a  sphere,  a  dihedral 
oriented  vertically,  and  a  trihedral  would  be 


'l 

0" 

'0 

0" 

'-1 

-f 

^  sphere 

_0 

1_ 

S  r= 

u  vertical, dihedral 

0 

1 

^  tilted, dihedral 

-1 

1 

(49) 


for  every  frequency. 
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At  this  point  the  necessary  PDMs  for  each  bistatic  path  are  known,  relative  to  the 
gain  of  the  VV  channel.  A  separate  calibration  compensation  step  is  performed  on  the 
bistatic  leg  to  solve  for  the  constant  amplitude  and  phase  factor.  An  object  is  measured 
which  has  a  known  W-polarized  return;  in  the  case  of  the  BICOMS  calibration  method, 
the  end-cap  of  a  large  squat  cylinder  is  used  to  determine  the  phase  factor,  and  the  same 
cylinder  is  measured  to  determine  the  amplitude  factor,  but  rotated  such  that  the  curved 
surface  is  oriented  forming  a  peak  specular  return  along  the  bistatic  path.  The  PO 
solution  for  a  circular  disk  is  used  as  a  theoretical  prediction  for  the  end-cap  cylinder. 
The  ratio  of  the  theoretical  and  measured  amplitudes  of  the  cylinder’s  curved  surface,  and 
the  difference  of  the  theoretical  and  measured  phases  for  the  end-cap  cylinder  provide  the 
necessary  constants  c  and  (f>  for  the  VV  channel, 


\M„ 


c  = 


\ cylinder 


v, cylinder 


</>  = 


ZM 


w,  end-cap 


ZS 


w, end -cap 


(50) 


and  as  such  there  exists  an  absolute  PDM  for  each  antenna  subsystem. 

This  technique,  while  providing  a  way  to  perform  the  general  Type-3  calibration 
technique  without  a  precise  theoretical  prediction  for  the  three  monostatic  reference 
objects,  has  the  disadvantage  that  it  requires  additional  objects  to  compensate  for  absolute 
amplitude  and  phase.  For  an  RCS  measurement  facility  such  as  BICOMS,  this  is  a 
significant  deficiency;  it  requires  more  time  in  mounting  the  additional  objects,  and 
system  drift  and  environmental  considerations  will  tend  to  degrade  the  performance  of 
the  calibration  over  that  time. 
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4  METHODOLOGY 


4.1  Purpose 

Since  the  purpose  of  this  work  is  primarily  to  assess  the  performance  of  several 
existing  and  proven  bistatic  calibration  techniques  the  techniques  are  evaluated  on  the 
basis  of  the  factors  listed  in  Section  1.  Again,  they  are: 

•  Sensitivity  of  calibration  to  reference  object  alignment  error 

•  Sensitivity  to  antenna  polarization  impurity 

•  Robustness  and  repeatability  of  the  technique 

•  Applicability  to  fixed-bistatic  or  true  bistatic  measurements 

•  Ease  and  efficiency  of  implementation 

•  Sensitivity  to  the  condition  of  the  reference  object  scattering  matrix 

These  are  the  major  concerns  for  the  accuracy  and  efficiency  of  a  given  calibration 
technique.  Of  course,  given  the  wide  variance  in  measurement  conditions  from  one 
measurement  range  to  another,  a  particular  calibration  technique  that  is  optimal  for  one 
range  may  be  unreasonable  for  another.  It  is  for  this  reason  that  calibration  techniques 
that  have  been  or  are  being  used  in  a  wide  variety  of  measurement  facilities  and  scenarios 
are  evaluated. 

In  addition  to  evaluating  these  techniques,  another  goal  of  this  research  is  to  provide 
documentation  of  the  specific  problems  and  benefits  of  using  a  particular  calibration 
object.  Many  concerns,  such  as: 

•  Ability  to  produce  precise  theoretical  predictions  for  the  object 

•  Alignment  sensitivity  of  the  object 

•  Object  signal-to-noise  ratio 
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•  Other  electromagnetic  scattering  characteristics  of  the  object 

are  all  important  information  in  selecting  an  object  that  will  yield  a  precise  calibration.  A 
reference  calibration  object,  for  instance,  nearly  always  requires  a  precise  calculation  of 
the  exact  RCS  amplitude  and  phase.  Reference  objects  that  are  more  complex  than  the 
canonical  set  of  spheres,  disks,  or  small  cylinders,  require  computationally  intense 
Moment  Method  simulations,  which  sometimes  are  impractical,  depending  on  the 
complexity  and  electrical  size  of  the  object. 

Objects  such  as  disks  and  long  cylinders,  though  often  easy  to  calculate  a 
prediction,  tend  to  be  very  sensitive  to  alignment  errors.  It  is  desired  that  reference 
objects  be  selected  such  that  they  are  rather  insensitive  to  alignment,  because  the  precise 
alignment  of  an  object  is  often  tedious  and  time-consuming,  and  also  tends  to  necessitate 
object  pedestals  that  are  large  and  interact  more  with  the  object. 

Objects  must  also  be  chosen  so  that  the  level  of  their  RCS  is  high  enough  to  be 
significantly  above  the  noise  and  clutter.  For  calibration  at  a  generic  bistatic  angle,  few 
objects  are  known  which  have  a  relatively  constant  RCS  for  all  bistatic  angles. 

Other  scattering  phenomenon,  such  as  creeping  waves,  can  also  adversely  affect  the 
suitability  of  a  reference  object.  If  an  object  is  chosen  with  a  large  creeping  wave 
component,  the  phase  of  the  RCS  will  not  be  a  suitable  reference  for  the  range  focus.  In 
addition,  any  surface  wave  on  the  object  could  be  retarded  when  it  comes  in  contact  with 
the  pedestal,  which  would  create  a  deviation  in  the  measurement  from  the  theoretical 
prediction. 
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4.2  Resources  Required 


In  order  for  the  factors  listed  in  Section  4.1  to  be  properly  evaluated,  the  work  was 
performed  in  a  facility  which: 


•  Has  the  capability  of  precisely  aligning  the  reference  obj  ects 

•  Has  very  good  antenna  polarization  purity  (at  least  30  dB) 

•  Measurements  in  the  facility  are  highly  repeatable,  necessitating  a  large  signal-to- 
noise  ratio 

•  Has  the  capability  to  perform  fixed-bistatic  or  true  bistatic  measurements 

•  Is  relatively  easy  to  mount  objects  and  configure  the  system 

•  Has  the  capability  to  perform  measurements  on  a  wide  variety  of  obj  ects 

The  European  Microwave  Signature  Laboratory  (EMSL)  was  chosen  because  it  fits 
most  of  these  criteria,  excepting  that  the  signal-to-noise  ratio  is  less  than  desirable  unless 
large  integration  is  used.  An  auxiliary  facility  is  also  available  for  the  duplication  of  the 
EMSL  measurements  at  the  bistatic  RCS  lab  of  the  Air  Force  Institute  of  Technology 
(AFIT)  at  Wright-Patterson  Air  Force  Base,  Ohio.  This  facility  does  not  have  the 
capability  to  collect  full-polarimetric  data,  nor  does  it  have  the  degree  of  accuracy  like  to 
the  EMSL.  These  factors,  among  others,  make  the  difference  between  the  facilities 
sufficient  enough  to  justify  its  use  in  extending  the  applicability  of  these  calibrations  to  a 
different  measurement  environment. 

As  well  as  an  accurate  and  versatile  measurement  facility,  it  is  also  necessary  to 
have  expertise  available  in  computational  scattering  and  RCS  measurements.  The 
theoretical  predictions  for  each  of  the  test  objects  are  a  major  portion  of  the  project. 
Expertise  in  the  Method  of  Moments,  as  well  as  Physical  Optics  and  the  Uniform 
Geometrical  Theory  of  Diffraction  was  necessary  in  order  to  obtain  all  the  theoretical 
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predictions  that  are  required.  The  proposed  object  test  matrix  is  varied  enough  such  that 
all  three  techniques  are  required. 

Experience  in  the  modeling  of  electromagnetic  scattering  as  well  as  RCS 
measurements  are  available  at  the  Technologies  for  Detection  and  Positioning  (TDP)  unit 
of  the  Space  Applications  Institute  (SAI)  at  JRC  as  well  as  AFIT. 

Not  least  of  all  the  requirements  for  this  project  was  time.  Collecting  the  raw  data 
for  the  calibration  test  objects  took  many  weeks  of  continuous  measurement.  The 
availability  of  range  time  in  the  EMSL  was  maximal. 

4.3  Approach 

4.3. 1  Assessment  of  Air  Force-wide  calibration  needs 

The  concern  of  the  Air  Force  in  the  field  of  bistatic  calibration  lies  naturally  in  its 
applicability  to  bistatic  RCS  measurement  facilities  currently  in  use  by  the  USAF. 
Facilities  currently  in  use  or  in  development  and  are  the  major  concern  of  this  work 
include  the  Bistatic  Coherent  Measurement  System  (BICOMS)  at  NRTF/RATSCAT  at 
Holloman  AFB,  NM,  the  dual  compact  range  under  construction  at  AFRL/SNHE  at 
Hanscom  AFB,  MA,  and  the  RCS  laboratory  at  the  Air  Force  Institute  of  Technology 
(AFIT)  at  Wright-Patterson  AFB,  OH.  None  of  these  facilities  are  suitable  for  a  thorough 
study  of  calibration  involving  the  specific  concerns  of  each  range,  so  another,  more 
versatile  facility  is  needed  to  complete  the  measurements,  namely  the  EMSL.  To  sum  up 
some  of  the  main  concerns  for  the  calibration  of  bistatic  radar  systems,  the  following  are 
of  considerable  importance: 
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•  Number  of  calibration  object  mounts  required 

•  Length  of  time  to  perform  calibration 

•  Ability  of  calibration  to  account  for  1st  and  2nd  order  cross-talk  effects  (full- 
polarimetric  calibration) 

•  Ease  and  practicality  of  producing  theoretical  predictions  for  reference  objects 

A  large  part  of  the  solution  to  these  problems  lies  in  the  careful  selection  of  the 
calibration  reference  objects. 

4.3.2  Calibration  object  selection 

There  is  a  range  of  objects  that  have  been  accepted  as  a  canonical  set  for 
calibration.  Monostatic  and  bistatic  calibration  objects  are  often  similar,  but  there  are 
objects  that  are  used  for  monostatic  calibration  that  are  not  suitable  in  a  bistatic  scenario. 
Some  popular  objects  for  monostatic  calibration  are  trihedrals,  dihedrals,  disks  and  plates, 
spheres,  and  cylinders.  All  objects  are  usually  metallic  with  a  brushed  or  polished  finish. 
For  bistatic  calibration,  trihedrals  are  not  suitable,  and  dihedrals,  though  in  some  cases 
practical,  are  not  as  useful  as  in  monostatic.  Disks,  plates,  and  cylinders  make  good 
calibration  objects  at  small  bistatic  angles,  with  varying  degrees  of  alignment  sensitivity 
among  the  three.  Due  to  their  complete  insensitivity  to  alignment,  and  the  available  exact 
analytic  solution  (the  Mie  Series),  the  sphere  is  a  universally  good  calibration  object, 
though  the  problem  of  low  RCS  at  any  given  bistatic  angle  is  an  ever-present  problem. 
Though  not  considered  a  canonical  object,  the  wire  mesh  has  also  proven  to  be  a  good 
object  for  polarimetric  calibration,  and  as  a  result  this  object  is  included  in  the  analysis. 

A  set  of  calibration  objects  was  chosen  based  on  their  use  in  bistatic  ranges  of  Air 
Force  interest,  and  additional  non-canonical  objects  were  chosen  for  comparison.  The 
bistatic  measurements  were  performed,  in  general,  at  a  relatively  small  bistatic  angle,  in 
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hopes  that  for  some  objects,  the  ease  and  efficiency  of  a  monostatic  calibration  can  be 
extended  to  a  bistatic  calibration  at  a  small  bistatic  angle.  Table  1  lists  the  objects  that 
were  tested. 


Table  1:  Test  Matrix 


Object 

Orientation 

Measurement  Type 

Metallic  sphere 

swept  bistatic 
monostatic 

long  cylinder 

d  =  50  mm,  h  =  200  mm 
axis  parallel  to  x-axis 
axis  45  deg  from  horizontal 
axis  45  deg  from  horizontal 
axis  67.5  deg  from  horizontal 
axis  67.5  deg  from  horizontal 

swept  bistatic 
bistatic:  4.6  deg 
monostatic 
bistatic:  4.6  deg 
monostatic 

Squat  cylinder 

d  =  215  mm,  h  =  108  mm 
axis  parallel  to  x-axis 
axis  45  deg  from  horizontal 
axis  45  deg  from  horizontal 
axis  22.5  deg  from  horizontal 
axis  22.5  deg  from  horizontal 

swept  bistatic 
bistatic:  4.6  deg 
monostatic 
bistatic:  4.6  deg 
monostatic 

Dihedral 

one  face:  1=300  mm,  w=1 50  mm 
Seam  oriented  vertical 

Seam  oriented  vertical 

Seam  45  deg  from  vertical 

Seam  45  deg  from  vertical 

Seam  22.5  deg  from  vertical 

Seam  22.5  deg  from  vertical 

bistatic:  4.6  deg 
monostatic 
bistatic:  4.6  deg 
monostatic 
bistatic:  4.6  deg 
monostatic 

Circular  disk 

20  cm  diameter 

Specular  orientation 

Specular  orientation 

bistatic:  4.6  deg 
monostatic 

wire  mesh 

19  cm  x  20  cm 

(19  20-cm  wires  spaced  1  cm  apart) 
Wires  oriented  vertically 

Tilted  45  deg  counter-clockwise 

Wires  oriented  vertically 

Tilted  45  deg  counter-clockwise 

bistatic:  4.6  deg 
monostatic 
bistatic:  4.6  deg 
monostatic 

Trihedral 

one  face:  a=102  mm,  c=145  mm 
Specular  orientation 

monostatic 
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For  each  object  and  orientation  listed  in  the  test  matrix,  measurements  were 
performed  for  each  of  the  two  Tx/Rx  combinations,  and  all  were  fully  polarimetric;  i.e. 
HH,  VH,  HV,  and  VV  pol  were  all  collected. 

4.3.3  Experimental  validation  of  measurement  facility 

In  order  for  a  valid  study  of  calibration  to  be  done,  the  idiosyncrasies  of  the 
measurement  facilities  must  be  first  understood.  The  design  of  the  chamber  is  the  most 
important  factor  to  consider.  Many  experimental  facilities  are  designed  with  a  specific 
puipose  in  mind.  In  most  cases,  an  RCS  measurement  facility  is  designed  around  a 
specific  object  or  set  of  objects  to  be  measured.  The  BICOMS  system  at  Holloman  AFB, 
NM,  for  instance,  was  designed  for  full-scale  polarimetric  measurements  of  aircraft 
(Alexander,  et  al.).  The  calibration  considerations,  namely  the  time  necessary  to  perform 
the  calibration,  for  a  range  such  as  this,  are  sometimes  much  different  than  those  of  an 
indoor  chamber.  Other  indoor  chambers  have  very  different  design  specifications,  as 
well.  Section  4.4  contains  a  discussion  of  the  design  of  the  EMSL  and  its  capabilities  for 
this  study. 

4.3.4  Validation  of  measurements  via  computational  methods 

To  ensure  that  the  data  taken  from  the  EMSL  are  reliable,  a  simple  validation  was 
performed  on  the  measurement  of  the  circular  disk.  The  PO  solution  for  the 
backscattering  from  a  metallic  disk  is  given  above  in  Equation  (13).  The  quasi¬ 
monostatic  measurement  of  the  disk  was  calibrated  using  the  EMSL  simple  polarimetric 
technique,  and  compared  against  the  PO  prediction.  As  seen  in  Figure  9,  the  RCS  error  is 
well  within  1  dB,  but  the  calibrated  phase  shows  a  marked  difference  in  slope  and  initial 
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phase.  This  is  due  to  the  fact  that  the  disk  was  not  located  at  exactly  the  phase  center 
defined  by  the  calibration.  The  deviation  here  shows  a  difference  in  position  of  about  Vi 
cm  (for  an  approximate  phase  difference  of  50°  at  10  GHz). 


RCS  Comparison  of  Disk  Measurement  to  PO  Prediction 


Phase  Comparison  of  Disk  Measurement  to  PO  Prediction 


Figure  9:  Validation  of  W-pol  Measurement  in  EMSL. 

4.3.5  Implementation  of  calibration  techniques 

All  of  the  calibration  techniques  listed  in  Section  4  were  implemented  as  MATLAB 
functions,  which  are  listed  in  Appendix  B. 

4,4  Characterization  of  the  EMSL 

Before  performing  the  measurements  on  the  calibration  test  objects,  the  EMSL  was 
thoroughly  evaluated  for  its  performance  of  the  criteria  listed  in  Section  4.2. 
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4. 4. 1  Design  of  the  EMSL 


The  EMSL  was  designed  as  a  versatile  radar  signature  measurement  facility,  with 
the  capability  to  stable  environmental  conditions,  as  well  as  various  operating  modes 
including  one-dimensional  scatterometer,  two-dimensional  SAR,  and  three-dimensional 
tomogram  measurements  (Appendix  C).  Figure  10  is  an  exploded  view  of  the  layout  of 
the  EMSL. 


The  European  Microwave  Signature  Laboratory  (EMSL) 


1.  Main  Door 

2.  Service  door 

3  Target 

4 .  Target  Positioner 

5.  Target  Positioner 
rail 

6.  Radar  Sleds 

7  Sled  rail 

S .  F ere  d  Antennae 

9.  Free  standing 
structure 
10.  Microwave 
absorber  cover 


Figure  10:  Physical  Layout  of  the  E1MSL 

The  distance  from  the  antennae  to  the  range  focus  is  10  m.  One  of  several  object 

positioners  of  varying  heights  can  be  selected  based  on  the  size  of  the  object  to  be 
measured.  To  increase  the  ease  of  mounting  heavy  objects,  the  object  positioner  can  be 
moved  outside  the  range  through  the  main  door,  and  a  crane  is  used  to  place  the  object 
upon  the  positioner.  For  small  objects  such  as  the  canonical  objects  used  in  this  work  for 


far-field  measurements,  a  mechanical  scissor-lift  can  be  brought  inside  the  range  so  that 
the  object  can  be  manually  placed  and  aligned  with  relative  ease.  The  interior  of  the 
chamber  is  covered  with  pyramidal  microwave  absorber  to  attenuate  the  reflections  from 
radiation  incident  upon  the  chamber  walls. 

4.4.2  Alignment  of  reference  objects 

The  EMSL  uses  lasers  mounted  on  the  antenna  sleds  to  align  the  objects.  The  laser 
is  located  at  the  center  of  the  sled,  and  points  to  the  center  of  the  range  focus  to  within 
about  0.1  degrees.  The  transmit  and  receive  antennae,  however,  are  located  at  a  position 
offset  from  the  laser,  such  that  the  angular  separation  between  the  lasers  of  both  sleds  is 
the  same  as  the  angular  separation  between  the  transmit  and  receive  antennae  on  their 
corresponding  sleds  (within  the  plane  of  sled  movement).  When  performing  a  laser 
alignment  for  an  object  such  as  the  calibration  disk,  the  laser  of  sled  B  is  activated  and 
the  disk  is  oriented  such  that  the  laser  reflects  off  the  disk  and  onto  the  center  of  sled  A. 
The  position  of  the  receive  sled  (usually  sled  A)  must  then  be  corrected  to  account  for  the 
offset  of  the  antennae  with  respect  to  the  laser  (Figure  11).  In  this  figure  the  dotted  red 
line  represents  the  laser  (transmitted  from  sled  B  to  sled  A),  and  the  dark  blue  line  is  the 
path  of  the  RF  energy  as  it  reflects  off  the  disk.  It  is  necessary  to  close  the  angular 
separation  between  the  sleds  by  1 .4  degrees  in  order  to  align  the  main  lobe  of  the  specular 
reflection  of  the  disk  onto  the  receive  antenna  of  sled  A. 
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SLED  B 


TARGET  (disk) 


Figure  11:  Alignment  of  Calibration  Disk  in  the  EMSL 


With  this  laser  alignment  procedure,  the  alignment  of  a  circular  metallic  disk  of  20 
cm  diameter  can  be  repeated  to  with  0.1  dB.  Figure  12  shows  two  independent 
alignments  of  this  disk,  separated  by  over  a  week’s  time,  and  with  completely 
independent  pedestal  placements.  The  measurements  are  both  gated  and  subtracted. 
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Independent  Measurements  of  Small  Calibration  Disk:  W  pol 


Figure  12:  Repeatability  Comparison  of  Calibration  Disk  Measurements 

The  0.2  dB  difference  between  these  measurements  is  likely  not  due  to  alignment 

error,  but  rather  system  drift.  If  alignment  error  was  the  source,  then  the  disparity 
between  the  measurements  would  increase  with  frequency,  as  the  electrical  size,  and 
therefore  the  alignment  sensitivity  of  the  disk,  increased.  From  the  results  of  this 
analysis,  it  can  be  safely  assumed  that  alignment  repeatability  will  yield  results  well 
within  the  0.5  dB  measurement  error  limit  that  is  advertised  in  the  EMSL  (Nesti,  et  al, 
1994). 

4.4.3  Antenna  polarization  purity 

Good  antenna  polarization  purity  is  an  important  characteristic  to  the  range  in 
which  the  test  measurements  are  being  taken.  Since  the  goal  is  to  produce  results  with 
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wide  applicability,  to  include  measurement  systems  with  poor  cross-polarization  purity  as 
well  as  those  with  good  cross-polarization  purity,  this  range  must  be  able  to  simulate  both 
of  these.  For  comparison,  one  can  simulate  poor  cross-polarization  purity  by  a  matrix 
rotation  of  the  measurement.  The  purity  of  the  EMSL  antennae  is  between  30  and  40  dB 
(Nesti,  et  al,  1994).  To  verify  this,  a  measurement  of  a  metallic  sphere  was  taken.  The 
sphere  has  no  theoretical  cross-polarized  RCS;  therefore  the  level  of  cross-polarized  radar 
return  relative  to  the  RCS  of  the  co-polarized  channels  represents  the  antenna  cross¬ 
polarization  impurity.  Figure  13  shows  a  measurement  of  such  a  sphere  from  5-15 
GHz,  which  is  the  frequency  range  of  interest  for  this  work. 


Cross  Polarization  Purity  of  EMSL  for  Monostatic  Sphere  Measurement 


Figure  13:  Antenna  Cross-Polarization  Purity 
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The  separation  between  the  co-polar  and  cross-polar  components  of  the 
measurements  is  about  30  dBsm.  As  mentioned  above,  the  cross-polar  RCS  that  is 
measured  is  representative  of  the  antenna  cross-polarization  impurity.  Therefore  it  can  be 
said  that  the  cross-polarization  purity  of  the  EMSL  antennae  is  about  30  dBsm. 

4. 4. 4  Signal-to-noise  ratio  and  repeatability 

The  EMSL  has  no  hardware  gating  capability.  As  a  result,  ambiguous  clutter 
(clutter  that  occurs  outside  the  range  cell)  in  the  measurement,  which  arise  as  a  result  of 
multiple  reflections  from  the  chamber  walls  and  mounting  apparatus,  as  well  as  antenna 
coupling  decrease  the  signal-to-noise  ratio  (SNR).  Increasing  the  integration,  or 
averaging  of  each  individual  measurement  can  increase  the  SNR.  Averaging  over  2N 
points  increases  the  SNR  by  a  factor  of  N  dBsm  relative  to  no  averaging.  However, 
averaging  beyond  1 024  points  becomes  time-prohibitive  and  decreases  the  quality  of  the 
background  subtraction  due  to  the  time  separation  between  the  object  measurement  and 
empty-room  measurement.  To  evaluate  the  repeatability  of  a  given  microwave 
measurement  and  excluding  the  factor  of  alignment,  independent  measurements  (i.e., 
separate  mountings,  and  separated  by  more  than  a  day’s  time)  are  taken  for  an  object  that 
is  completely  insensitive  to  alignment.  The  metallic  sphere  is  such  an  object.  Figure  14 
shows  two  independent  measurements  and  the  coherent  difference  of  the  two.  For  these 
measurements,  2048  averaging  points  were  used.  The  data  were  calibrated  with  identical 
distortion  matrices.  This  amount  of  averaging  is  excessive  for  most  bistatic 
measurements  in  the  EMSL,  as  it  degrades  the  quality  of  the  background  subtraction. 
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Sphere  Measurements  and  Their  Difference:  W  pol,  2048-averaging 


Figure  14:  Repeatability  of  a  Measurement  of  the  30.5  cm  Metallic  Sphere,  a=120° 

The  difference  between  the  measurements  of  Figure  14  always  increases  with 

frequency.  This  is  due  to  the  increased  SNR  at  high  frequencies,  as  a  result  of  the  loss  in 
the  coaxial  cables.  At  the  highest  frequency  of  interest  (15  GHz),  the  measurement  can 
be  repeated  to  within  a  difference  of -12  dB. 

One  can  get  a  better  idea  of  the  actual  SNR  of  the  range  by  eliminating  other  factors 
such  as  mounting  and  time  separation.  Two  dependent  measurements  (i.e.  a  single 
mount,  very  little  time  separation  between  the  two)  are  taken  and  the  results  of  the  two 
are  subtracted.  Figure  15  shows  these  measurements.  The  difference  between  the  two 


50 


represents  the  noise  floor  for  512-point  averaging. 

Noise  Floor  Derived  from  Consecutive  Sphere  Measurements 


Figure  15:  Noise  Floor  Derived  from  Sphere  Measurement  (512-averaging) 

From  the  figure,  it  can  be  seen  that  the  noise  floor  does  not  exceed  -85  dBsm.  The 

noise  floor  of  the  system  (that  is,  without  averaging)  can  be  calculated  by  subtracting  the 
effect  of  averaging.  Given  that  this  measurement  uses  (2  9  =)  512-point  averaging,  9  dB 
should  be  added  making  the  system  noise  floor  about  -76  dB. 

4.4.5  Ease,  efficiency,  and  versatility 

Pedestal  design  is  the  most  important  factor  in  the  ease  and  efficiency  of  mounting 
an  object  for  measurement.  If  pedestals  can  be  designed  such  that  objects  can  be 
switched  on  the  object  column  without  the  need  of  any  other  modifications  to  the  setup, 
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the  time  it  takes  to  move  the  objects  is  minimal.  Because  of  the  EMSL’s  large  size, 
mechanical  lifts  are  required  to  facilitate  the  movement  of  objects,  and  as  a  result  safety 
becomes  a  major  concern.  Measurement  setup  is  rather  efficient  with  two  or  more  people 
operating  the  range.  After  the  measurement  is  set  up,  the  EMSL  uses  an  automated  data 
acquisition  and  control  system  to  allow  it  to  perform  a  complete  measurement  without 
human  intervention. 

The  EMSL  was  designed  to  measure  a  large  variety  of  objects.  The  object  column 
can  support  objects  of  up  to  six  tons  in  weight.  For  far-field  measurements,  the 
maximum  linear  dimension  of  the  object  must  not  be  more  than  about  27  cm.  For  a 
direct-illumination  range  (one  without  reflectors)  this  is  a  very  sizeable  quiet  zone,  and 
can  easily  support  a  wide  variety  of  canonical  objects.  ROHACELL  51  foam  is  used  as 
mounting  material.  The  foam  is  much  more  rigid  than  commercial  Styrofoam,  and  has  a 
relative  dielectric  constant  of  sr  =  1.071  (Weinberger  and  Nesti,  1994).  The  rigidity  of 
this  foam  permits  the  use  of  heavy  objects,  and  the  foam  can  also  be  machined  into 
mechanical  parts  such  as  pivots  and  screw  threads  to  make  object  alignment  more 
efficient. 

4.5  Measurement 

For  each  measurement  of  Table  1 ,  a  systematic  alignment  error  was  introduced 
(with  the  exception  of  the  sphere  measurement).  As  shown  in  Figure  16  for  the  small 
disk,  alignment  error  was  introduced  from  0°  to  2°,  in  increments  of  0.5°. 


52 


Circular  Disk  (VV-pol):  Measurement  with  systematic  misalignment 


frequency  (GHz) 

Figure  16:  Measurement  of  circular  disk  with  systematic  misalignment 

Depending  on  the  object,  this  alignment  error  was  added  in  azimuth  or  elevation. 

The  choice  of  which  axis  to  rotate  the  object  was  made  in  accordance  to  which  would 
produce  the  most  destructive  error  in  the  measurement.  In  the  case  of  the  horizontally 
cylinders,  for  instance,  an  alignment  error  in  the  elevation  plane  would  not  produce  an 
appreciable  error,  given  that  the  measurement  is  in  the  quasi-monostatic  region.  An 
azimuthal  alignment  error,  therefore,  was  added.  For  the  dihedral  measurements,  an 
elevation  error  was  used  since  the  dihedral  in  the  vertical,  or  near-vertical  orientations  has 
a  relatively  flat  RCS  vs.  bistatic  angle  when  the  receiver  is  scanned  in  the  azimuth  plane. 
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5  RESULTS 


5.1  Dihedral  (Monostatic) 

In  a  monostatic  configuration,  a  dihedral  provides  an  excellent  reference  object  for 
the  purpose  of  cross-polarization  calibration.  The  theoretical  prediction  is  not  easy  to 
realize;  it  is  necessary  to  generate  a  Moment  Method  simulation.  However,  cross¬ 
polarization  calibration  can  be  accomplished  with  two  measurements  of  the  same 
dihedral  in  different  orientations,  without  knowledge  of  the  exact  theoretical  RCS 
prediction,  as  done  in  the  EMSL  simple  polarimetric  calibration  technique  (3.2). 

The  scattering  matrix  of  the  dihedral  oriented  with  the  seam  parallel  to  the  vertical 
axis  (the  direction  of  vertical  polarization)  looks  like: 


S  =  S 


vert 


0 
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(51) 


For  a  dihedral  composed  of  two  identical  faces  at  right  angles  to  each  other,  with  face 
dimensions  a  and  b,  the  RCS  (Sverl)  is: 
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The  scattering  matrix  of  the  same  dihedral  rotated  in  the  plane  defined  by  the  direction  of 
the  wave  propagation  can  be  found  by  a  rotation  of  the  scattering  matrix  of  Eq.  (51), 
explained  in  (Mott,  1992).  At  a  rotation  through  angle  0,  the  scattering  matrix  becomes: 


and  if  0  =  22.5°,  then, 
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The  advantage  of  using  the  dihedral  in  an  orientation  such  that  the  scattering  matrix  is  of 
the  form  of  Eq.  (54)  is  twofold:  first,  the  absolute  amplitude  need  not  be  known  if  the 
amplitude  in  the  0  =  0°  orientation  is  known  (Svert  is  known);  second,  the  elements  of  the 
scattering  matrix  each  have  equal  amplitude,  yielding  strong  cross-polarized  response,  as 
well  as  making  it  linearly  independent  from  all  co-polarizing  scatterers. 


The  disadvantage  of  using  the  dihedral  as  a  calibrator  is  its  alignment  sensitivity. 
For  measurement  systems  which  use  two  closely-separated  antennas  for  monostatic 
measurements  (such  as  the  EMSL),  this  antenna  separation  can  be  significant  enough  so 
that  the  RCS  of  each  polarization  channel  is  not  balanced  as  Eq.  (51)  and  (54)  would 
indicate.  In  the  case  of  the  EMSL,  the  antennae  on  each  sled  are  separated  by  about  1 .7°. 
Figure  17  is  the  MoM  prediction  of  a  measurement  in  this  configuration.  It  can  be  seen 
that  the  difference  in  the  co-polarized  channels  reflects  a  difference  of  up  to  1  dB  at  lower 
frequencies.  The  cross-polarized  channels  are  non-zero,  but  negligible. 


55 


MoM  Simulation  of  Vertical  Dihedral  (Monostatic  Sled  A) 


Figure  17:  MoM  simulation  of  an  EMSL  monostatic  measurement  of  the  vertical 

dihedral 

These  results  indicate  that  the  dihedral  is  best  suited  for  monostatic  cross¬ 
polarization  calibration  of  a  radar  system  that  has  a  single  Tx/Rx  antenna,  rather  than  two 
closely  separated  antennae. 

For  a  Type-2  calibration,  the  ratio  of  the  RCS  of  the  vertical  dihedral  to  the  tilted 
dihedral  is  important.  Using  Eq.  (51)  and  (54),  one  can  predict  the  trend  vs.  frequency  of 
this  ratio. 


c 

°  vertical 

C 

°  tilted 


1  0 
0  1 


(55) 


If  a  dihedral  measurement  such  as  this  were  to  be  used  in  a  Type-2  calibration  (one 
for  which  the  exact  scattering  matrices  of  the  cross-polarizing  reference  objects  is  not 
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known),  the  ratio  of  the  RCS  the  co-polarizing  orientation  to  the  RCS  in  the  cross- 
polarizing  orientation  (in  the  case  of  the  dihedral,  0°  and  22.5°  from  vertical, 
respectively),  the  trend  of  RCS  vs.  frequency  shown  in  Figure  18  results: 


RCS  ratios  of  Vertical  Dihedral  to  Dihedral  Tilted  22.5  deg  cw 


Figure  18:  Calibrated  RCS  ratios  of  measurement  of  dihedral  in  two  orientations — 

performed  in  the  EMSL 

The  calibrated  ratios  of  this  figure  would  be  relatively  flat  in  a  true  monostatic 
scenario.  The  rippling  due  to  the  small  separation  between  the  Tx  and  Rx  antennas  in  the 
EMSL  makes  these  measurements  inappropriate  for  cross-polarization  calibration  in  a 
Type-2  technique,  since  the  RCS  of  the  tilted  dihedral  relative  to  the  RCS  of  the  vertical 
dihedral  for  each  polarization  channel  is  a  smooth  trend  vs.  frequency,  as  indicated  by 
Eq.  (54). 
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For  a  Type-3  calibration,  a  dihedral  measured  with  a  similar  angular  separation 
between  the  antennas  would  be  perfectly  suitable.  It  is  only  necessary  for  this  type  of 
calibration  that  the  depolarizing  object  has  a  high  cross-polarized  RCS — a  requirement 
that  is  met  in  the  tilted  dihedral,  even  at  a  small  bistatic  angle. 

5.2  Dihedral  (Bistatic) 

As  one  might  predict,  the  dihedral  measurement  at  a  larger  bistatic  angle  (this  time, 
using  Tx  and  Rx  antennae  on  separate  sleds,  creating  a  bistatic  angle  of  about  6°),  creates 
an  even  larger  deviation  from  the  ideal  monostatic  response.  The  tilted  dihedral  still 
yields  a  large  cross-polarized  RCS,  which  makes  it  suitable  for  a  Type-3  calibration,  but 
the  vertical  dihedral  begins  to  return  higher  RCS  in  the  cross-polarized  channels.  This 
phenomenon  creates  problems  for  the  calibration,  since  it  is  no  longer  a  good 
approximation  to  set  the  cross-polar  RCS  to  zero  in  the  theoretical  prediction;  it  instead 
must  be  solved  for.  But  the  low  levels  of  the  cross-polarized  RCS  make  it  difficult  for  an 
iterative  MoM  solver  code  to  calculate  with  accuracy  (see  Figure  17,  specifically  the  HV- 
pol  channel).  This  effect  will  not  just  degrade  the  cross-polarization  calibration 
capability  of  a  Type-3,  but  will  degrade  the  overall  performance  of  the  calibration. 

5.3  Long  Cylinder 

The  long  cylinder  (that  is,  a  cylinder  which  has  a  length/diameter  ratio  of  greater 
than  two)  offers  a  higher  specular  RCS  than  a  sphere  of  the  same  diameter,  but  suffers  the 
problem  of  alignment  sensitivity — more  so  than  is  true  of  the  squat  cylinder.  Figure  19 
shows  the  measurement  of  the  long  cylinder  with  the  axis  of  rotation  oriented 
horizontally,  or  parallel  to  the  E-field  of  the  H-polarization  channel.  The  results  of  the 


bistatic  measurement  (using  both  antenna  sleds)  are  similar  to  the  monostatic 
measurement  (using  a  single  sled),  with  the  exception  that  the  SNR  of  the  bistatic 
measurement  is  higher,  which  is  a  result  of  poor  amplification  characteristics  of  the 
antennas  on  sled  A.  The  levels  of  the  cross-polarized  channels  (theoretically  zero)  give 
some  indication  of  the  cross-polarization  impurity  of  the  measurement. 

Horizontal  long  cylinder  bistatic  (5  deg  separation)  Horizontal  long  cylinder:  monostatic  (sled  B) 


frequency  (GHz)  frequency  (GHz) 


Figure  19:  EMSL  measurements  (Type-2  calibrated)  of  the  long  cylinder 

Figure  20  compares  the  co-polarized  channels  of  this  measurement  with  the 

predicted  MoM  solution.  The  agreement  of  the  co-polarized  channels  is  within  1  dB  over 
5-15  GHz. 
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Horiz.  long  cylinder:  bistatic 


Horiz.  long  cylinder:  monostatic 
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Figure  20:  EMSL  measurements  of  Figure  19  vs.  MoM  simulations 

A  possible  reason  for  using  a  long  cylinder  as  a  calibration  reference  object  would 

be  the  ability  to  produce  a  large  cross-polarized  RCS  component  when  the  object  is 
rotated.  A  thin  cylinder,  when  tilted  within  the  plane  defined  by  the  incident  plane  wave, 
when  the  axis  of  rotation  of  the  cylinder  is  parallel  to  that  plane,  will  produce  cross- 
polarized  RCS  return,  much  like  the  dihedral,  but  has  little  specular  return.  By  using  a 
fatter  cylinder  (this  long  cylinder  has  a  diameter  of  50  mm  and  a  length  of  200  mm),  the 
specular  RCS  is  increased,  at  the  cost  of  the  amplitude  of  the  cross-polarized  RCS  that 
can  be  returned.  Figure  21  shows  that  the  cross-polarized  RCS  of  such  a  cylinder  rotated 
45°  in  the  plane  of  the  incident  plane  wave  produces  a  cross-polar  RCS  below  about  -30 
dB,  which,  for  the  sensitivity  of  the  EMSL,  is  not  sufficient  for  cross-polarization 
calibration. 
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Long  Cylinder  tilted  45  deg:  monostatic 


frequency  (GHz) 

Figure  21 :  EMSL  measurement  of  long  cylinder,  tilted  45° 


5.4  Squat  Cylinder 

The  squat  cylinder  (a  cylinder  with  a  diameter/length  ratio  of  greater  than  two)  has 
a  higher  monostatic  RCS  than  the  long  cylinder  with  the  same  maximum  linear 
dimension.  This  is  important  given  that  the  maximum  linear  dimension  of  the  object  is 
the  limiting  factor  for  the  maximum  size  of  the  object  required  for  far-field 
measurements.  This  higher  RCS  in  monostatic  and  quasi-monostatic  measurements  is 
combined  with  the  fact  that  the  squat  cylinder  provides  much  less  cross-polarized  RCS 
when  tilted  off  the  axis  of  rotation,  as  seen  in  Figure  22,  in  which  case  the  RCS  of  the 
squat  cylinder  is  approximately  2  dB  higher  than  in  Figure  21,  and  the  cross-polarized 
RCS  component  is  low  enough  to  be  indistinguishable  from  the  ambient  noise. 
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Short  Cylinder  tilted  45  deg:  monostatic 


Figure  22:  EMSL  measurement  of  short  cylinder,  tilted  45° 

The  squat  cylinder  is  used  successfully  in  the  calibration  of  the  BICOMS  system 

(Tuley,  1995),  and  takes  advantage  of  the  flat  end-cap  of  the  cylinder  for  bistatic  phase 
calibration.  This  particular  usage  illustrates  each  of  the  strengths  of  this  cylinder  as  a 
bistatic  reference  object.  The  RCS  is  high  (higher  than  a  sphere  of  the  same  diameter)  in 
the  co-polarized  channels,  is  largely  insensitive  to  alignment  error  (illustrated  in  Figure 
23  through  the  comparison  of  the  bistatic  and  monostatic  measurements),  and  provides  a 
flat,  circular  surface  (the  end-cap)  to  provide  an  RCS  return  with  a  distinct  phase  peak. 
The  MoM  prediction  of  the  squat  cylinder  in  both  monostatic  and  bistatic  configurations 
(shown  in  Figure  24)  also  demonstrate  that  the  squat  cylinder  is  predicted  more  easily 
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than  for  the  long  cylinder  (using  a  Combined-Field  Integral  Equation  (CFIE)  and  an 
iterative  matrix  solver). 


Horizontal  short  cylinder:  bistatic  (5  deg  sep.)  Horizontal  short  cylinder:  monostatic  (sled  B) 


Figure  23:  EMSL  measurements  (Type-2  calibrated)  of  the  short  cylinder 
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Figure  24:  EMSL  measurements  of  Figure  23  vs.  MoM  simulation 
5.5  Trihedral 

The  high  RCS  return  of  the  dihedral  is  the  primary  reason  it  is  a  popular  monostatic 
calibration  reference  object.  The  cross-polarized  return  is  non-zero  in  this  case  (Figure 
25)  because  the  edges  of  the  trihedral  are  manufactured  such  that  two  of  the  three  are 
aligned  diagonally  with  respect  to  the  horizontal  and  vertical  directions.  The  diffraction 
from  these  edges  reflects  a  cross-polarized  RCS. 

The  experiment  published  in  (Peters,  1962)  provides  an  analysis  of  this  type  of 
trihedral  (called  a  triangular  dihedral)  as  a  bistatic  calibration  device.  The  RCS  of  the 
trihedral  becomes  highly  rippled  and  dramatically  lower  as  the  bistatic  angle  is  increased 
past  about  6°.  The  high  degree  of  rippling  makes  the  trihedral  extremely  difficult  to  align 
such  that  the  RCS  is  predictable.  From  this  it  can  be  stated  that  the  triangular  trihedral  is 
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not  suitable  at  all  for  bistatic  calibration— only  for  monostatic  co-polar  amplitude 
calibration. 


Trihedral:  monostatic  (sled  B) 


Figure  25:  EMSL  measurement  of  the  trihedral 

Figure  25  displays  a  monostatic  measurement  of  the  trihedral.  The  cross-polarized 

RCS  is  non-zero,  but  low  enough  to  be  indiscernible  from  the  clutter  and  polarization 
channel  cross  talk.  The  cross-polarized  RCS  levels  are  the  result  of  edge  diffraction  from 
the  forward  edges  that  are  oriented  diagonally  with  respect  to  the  antenna.  The  low  levels 
of  the  cross-polar  RCS  make  it  hard  to  predict,  and  thus  unsuitable  for  any  Type-3 
calibration,  because  an  accurate  theoretical  prediction  cannot  be  produced.  The  use  of  a 
square  dihedral,  where  the  edges  are  all  parallel  or  perpendicular  with  respect  to  the 
antenna,  would  alleviate  this  problem. 
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Figure  26:  EMSL  measurement  of  Figure  25  vs.  MoM  simulation 

The  MoM  simulation  on  the  trihedral  produces  the  desired  accuracy  (within  1  dB), 

as  seen  in  Figure  26.  The  lobing  in  the  frequency  domain  makes  the  trihedral  sensitive  to 
misalignment,  however  the  lobing  is  shallow  enough  so  that  the  co-polar  RCS  remains 
highly  predictable. 


5.6  Circular  Disk 

The  combination  of  high  co-polar  RCS  and  a  readily  available  theoretical  prediction 
(via  Physical  Optics)  make  the  circular  disk  an  obvious  choice  for  monostatic  or  quasi¬ 
monostatic  calibration.  The  disk  is  extremely  alignment-sensitive,  however,  and 
alignment  verification  can  add  precious  time  to  a  calibration  procedure.  In  the  EMSL, 
which  uses  a  precise  laser  alignment  system,  the  circular  disk  has  shown  to  be  a  very 
valuable  object.  The  cross-polar  RCS  return  due  to  diffraction  from  the  disk’s  edge  is 
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very  low  relative  to  the  co-polar  RCS,  giving  the  scattering  matrix  a  nearly  pure  co-polar 
response.  This  response  is  displayed  in  Figure  27. 


Small  disk:  bistatic  (5  deg  separation) 


frequency  (GHz) 

Figure  27:  EMSL  measurement  of  the  small  disk  vs.  PO  prediction 
5.7  Sphere 

The  rotational  symmetry  of  the  sphere  about  any  axis  alleviates  almost  all 
alignment  problems  normally  associated  with  a  calibration  measurement.  Since  the 
sphere  needs  to  be  fairly  large  to  produce  an  RCS  high  enough  to  calibrate  with,  object 
positioning  can  become  a  problem.  A  large  sphere  measured  at  a  small  bistatic  angle  is 
also  necessary  to  mitigate  the  creeping  wave  contribution.  The  creeping  wave  for  a 
sphere  that  occupies  all  or  exceeds  the  dimensions  of  the  quiet  zone  will  not  interact  with 
the  specular  return  in  a  manner  characteristic  of  a  far-field  measurement.  The  magnitude 


67 


of  this  creeping  wave  decreases  in  proportion  to  the  size  of  the  sphere,  and  is  further 
reduced  when  measured  at  small  bistatic  angles,  leaving  only  the  specular  RCS  return, 
which  is  contributed  primarily  from  only  the  front  of  the  sphere,  allowing  a  far-field  near¬ 
monostatic  measurement  to  be  performed  on  a  sphere  for  which  the  dimensions  exceed 
the  dimensions  of  the  chamber’s  quiet  zone. 


Sphere:  bistatic  (5  deg  separation) 
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Figure  28:  EMSL  measurement  of  the  30.5  cm  diameter  sphere 
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Sphere:  bistatic  (5  deg  separation) 


Figure  29:  EMSL  measurement  of  Figure  28  vs.  Mie  series  prediction 

As  seen  from  Figure  28  and  Figure  29,  the  RCS  of  the  sphere  is  relatively  constant 

versus  frequency  in  a  near-monostatic  configuration.  This  fact  brings  to  light  another 
possible  advantage  of  using  the  sphere  as  a  calibrator.  As  seen  in  the  figures,  the  signal- 
to-noise  ratio  of  the  measurement  is  not  large  enough  to  be  able  to  discern  the  lobing 
pattern  of  the  RCS  in  the  frequency  domain  (for  both  polarizations).  Essentially,  for  this 
bandwidth,  the  lobing  pattern  is  inconsequential.  In  creating  a  theoretical  prediction, 
then,  it  would  be  just  as  useful  to  use  a  single  complex  constant  for  this  entire  bandwidth, 
rather  than  the  vector  that  represents  the  true  RCS.  This  would  save  memory  and 
computational  time  in  calibrating  with  this  sphere  at  this  bandwidth. 
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5.8  Wire  Mesh 


The  wire  mesh  is  a  simple  arrangement  of  parallel  wires,  in  this  case  separated  from 
each  other  by  a  distance  of  1  cm.  The  mesh  consists  of  19  wires,  each  of  length  20  cm.  It 
acts  as  a  sort  of  polarization  filter,  in  that  only  the  polarization  component  of  the  incident 
wave  that  is  parallel  to  the  length  of  the  wires  will  be  reflected — all  other  components  are 
transmitted  through  the  mesh.  This  assumes  that  the  wires  are  thin  enough  such  that  the 
specular  return  from  the  other  polarization(s)  is  negligible. 


Vertical  wire  mesh:  bistatic  (5  deg  separation) 


Figure  30:  EMSL  measurement  of  wire  mesh  (vertical  orientation) 

Under  the  assumption  that  the  wires  are  both  long  and  thin  with  respect  to  the 

incident  wavelength,  and  further  that  the  main  lobe  of  the  RCS  is  wide  relative  to  the 
bistatic  angle  being  considered,  the  scattering  characteristics  as  a  function  of  the 
orientation  of  the  wires  is  approximated  by  the  relationship  of  Eqn’s  (19)  and  (20). 
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Tilted  wire  mesh:  bistatic  (5  deg  separation) 


frequency  (GHz) 

Figure  31 :  EMSL  measurement  of  wire  mesh  (oriented  45°  counter-clockwise) 

In  the  measurement  of  the  mesh  for  which  the  wires  are  aligned  vertically  (shown  in 

Figure  30),  the  high,  flat  response  of  the  VV  channel  is  immediately  apparent.  Like  the 
sphere,  this  is  also  a  good  characteristic  for  a  calibration  object.  The  large  HH-pol 
component  in  this  measurement  is  unexplained,  although  it  could  be  a  specular  response 
or  a  traveling  wave  effect.  The  effect  of  the  Styrofoam  substrate  on  the  measurement  of 
the  wires  has  not  been  considered,  and  may  be  necessary  for  further  work  on  this 
calibration  object. 

Figure  31  shows  the  same  wire  mesh  measured  at  an  orientation  45°  from  vertical. 
As  predicted,  the  RCS  of  each  polarization  channel  is  nearly  identical,  to  at  least  within  1 
dB.  The  ratio  of  the  RCS  versus  frequency  for  the  vertical  and  tilted  orientations. 
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however,  is  not  constant,  which  becomes  apparent  by  a  visual  comparison  of  each  plot. 
Since  this  object  is  used  only  as  a  cross-polarization  compensation  object  for  a  Type-2 
calibration,  however,  this  ratio  need  not  be  as  precise  as  the  co-polarized  amplitude 
calibrator,  since  the  goal  is  not  to  create  an  amplitude  correction  factor  similar  to  the  co- 
polar  channel  compensation,  but  to  isolate  the  distortion  terms  of  the  antenna’s 
polarization  distortion  matrix. 

The  mesh  of  long,  thin  cylinders  (wires)  also  gives  this  object  exceptional  bistatic 
scattering  characteristics.  The  wide  main  lobe  of  the  RCS  response  of  each  polarization 
channel  allows  it  to  be  used  as  a  good  polarizing/depolarizing  object  for  bistatic  angles 
wider  than  2°.  This  cannot  be  said  of  the  dihedral,  nor  any  other  object  with  a  high, 
predictable  cross-polarized  RCS.  In  the  Type-2  calibrations  performed  in  the  EMSL,  the 
calibration  objects  are  measured  at  a  bistatic  angle  of  5°,  and  the  resulting  cross¬ 
polarization  isolation  is  quite  good,  as  will  be  demonstrated  in  the  next  section. 
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5.9  Initial  Calibration  Analysis 


For  each  of  the  four  calibration  techniques  described  in  Section  3,  an  initial 
implementation  was  done  as  a  validation  of  the  calibration  routines  and  measurements. 


Figure  32:  Metallic  sphere  calibrated  using  basic  Type-1  calibration  technique 
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Metallic  Sphere:  quasi-monostatic,  (type-2  cal) 
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Figure  33:  Metallic  sphere  calibrated  using  simple  polarimetric  calibration 

technique 


Metallic  Sphere:  quasi-monostatic,  (BICOMS  calibration) 


Frequency  (GHz) 

Figure  34:  Metallic  sphere  calibrated  using  bistatic  via  monostatic  calibration 

technique 
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Metallic  Sphere:  quasi-monostatic,  (type-3  cal) 


Frequency  (GHz) 


Frequency  (GHz) 

Figure  35:  Metallic  sphere  calibrated  using  the  general  dual-antenna  calibration 

technique 

The  calibrated  measurements  of  in  Figure  32  -  Figure  35  show  unacceptably  high 
rippling  in  the  region  of  6-7  GHz.  This  is  due  to  tight  (software)  range  gating  of  the 
measurements.  When  the  software  range  gate  that  is  applied  to  the  calibration  is  of  a 
different  size  or  phase  than  the  gate  performed  on  the  measurement  itself,  this  rippling 
becomes  prevalent.  Provided  that  the  gates  are  consistent  with  each  other  in  every 
measurement,  this  effect  can  be  significantly  reduced.  This  initial  validation  proves  each 
of  the  calibration  routines  (implemented  in  MATLAB)  to  be  legitimate. 

5.10  Basic  Type-1  Calibration  Technique 

This  simple  technique  is  analyzed  mainly  as  a  reference  for  the  calibration  methods 
under  comparison.  An  error  analysis  on  the  Type-1  technique  is  another  term  for  an  error 
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analysis  of  the  calibration  object  prediction  itself,  which  one  must  find  out  in  order  to 
separate  the  shortcomings  of  the  calibration  and  the  shortcomings  of  the  object 
prediction.  The  test  object  on  which  the  calibration  is  used  for  the  analysis  here  and  in 
the  following  sections  is  a  bistatic  (5  degrees  antenna  separation)  of  the  metallic  sphere. 
These  statistics  are  calculated  only  on  the  highest  1  GHz  band  of  the  measurement 
bandwidth  (13-14  GHz)  in  order  to  highlight  the  effect  of  misalignment  on  higher 
frequencies. 


Table  2:  Average  alignment  error  statistics  for  Type-1  calibration  (small  disk) 


Object  Misalignment 

W-error 

HH-error 

VH  level 

HV  level 

or 

0.17  ±0.01  dB 

0.10  +  0.01  dB 

-inf 

-inf 

i° 

0.32  ±  0.04  dB  ; 

0.35+0.01  dB 

-inf 

-inf 

2° 

1.24  ±  0.04  dB 

1.28  +  0.01  dB 

-inf 

-inf 

Upon  implementing  a  Type-1  calibration  with  the  small  disk,  the  statistics  of  Table 
2  result.  The  cross-polar  levels  are  zero  (-infinity  in  dBsm)  because  the  cross-polar  RCS 
of  the  disk  prediction  was  approximated  as  zero. 


By  implementing  the  Type-1  calibration  with  a  sphere,  one  eliminates  the  need  for 
precise  alignment,  at  the  cost  of  average  error  in  the  calibration,  as  seen  in  Table  3. 


Table  3:  Error  statistics  for  Type-1  calibration  (metallic  sphere) 


Object  Misalignment 

W-error 

HH-error 

VH  level 

HV  level 

N/A 

0.79  ±  0.15  dB 

0.28  ±  0.03  dB 

-inf 

-inf 

In  order  to  obtain  reasonable  values  for  the  cross-po] 

larized  calibrated  RCS,  an 

object  must  be  used  that  has  high  RCS  for  all  four  polarization  channels.  The  dihedral, 
tilted  22.5°  from  the  vertical  orientation  is  such  an  object. 
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Table  4:  Error  statistics  for  Type-1  calibration  (dihedral  tilted  22.5°) 


Object  Misalignment 

W-error 

HH-error 

VH  level 

HV  level 

0° 

0.38  ±  0.04  dB 

0.95  ±  0.01  dB 

-39.3  dBsm 

-39.3  dBsm 

1° 

2.76  ±  0.10  dB 

3.68  ±  0.06  dB 

-36.8  dBsm 

-36.0  dBsm 

2° 

7.10  + 0.13  dB 

8.46  ±  0.43  dB 

-32.5  dBsm 

-31.5  dBsm 

As  expected,  a  misalignment  in  the  tilted  dihedral  has  severe  consequences  in  the 
calibration.  However,  this  calibration  object  does  calibrate  all  four  polarization  channels 
rather  well,  offering  a  cross-polarization  isolation  of  about  28  dB  with  a  nominal  0° 
misalignment.  This  value  can  be  obtained  approximately  by  taking  the  co-polar  RCS  of 
the  sphere  in  the  region  of  13-14  GHz  (about  -11  dBsm)  and  subtracting  the  cross- 
polarized  levels  (about  -39  dBsm).  It  is  worthy  to  note  that  the  HH-polarization  error  in 
the  calibration  of  Table  4  is  consistently  higher  than  the  W-polarization  error.  This  is 
most  likely  due  to  errors  in  the  MoM  simulation  of  the  dihedral. 

5.11  EMSL  Simple  Polarimetric  Calibration  Technique 

The  first  object  set  used  for  the  analysis  of  this  calibration  is  the  typical  set  used  in 
the  EMSL  procedure.  For  the  co-polar  amplitude  reference  object  (Object  1),  the  small 
disk  is  used.  For  the  cross-polar  reference  objects  (Objects  2  and  3,  respectively),  the 
wire  mesh  is  used  in  the  vertical  orientation  and  subsequently  with  the  wires  aligned  45° 
counter-clockwise  from  vertical.  The  next  table  displays  the  results  of  the  calibration 
error  analysis  for  this  object  set. 
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Table  5:  Error  statistics  for  EMSL  calibration  (Tl:  small  disk,  T2:  vertical  wire 

mesh,  T3:  tilted  wire  mesh) 


Object  Misalignment 

W-error 

HH-error 

VH  level 

HV  level 

Tl=(f,  T2=(f,  T3=(f 

0.17  ±  0.01  dB 

-42.0  dBsm 

-42.0  dBsm 

Tl=l°,  T2=l°,  T3-10 

0.32  ±  0.04  dB 

0.35  ±0.01  dB 

-41.4  dBsm 

-41.4  dBsm 

Tl  =2°,  T2=0°,  T3=& 

1.24  ±  0.04  dB 

1.2810.01  dB 

-40.8  dBsm 

-40.8  dBsm 

Tl-(f,  T2=2°,  T3=(f 

0.17  +  0.01  dB 

0.1010.01  dB 

-43.7  dBsm 

-43.7  dBsm 

T1=0P,  T2-(f,  T3=2° 

0.17  ±0.01  dB 

0.1010.01  dB 

-39.7  dBsm 

-39.7  dBsm 

Tl=(f,  T2-20,  T3=2° 

0.1710.01  dB 

0.1010.01  dB 

-41.3  dBsm 

-41.4  dBsm 

Tl=2°,  T2=2°,  T3-T 

1.24  + 0.04  dB 

1.2810.01  dB 

-40.1  dBsm 

-40.2  dBsm 

For  the  “perfect”  alignment  of  the  calibration  objects,  the  Simple  Polarimetric 
calibration  offers  a  cross-polarization  isolation  of  about  3  dB  more  than  the  Type-1 
calibration  shown  in  Table  4  (-42  dBsm  compared  with  -39  dBsm  cross-polarization 
RCS  levels). 

For  the  Type-2  calibration  as  performed  in  the  EMSL,  the  object  set  used  in  Table  5 
is  the  only  appropriate  set  that  can  be  used.  The  wire  mesh  is  the  only  known  object  to 
return  high  cross-polarized  RCS  that  is  strongly  correlated  with  the  co-polar  RCS 
(approximated  by  a  cos2  relationship)  at  a  bistatic  angle  of  more  than  2°.  The  small  disk, 
though  very  sensitive  to  alignment,  returns  the  largest  bistatic  co-polar  RCS  of  any 
object,  which  reduces  the  signal-to-noise  ratio  in  the  calibration.  The  EMSL  technique 
could  be  modified,  however,  to  use  the  tilted  dihedral  as  a  cross-polarization  reference 
object,  if  the  calibration  routine  is  adjusted  to  use  a  MoM  prediction  rather  than  an 
approximate  prediction  based  on  the  vertical  wire  mesh.  This  would  eliminate  one  object 
in  the  calibration — only  the  small  disk  and  the  tilted  dihedral  (in  a  single  orientation) 
would  need  to  be  used. 
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Table  6:  Error  statistics  for  modified  EMSL  calibration  (Tl:  small  disk,  T2: 

dihedral  tilted  22.5°) 


Object  Misalignment 

W-error 

HH-error 

VH  level 

HV  level 

Tl=(f,  T2=0° 

0.17  +  0.01  dB 

0.10  +  0.01  dB 

-39.3  dBsm 

-38.7  dBsm 

ta 

ii 

***** 
w  © 

Kj 

II 

o 

0.32  ±  0.04  dB 

0.35  ±0.01  dB 

-36.8  dBsm 

-36.0  dBsm 

Tl=2°,  T2=0° 

1.24  ±  0.04  dB 

1.28  ±0.01  dB 

-39.3  dBsm 

-38.7  dBsm 

Tl=(f,  T2-2° 

0.17  ±0.01  dB 

0.10  ±0.01  dB 

-32.5  dBsm 

-31.5  dBsm 

TJ=2°,  T2=2° 

1.24  ±  0.04  dB 

1.28  ±0.01  dB 

-40.1  dBsm 

-40.2  dBsm 

,1  ii  A  f* 

As  seen  from  Table  6,  the  cross-polarization  isolation  is  about  3  dB  less  than  that  of 
the  normal  EMSL  calibration  (Table  5),  and  the  co-polarized  RCS  error  is  identical.  This 


degradation  is  weighed  against  the  advantage  of  measuring  one  less  calibration  object.  In 
effect,  this  modified  Type-2  calibration  is  identical  in  procedure  to  a  Type-1  calibration, 
with  the  exception  that  two  objects  are  used:  one  for  co-polar  amplitude  calibration,  and 
the  other  for  cross-polar  amplitude  calibration.  Comparing  this  calibration  with  the 
Type-1  using  the  tilted  dihedral  (Table  4)  one  observes  that  the  use  of  the  disk  enhances 
the  co-polarized  RCS  error,  while  the  cross-polarization  isolation  is  the  same. 


5.12  Generalized  Dual-Antenna  Calibration  Technique 

This  technique  (abbreviated  as  the  GCT  in  (Whitt,  et  al,  1991),)  was  expected  to 
yield  the  best  results  for  ideal  alignment  conditions  and  given  accurate  theoretical 
solutions.  The  effect  of  misalignment  of  one  or  more  of  the  reference  objects,  however, 
was  unknown.  Figure  36  displays  the  co-polarized  error  in  a  measurement  of  the  metallic 
sphere  calibrated  with  the  GCT,  using  what  is  considered  to  be  a  near-optimal  object 
set— the  three  objects  used  were  the  small  disk  (object  1),  the  vertical  dihedral  (object  2), 
and  the  dihedral  tilted  at  22.5°  from  vertical  (object  3).  The  scattering  matrices  for  each 
of  these  objects  have  the  form: 
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(58) 


T,= 


'1  o' 

(56) 

T 2  = 

'-1  O' 

(57) 

T3  = 

'-1  -1' 

0  1 

_  0  1_ 

1  1  _ 

Every  scattering  matrix  meets  the  requirement  of  mutual  linear  independence. 


Co-polar  error  in  GCT  calibration  of  sphere:  "best"  alignment 


Figure  36:  Error  in  GCT  calibration  using  small  disk,  vertical,  and  tilted  dihedral 

For  this  calibration,  the  theoretical  prediction  for  the  small  disk  was  approximated 

with  PO,  and  the  predictions  for  the  dihedral  was  generated  using  an  MoM  code,  for 
which  the  scattering  geometry  in  the  simulation  reflected  the  true  geometry  of  the  dual¬ 
sled  measurement  in  the  EMSL. 

Using  the  same  object  set,  but  with  a  systematic  misalignment  of  2°,  the  following 
error  versus  frequency  plot  was  generated: 
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Co-po)ar  error  in  GCT  calibration  of  sphere:  2  deg  alignment 


Figure  37:  Error  in  GCT  calibration  with  2°  misalignment  associated  with  each 

reference  object 

As  one  might  expect,  when  misalignment  is  introduced  into  the  calibration,  the 
calibrated  measurement  will  tend  to  overestimate  the  true  prediction,  and  the 
overestimation  will  increase  with  frequency.  Table  7  displays  the  mean  and  variance  of 
the  RCS  error  between  the  calibrated  measurements  and  the  prediction  of  the  sphere. 
These  statistics  are  performed  only  on  the  highest  1  GHz  band  of  the  measurement 
bandwidth  (13-14  GHz)  in  order  to  highlight  the  effect  of  misalignment  on  higher 
frequencies. 
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Table  7:  Average  alignment  error  statistics  for  GCT  calibration  (Tl:  small  disk, 
T2:  vertical  dihedral,  T3:  tilted  dihedral) 


Object  Misalignment 

W-error 

HH-error 

VH  level 

HV  level 

77=0°,  T2=(f,  T3=(f 

0.17  ±0.01  dB 

0.10  ±0.03  dB 

-49.2  dBsm 

-51.5  dBsm 

Tl=l°,  T2=l°,  T3=l° 

0.32  ±  0.04  dB 

0.35  ±  0.01  dB 

-48.8  dBsm 

-51.5  dBsm 

Tl  =2°,  T2=(f,  T3—(f 

1.24  ±  0.04  dB 

1.28  ±0.01  dB 

-47.8  dBsm 

-50.6  dBsm 

T1=0°,  T2-2°,  T3=(f 

0.17  ±0.01  dB 

0.10  ±0.01  dB 

-51.4  dBsm 

Tl-0°,  T2=(f,  T3=2° 

0.17  ±0.01  dB 

0.10  ±0.01  dB 

T1=0°,  T2=2°,  T3=2° 

0.17  ±0.01  dB 

0.10  ±0.01  dB 

-49.2  dBsm 

Tl=2°,  T2=2°,  T3=2° 

1.24  ±  0.04  dB 

1.27  ±0.01  dB 

-47.8  dBsm 

-50.5  dBsm 

The  cross-polarization  isolation  for  this  case  is  10  d] 

B  greater  than  the  Type-1 

calibration  (Table  4)  and  about  7  dB  greater  than  the  Simple  Polarimetric  Calibration 
(Table  5).  If  the  object  set  is  re-ordered,  so  that  the  vertical  dihedral  is  now  Tj,  the  tilted 
dihedral  is  T2,  and  the  small  disk  is  T3,  the  following  error  statistics  result: 


Table  8:  Error  statistics  for  GCT  calibration  (Tl:  tilted  dihedral,  T2:  small  disk, 

T3:  vertical  dihedral) 


Object  Misalignment 

W-error 

HH-error 

VH  level 

HV  level 

T1=0P,  T2=0°,  T3=(f 

0.24  ±  0.02  dB 

-0.06  ±  0.02  dB 

-33.9  dBsm 

-51.6  dBsm 

Tl=l°,  T2=l°,  T3=l° 

-0.44  ±  0.03  dB 

-0.73  ±  0.03  dB 

-30.6  dBsm 

-51.3  dBsm 

Tl=2°,  T2=(f,  T3=(f 

0.24  ±  0.02  dB 

-0.06  ±  0.02  dB 

-33.7  dBsm 

-51.9  dBsm 

Tl-(f,  T2=2°,  T3=(f 

-2.76  ±  0.06  dB 

-3.05  ±  0.06  dB 

-28.8  dBsm 

-48.6  dBsm 

T1=0°,  T2=0°,  T3=2° 

0.22  ±  0.03  dB 

-0.08  ±  0.02  dB 

-28.3  dBsm 

-5 1 .4  dBsm 

Tl=2°,  T2=(f,  T3=2° 

0.22  ±  0.03  dB 

-0.08  ±  0.02  dB 

-28.4  dBsm 

Tl  =2°,  T2=2°,  T3=2° 

-2.77  ±  0.07  dB 

-3.07  ±  0.07  dB 

-25.1  dBsm 

Jj 

The  cross-polar  levels  are  (specifically,  of  the  VH-channel)  are  considerably 
lower  for  the  calibration  of  Table  7.  Realizing  this,  as  well  as  observing  that  the  penalty 
resulting  from  the  misalignment  of  T2  in  Table  8  is  more  than  twice  the  penalty  of  the 
misalignment  of  Ti  in  Table  7,  it  is  apparent  that  though  the  objects  sets  are  identical,  the 
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order  in  which  the  objects  are  implemented  in  the  calibration  is  important.  Only  Object  1 
is  used  in  the  absolute  amplitude  calibration  in  the  GOT  method.  Therefore  for  any  given 
object  set,  the  most  accurate  calibration  will  be  that  one  which  uses  the  object  that  has  the 
most  precise  theoretical  prediction  and  uses  it  as  Object  1.  In  this  case,  the  most  precise 
prediction  accompanied  the  small  disk,  and  as  a  result,  the  first  calibration  (Table  7)  used 
the  disk  as  Object  1  and,  as  a  result,  was  most  accurate. 

For  the  measurement  of  the  metallic  sphere,  good  cross-polarization  isolation  will 
be  demonstrated  in  very  low  cross-polarized  RCS  levels.  The  first  calibration  represents 
the  better  cross-polarization  isolation.  In  order  to  optimize  the  cross-polar  isolation  in  the 
measurement,  one  must  determine  which  of  objects  2  and  3  have  the  most  precise 
prediction,  and  which  of  Eqn’s  (36)  and  (39)  should  be  used  in  the  calibration  routine  to 
offer  the  best  cross-polarization  isolation. 

5.13  BICOMS  Full-Polarimetric  Calibration  Technique 

The  BICOMS  technique  is  a  modified  implementation  of  the  GCT.  This  method 
solves  for  the  monostatic  subsystem  distortion  characteristics  of  each  leg  of  the  bistatic 
measurement  path,  and  uses  this  complete  distortion  model  to  solve  for  the  distortion 
along  the  whole  bistatic  path.  Since  the  PDM’s  of  the  antennae  are  solved  for  using 
monostatic  measurements,  it  is  not  necessary  to  address  the  difficult  problem  of  finding  a 
bistatic  cross-polarization  calibration  reference  object.  This  also  means,  though,  that 
more  calibration  measurements  must  be  taken,  and  the  operational  efficiency  of  the 
measurement  facility  is  affected. 
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Four  objects  are  used  for  the  analysis  of  this  calibration.  In  the  actual  BICOMS 
procedure,  a  short  cylinder  is  used  as  one  of  the  three  calibration  reference  objects  that 
determine  the  system  distortion,  and  the  same  cylinder  is  also  rotated  90°,  and  the  end- 
cap  of  the  cylinder  is  measured  to  obtain  an  absolute  amplitude/phase  correction.  The 
remaining  two  calibration  objects  are  a  horizontal  dihedral,  and  a  dihedral  oriented  67.5° 
from  horizontal.  In  the  interest  of  consistency  with  the  calibrations  of  the  previous 
section,  this  analysis  uses  a  vertical  dihedral  and  a  dihedral  oriented  22.5°  from  vertical. 
This  does  not  have  any  bearing  on  the  performance  of  the  calibration.  Table  9  lists  the 
error  statistics  for  this  particular  calibration.  The  cross-polarization  isolation  achieved 
with  this  calibration  is  similar  to  that  achieved  in  the  optimum  GCT  calibration  (Table  7). 


Table  9:  Error  statistics  for  BICOMS  calibration  (TO:  small  disk,  Tl:  short 
cylinder,  T2:  tilted  dihedral,  T3:  vertical  dihedral) 


Object  Misalignment 

W-error 

HH-error 

VH  level 

HV  level 

T0=(f,  Tl—(f,  T2=(f, 
T3=(f 

0.12  ±  0.03  dB 

-0.06  ±  0.01  dB 

-47.0  dBsm 

-49.9  dBsm 

T0—1°,  Tl—(f,  T2=(f, 
T3=0° 

-0.17  ±  0.04  dB 

-0.34  ±  0.01  dB 

-46.7  dBsm 

-49.6  dBsm 

T0=2°,  Tl-(f,  T2=(f, 
T3-(f 

-1.09  ±  0.04  dB 

-1.26  ±0.01  dB 

-45.8  dBsm 

-48.7  dBsm 

T0=(f,  Tl-1°,  T2=l°, 
T3=l° 

0.08  ±  0.02  dB 

-0.06  ±  0.01  dB 

-45.7  dBsm 

-51.0  dBsm 

T0=(f ,  Tl  =2°,  T2=2°, 
T3=2° 

0.17  + 0.02  dB 

-0.06  ±0.01  dB 

-45.7  dBsm 

-50.6  dBsm 

Objects  T1-T3  are  all  monostatic  measurements  that  cl 

laracterize  the  polarimetric 

distortion  involved  in  each  leg  of  the  bistatic  path.  Since  objects  T1-T3  are  used  only  to 
solve  for  the  PDM’s  and  not  the  absolute  amplitude  of  the  calibrated  data,  they  need  not 
necessarily  be  accompanied  by  a  precise  MoM  prediction.  If  the  scattering  matrix 
structure  of  these  objects  is  known  (which  it  is,  according  to  Eqn’s  (56)  -  (58),  these 
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simple  binary  matrices  can  be  used  in  place  of  the  theoretical  prediction.  The  only  object 
that  needs  an  exact  prediction  is  the  amplitude  reference  object  (TO).  The  calibration 
implementation  shown  in  Table  9  was  performed  in  this  manner  with  near-identical 
results. 
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6  CONCLUSIONS 


6.1  Calibration  Object  Selection 

It  has  been  shown  that  the  sensitivity  of  a  particular  calibration  object  to  alignment 
error,  though  significantly  affecting  the  error  in  the  absolute  amplitude  calibration,  has  a 
less  pronounced  effect  in  the  determination  of  the  associated  antenna  PDM’s.  This  has 
been  demonstrated  with  the  wire  mesh  and  dihedral  specifically,  where  a  2°  misalignment 
in  these  objects  affects  a  loss  in  cross-polarization  isolation  of  about  2  dB. 

For  full-polarimetric  (Type-3)  calibration,  it  is  also  clear  that  selecting  an  object 
with  a  simple  scattering  matrix  offers  great  advantages  in  the  way  of  computational 
efficiency  and  accuracy.  Binary  scattering  matrix  structures,  such  as  the  vertical 
dihedral,  sphere,  and  disk  (among  many  others)  offer  the  capability  to  predict  the  antenna 
PDM’s  more  precisely  than  if  these  objects  were  predicted  with  a  complex  computational 
method,  such  as  MoM.  In  the  amplitude  correction  of  a  Type-3  calibration,  the  same 
considerations  apply  as  to  that  of  a  Type-1  calibration:  namely  that  objects  with  high 
RCS  and  are  simultaneously  insensitive  to  alignment  error  are  optimal. 

It  has  further  been  shown  that  the  parallel  wire  mesh  offers  promise  as  a  versatile 
bistatic  cross-polarization  calibration  reference  object.  It  offers  high  cross-polarization 
RCS  at  a  much  greater  bistatic  angular  range  than  any  other  object  evaluated  here. 
Further  investigation  into  the  suitability  of  this  object  for  a  polarimetric  calibration  would 
include: 

•  Finding  accurate  predictions  for  the  mesh  using  a  numerical  electromagnetic  code 
(NEC), 
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•  Exploring  the  optimum  structure  of  the  mesh,  considering  variables  such  as  wire 
length  and  spacing,  and  selection  of  the  foam  mounting  substrate, 

•  The  bistatic  angular  range  at  which  the  wire  mesh  gives  good  cross-polarization 
isolation  in  a  Type-2  or  Type-3  calibration, 

•  The  effect  of  segmenting  the  wire  mesh  in  order  to  mitigate  traveling  wave  and 
grating  lobe  effects  in  the  RCS  return. 

6.2  Performance  of  Type-1,  Type-2,  and  Type-3  Calibrations 

For  the  measurement  parameters  of  the  EMSL  (antenna  polarization  purity,  gain, 
signal-to-noise  ratio,  etc),  it  has  been  shown  that  the  gain  in  cross-polarization  isolation 
in  using  a  Type-3  technique  as  opposed  to  a  Type-1  technique  is  about  10  dB.  The 
absolute  error  in  the  co-polar  channels  is  comparable  between  comparisons,  provided  that 
the  alignment  error  and  signal-to-noise  ratio  are  kept  relatively  constant.  In  general,  for 
any  type  of  calibration,  the  sources  of  absolute  amplitude  error  and  degradation  of  cross¬ 
polarization  isolation  are  separate.  Each  calibration  technique  evaluated  here  follows  the 
general  paradigm  that  two  objects  with  highly  independent  scattering  matrices  are 
necessary  to  characterize  polarization  distortion,  and  a  single  object  with  high  co-polar 
RCS  is  necessary  to  accurately  determine  the  amplitude  and  phase  constant  applied  to 
these  PDMs. 

Given  any  particular  calibration  object,  such  as  those  listed  in  Table  1,  it  is  not 
difficult  to  predict  the  effect  of  misalignment  on  the  overall  error  of  the  calibration. 
There  is  a  1:1  correspondence  between  the  dB  error  per  degree  of  misalignment,  and  the 
average  amplitude  error  of  the  co-polar  channels  in  the  calibration. 

These  results  offer  fair  insight  into  the  capabilities  of  these  various  types  of 
calibrations  relative  to  each  other.  However,  the  results  are  only  applicable  in  the  context 


87 


of  the  EMSL  measurement  environment.  Further  study  into  this  subject  should  include  a 
similar  analysis  in  varying  measurement  environments. 
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APPENDIX  A:  VALIDATION  OF  MBETs  AND  MoM 
COMPUTATIONS 

In  order  to  experimentally  validate  the  accuracy  of  the  MoM  computations 
performed  on  the  Fast-Illinois  Solver  Code  (FISC),  a  single-frequency  study  was  done  for 
two  objects:  one  simple  object  and  another  rigorously  complex  object,  consisting  of  an 
assembly  of  canonical  shapes,  including  a  comer  reflector,  open  and  closed  cylinders, 
and  a  shadowing  plate.  In  addition,  Falconer’s  monostatic-to-bistatic  scattering 
conversion  formulae  for  far- field  and  near-field,  (Falconer,  “Extrapolation  of...”,  1988) 
and  (Falconer,  “Near-Field  Statement...”,  1988),  are  evaluated  for  their  ability  to  predict 
bistatic  scattering  phenomena. 


6.3  MBET  Theory 

Kell’s  commonly  cited  MBET  (Kell,  1965)  relates  the  bistatic  RCS  of  a  point 
scatterer  to  it’s  monostatic  RCS  at  the  bisector  of  the  bistatic  angle  and  at  a  reduced 
frequency.  Kell’s  MBET  is  stated  as: 


<*b&,OrJ)= 


(59) 


where 

<jb  -  Estimated  Bistatic  RCS  in  dBsm 
<jM  -  Monostatic  RCS  in  dBsm 
0=  Equivalent  Monostatic  Azimuth  Angle 
Or  =  Bistatic  Receiver  Azimuth  Angle 
0r=  Bistatic  Transmitter  Azimuth  Angle 
/  =  Transmitter  Frequency 
fM  =  Equivalent  Monostatic  Frequency 
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Based  on  the  Physical  Optics  (PO)  model  for  bistatic  scattering  of  point  scatterers  at 
small  bistatic  angles,  the  equivalent  monostatic  illumination  angle  is  approximated  by  the 
bisector  of  the  bistatic  angle  of  the  measurement. 

By  modeling  the  object  as  a  set  of  discrete  scattering  centers,  Kell  states  that  if  for 
any  given  aspect  angle,  the  amplitude  and  phase  of  the  scattering  centers  are  insensitive 
to  the  bistatic  angle  jB  over  the  range  of  interest,  then  the  monostatic  cross-section 
measured  at  the  bisector  of  the  bistatic  angle  and  at  a  frequency  lower  by  a  factor  of 
cos(/?/2 )  is  equal  to  the  bistatic  cross-section. 

The  range  of  angles  at  which  this  assumption  is  valid  is  dependent  upon  the  object’s 
geometry.  Classes  of  objects  for  which  specular  reflection  is  the  dominant  scattering 
mechanism  are  best  suited  to  this  approximation.  The  maximum  bistatic  angle  for  which 
the  approximation  is  valid  is  also  highly  dependent  on  the  beamwidth  of  the  pattern  of  the 
dominant  scattering  centers. 


Another  commonly  cited  MBET,  developed  by  Crispin  (Crispin  and  Siegel,  1968), 
is  similar  to  Kell’s,  but  simplified  to  exclude  the  adjustment  of  the  monostatic  equivalent 
frequency. 


(60) 


The  advantage  of  this  expression  of  the  MBET  can  be  immediately  seen  in  that  it 
requires  no  RCS  measurements  at  frequencies  different  from  the  bistatic  measurement. 

Falconer’s  adaptation  to  the  above  theorems  arrives  at  a  different  relationship  between 
the  bistatic  angle  and  equivalent  monostatic  aspect  angle.  Falconer’s  theorem  states: 
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The  above  MBET  is  actually  valid  in  both  the  near  and  far-fields.  Falconer's 
MBET,  similar  to  Crispin's,  approximates  the  bistatic  return  by  using  the  monostatic 
signal  at  the  same  transmission  frequency  but  at  a  reduced  look  angle.  While  the 
structure  of  the  formula  appears  similar  to  Crispin’s  formulation  of  the  MBET,  Falconer 
employed  the  discrete  scattering  center  model  of  the  object’s  radiation  pattern  for  his 
formulation,  the  same  formulation  that  resulted  in  Kell’s  MBET. 

This  work  focuses  on  a  near-field  MBET  derived  by  Falconer  and  based  on  a 
similar  theorem  stated  by  Kell.  To  arrive  at  the  MBET,  the  vector  potential  is  written  in 
scalar  form  by  approximating  the  scattering  surface  to  be  planar  and  the  angles  defining 
the  bistatic  path  to  be  paraxial,  i.e. 

A(xr )  =  J  J(x)G(x,  xr  )dx  (62) 

where  A(xr )  is  the  vector  potential  observed  at  the  point  xr ,  induced  by  the  current  J(x) , 
and  G(x,xr)  is  the  Green’s  function  for  the  particular  medium  and  frequency. 

The  magnetic  field  incident  upon  the  surface  can  be  written  as 

H(x)  =  H0elir’  (63) 

where  r'-\x-x\  and  x,  is  the  location  of  the  transmitting  antenna. 

The  Physical  Optics  current  is 

J(x)  =  2HJ(x)e,kr'  (64) 
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where  T(x)  is  the  complex  reflectivity  of  the  scatterer,  assumed  to  be  relatively 
independent  of  frequency. 

This  is  then  substituted  into  the  vector  potential  yielding 


A(xr)  =  2H0jT(x)eltrG(x,xr)dx 


(65) 


The  vector  potential  is  related  to  the  E-field  by 


E  -  -\20mkA 


(66) 


in  free  space.  The  E-field  at  the  observation  point  then  becomes 


4/rr 


(67) 


which  for  the  two  dimensional  case  becomes 


E{x,)  =  Ea\ 


•y/-  il  XrT(x)eil 
4ttt 


~dx 


(68) 


The  Green’s  function  is  approximated  with  a  Fresnel  phase  approximation  for  each 
of  the  near-field  and  far-field  cases.  For  the  near-field  (valid  when 
0.6(£>3  IX)'11  <R< 2D1  / A )  is 


r  &  R-x sin<?r  -zcosOr  +(x2  +z2)/2R 


(69) 


and  for  the  far-field  (when  R  >  2D 2 IX), 

r  aR-xsmOr -zcos0r  (70) 

The  above  expressions  are  inserted  into  the  Green’s  function  exp(dr/4w)  in  the 
numerator  of  the  exponent,  and  the  denominator  is  approximated  by  4m-  ~  4nR . 

The  scalar  near-field  scattering  amplitude  can  now  be  expressed  as 
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(71) 


Eb{dr,8],k,R)  =  E^-HMenkR\  T(x,z ) 

•  exp[f'A:(x2  +z2)/ R-ikx(sindr  +  sin  0, ) 

-  ikz(  cos  9r  +  cos  0,  )]dx 

In  the  development  of  the  near-field  MBET  (heretofore  designated  as  F2)  from  the 
above  expression,  the  exponent  is  simplified  using  trigonometric  identities,  and  the 
reduced  monostatic  frequency  defined  by  Kell  (k'=  kcos0U2 ),  where  Om  is  half  the  bistatic 
angle.  In  addition,  Falconer  introduces  a  reduced  range  R'=Rcos0l/2.  The  monostatic 
amplitude  is  observed  at  these  reduced  frequency  and  range  parameters  that  altogether 
yields  the  scattering  amplitude  of  the  following  form,  where  0  =  0r=0,: 

, _  (72) 

T(x>z) 

•  expH'£'  (x2  +z2)/  R-i2k'xsm.6- i2k'zcos0]dx 


Recognizing  that  X'R'=  XR  and  k'/R'=k/R,  we  get 

Em(6,k',R')  =  Ej-i/XRenk'R  \  T(x,z ) 

•  exppAfx2  +  z2)!  R -ilk'  x{smO  +  cos  0)]dx 


(73) 
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A  comparison  of  this  monostatic  scattering  amplitude  to  the  bistatic  scattering 
amplitude  reveals  a  relationship  defined  by: 

Eb  (0r,0„k,R)  =  enk*-ik'R'Em  (< 9a;k',R' )  (74) 

where  0a  is  the  bisector  of  the  bistatic  angle  (0a  =  (9r  +  0,)/2).  This  is  the  expression  for 

Falconer’s  near-field  MBET,  which  requires  adjustment  of  the  monostatic  frequency  and 
range  in  order  to  approximate  the  bistatic  cross-section. 


Figure  38:  Falconer’s  MBET  geometry.  R  is  observation  sphere  radius,  a  and  p 
observation  angles,  both  presumed  small,  D  represents  object’s  nominal  diameter, 
and  H  its  nominal  height  relative  to  the  illumination  direction.  (Falconer, 

“Extrapolation  of...”,  1988) 

Substituting  the  appropriate  expressions  for  k’  and  R  \  Falconer’s  extended  near¬ 
field  MBET  (F2)  is  stated  completely  in  terms  of  frequency  and  range  by: 


E,(er,eR,f,k)=e 


(75) 


where 
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,  Monostatic  Angle 


„  (eT+e' \ 

o=  — — - 

l  J 

fu  =  /  cos^^-j ,  Reduced  frequency 
Rm  =  /?cos(^r~^j  ,  Reduced  Range 

Falconer’s  CW  version  of  the  F2  MBET  is  valid  for  both  near  field  and  far  field 
scenarios.  A  different  angle  parameter,  df  is  introduced  which  is  defined  as 

.  „  sin  9  +  sin  9,  (76) 

sin  9 ,  = - - - - 

f  2 

and  the  original  relationship  between  the  bistatic  and  monostatic  scattered  fields  is 


Eb(9,AXR)  =  E0J^illRe™\  T(x,z ) 

•  exp[i/:(x2  +  z 2 )  /  R  -  ilkx  sin  6f 
-  ilkz{cos  0r  +  cos  0t  )}dx 


(77) 


If  dr  and  d,  are  small,  then 

cos^7  «cos  6t  «  cos  0j  =1 


(78) 


and  it  follows  that 

Et{6r,e,\k,R)  =  Em{fi/XR)  (79) 

Therefore  the  frequency  and  range  of  the  bistatic  measurement  are  the  same  as  the 
monostatic  equivalent  for  small  bistatic  angles.  We  designate  this  MBET  as  FI . 

The  region  of  applicability  of  Falconer’s  CW  MBET  is  restricted  by  the  object’s 
nominal  thickness  at  the  desired  observation  angle.  As  seen  in  Figure  38,  the  thickness  H 
is  defined  as  the  linear  extent  of  the  object  in  the  direction  of  the  observation  angle.  It  is 

appropriate  to  apply  the  CW  MBET  when  the  bistatic  angle  satisfies  db  <  sjxiH  . 
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6.4  Methodology 


The  following  presents  a  five  phased  research  project  consisting  of:  1)  the 
collection  of  far-field  empirical  data,  2)  the  evaluation  of  Falconer’s  far-field  MBET  with 
the  empirical  data,  3)  the  validation  of  the  computer  prediction  code  with  empirical  data, 
4)  the  generation  of  near-field  and  far-field  scattering  with  the  computer  code,  and  5)  the 
evaluation  of  both  Falconer’s  MBETs  with  computer  generated  scattering  data.  Objects 
A  and  B,  shown  in  Figure  2,  were  chosen  for  their  ability  to  highlight  various  scattering 
phenomena  at  different  look  angles,  such  as  traveling  waves,  creeping  waves,  multi¬ 
bounce,  and  shadowing  effects.  Co-polarized  (W,  HH)  monostatic  and  bistatic  waterline 
pattern  cut  and  imaging  data  was  collected  according  to  the  matrix  in  Table  1  at  the 
European  Microwave  Signature  Laboratory  (EMSL)  of  the  EC  Joint  Research  Center  in 
Ispra,  Italy. 


Object  A 


280.3 


HH* 
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For  the  purposes  of  this  chapter,  all  measurements  and  analysis  occur  in  the 
horizontal  (x-y)  plane  defined  relative  to  the  measurement  object.  Specifically,  a 
Cartesian  coordinate  system  is  defined  such  that  the  z-axis  is  tangent  to  the  long  flat  plate 
dimension  of  object  A  and  the  cylinder  axes  of  object  B.  The  x-axis  is  then  defined  to  be 
the  plate  normal  of  objects  A  and  B,  where  the  normal  points  away  from  the  cylinders  in 
object  B.  The  bistatic  arrangement  is  defined  by  the  transmitter  illumination  angle,  0t, 
and  receiver  look  angle,  0r,  where  these  are  restricted  to  the  (x-y)  plane  (i.e.  0  =  90°  in 
spherical  coordinates).  Therefore,  the  bistatic  angle  becomes  0t-0r.  Polarization  is 
defined  by  the  E-field  vector  orientation.  For  example,  VV-polarization  indicates  both 
the  incident  and  scattered  E-field  polarizations  are  z  directed.  Similarly,  HH-polarization 
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indicates  the  E-field  vectors  lie  in  the  (x-y)  plane  and  are  orthogonal  to  the  propagation 
direction.  Only  co-polarized  data  is  analyzed  in  this  report  (i.e.,  same  transmit  and 
receive  polarizations).  No  measured  bistatic  data  exists  for  0r  angles  within  4°  of  the 
0T  angle.  In  the  azimuth  plane,  bistatic  estimated  data  is  generated  from  the  monostatic 
far  field  measurements  through  Falconer’s  CW  monostatic-to-bistatic  equivalence 
theorem  given  by  (Falconer,  “Near-Field  Statement  of. . .”,  1988). 

Simulation  data  are  generated  for  the  far-field  and  near-field  test  set  by 
commercially  available  electromagnetic  field  modeling  code  which  computes  scattering 
from  an  object  using  the  method  of  moments  (MoM). 

6.5  Results 

6. 5. 1  Far-field  MoM  code  and  MBET  evaluation 

As  an  initial  validation  of  the  scattering  code  and  the  MBETs,  a  bistatic 
measurement  of  Object  A  with  the  transmitter  positioned  at  zero  degrees  of  azimuth  and 
elevation  was  compared  with  the  MoM  simulation  and  also  plotted  against  the  MBET 
approximations. 
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Figure  40:  Object  A  Measured  vs.  Predicted  Bistatic  RCS,  W  Polarization,  8  HGz, 

0°  Transit  Angle 


Figure  40  shows  excellent  agreement  between  the  measurements  and  MoM 
simulation  data.  While  all  MBETs  perform  well  up  to  about  a  30°  bistatic  angle,  Kell’s 
and  Falconer’s  CW  MBET  (FI)  show  continued  agreement  out  to  about  55°.  Note  that 
Kell’s  MBET  is  limited  by  the  frequency  of  the  monostatic  data  collected. 


6.5.2  Near-field  MBET  evaluation  (fixed-angle  bistatic  RCS) 

A  primary  region  of  interest  for  the  evaluation  of  the  near-field  MBETs  is  the  center 
of  the  radiating  near-field.  Location  of  objects  in  the  middle  of  their  radiating  near  field, 
defined  in  (Gabig  et  al,  2000),  is  governed  by 


(80) 


where  X  refers  to  the  transmitting  signal's  wavelength,  and  /  designates  the  longest 
dimension  of  the  object  bisected  by  the  waterline  cut.  For  Object  A,  /  =  240mm,  and  for 
Object  B,  /  =  380mm.  We  analyzed  the  objects  at  8  GHz,  which  is  X  =  37.5mm.  The 
variable  R  designates  the  limits  of  the  near-field  radiating  region.  Previously  we 
calculated  the  bistatic  return,  and  its  corresponding  monostatic  return  for  Object  B  at  the 
region  i?’s  midpoint. 

Maintaining  a  constant  angular  separation  between  the  transmitter  and  receiver 
antenna  of  10°  and  rotating  this  system  around  the  object,  we  generated  the  data  in  Figure 
4.  This  data  illustrates  reasonable  null  alignment  between  the  computer  generated  near¬ 
field  bistatic  data  and  Falconer’s  F2  MBET  approximation,  where  the  main  scattering 
mechanism  is  specular  return  from  the  flat  shadowing  plate.  Additionally,  the  MBET 
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continues  to  track  the  nulls  reasonably  well  after  the  transmitter  passes  90°  and  begins 
illuminating  more  complex  components  of  Object  B. 

The  initial  results  obtained  and  illustrated  in  Figure  41  would  tend  to  indicate  that 
the  F2  MBET  performs  well  throughout  a  range  of  observation  angles.  To  determine  a 
more  quantitative  limit  on  the  antenna  separation  that  would  produce  results  that  are 
approximated  well  by  the  MBETs,  we  ran  a  series  of  near-field  MOM  simulations  on 
both  Objects  A  and  B,  for  illumination  angles  of  0°,  45°,  and  135°.  The  simulations  focus 
again  on  the  center  of  the  near-field  for  both  objects.  For  comparison,  simulations  for 
both  objects  were  also  taken  in  the  far-field. 

6.5.3  Near-field  MBET  evaluation  (swept  angle  bistatic  RCS) 

Following  are  figures  which  display  the  RCS  and  phase  errors  between  the  bistatic 


Object  B  8  GHz  HH  BSA  =  10  [Deg] 
BIS  [dBsml]  vs  MBET  [dBsm2] 


Figure  41 :  Object  B  MoM  vs.  F2  MBET  Fixed  Angle 
Bistatic  RCS,  HH  Polarization,  8  GHz,  10°  Bistatic  Angle. 
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“truth”  data  generated  by  MOM  and  the  F2  MBET  approximation  applied  to  monostatic 
data  generated  with  the  same  software.  The  frequency  and  range  parameters  used  in  the 
monostatic  calculations  were  adjusted  corresponding  to  the  appropriate  values  as 
determined  by  the  F2  MBET  and  quantized  to  1 1  levels  (in  order  to  save  computational 
time).  A  sliding  average  window  spanning  10°  was  applied  to  the  error  versus  angle  data 
in  order  to  smooth  the  trend  such  that  a  stable  threshold  might  be  established.  The 
threshold  of  tolerable  error  was  taken  to  be  +/-  1  dB  or  +/-  22.5°,  corresponding  to  the 
error  limits  set  forth  in  the  development  of  the  MBET  itself. 

For  an  incidence  angle  of  zero  degrees,  as  shown  in  Figure  5,  one  would  expect  the 
MBETs  to  perform  similarly  well  for  both  objects,  since  the  antenna  illuminates  no 
complex  features.  Non-specular  effects  largely  shadowed  in  Object  B,  however,  are  seen 
to  moderately  limit  the  range  of  validity  of  both  MBETs. 


Objects  A  &  B:  Tx  =  0  deg.  radiating  near-field,  W  pol 


Rx  position  (deg) 

Figure  42:  MBET  error  plot  for  a  0°  transmitter  angle.  Thin  line:  Object  A,  Thick 
line:  Object  B.  Dashed  lines  are  for  the  FI  MBET  and  solid  lines  are  for  the  F2 

MBET. 
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For  incidence  angles  that  illuminate  the  complex  features  of  Object  B,  such  as  45° 
and  135°,  shown  in  Figure  43  and  Figure  44  respectively,  both  MBETs  have  an 
expectedly  narrower  region  for  which  they  are  valid. 


Table  11:  Angular  Size  of  Valid  Regions  for  Each  of  Falconer’s  MBETs 


MBET  FI:  far-field 

MBET  F2:  radiating  near- 
field 

Tx: 

0° 

45° 

135° 

0° 

45° 

135° 

Object  A 

106° 

180° 

25° 

19° 

65° 

71° 

Object  B 

86° 

19° 

14° 

8° 

19° 

9° 

Objects  A  4  B:  Tx  =  45  deg,  relating  oear-fietd,  W  pot 


Rx  position  (deg) 


Figure  43:  MBET  error  plot  for  a  45°  transmitter  angle. 


Objects  A  4  B:  Tx  =  135  deg,  radiating  near-tieM,  W  pol 


Figure  44:  MBET  error  plot  for  a  135°  transmitter  angle. 


The  worst  case  for  Falconer’s  near- field  MBET  is  the  case  in  which  the  complex 
features  of  Object  B  are  illuminated  fully  which  corresponds  well  to  the  illumination 
angle  of  135°  (Figure  44).  From  the  figure  the  MBET  will  accurately  predict  the  RCS  in 
a  region  of  less  than  10°  of  antenna  separation. 

Following  we  present  data  that  displays  the  RCS  and  phase  error  between  the  true 
bistatic  MOM  simulation  and  the  FI  MBET  approximation  applied  to  the  monostatic 
data.  The  same  two  objects  are  used,  with  identical  antenna  configurations,  with  the 
exception  that  all  simulations  are  performed  in  the  far- field. 

Beyond  the  direct  comparison  with  the  data  approximated  by  this  MBET  against  the 
F2  MBET,  it  is  interesting  to  note  the  incredibly  good  agreement  between  the  truth  and 
the  approximation  for  Object  A  at  0°  and  45°  incidence  as  shown  in  Figure  45  and  Figure 
46.  The  phase  error,  also,  tends  to  be  very  low  in  all  cases  using  the  FI  MBET. 


Objects  A  &  B:  Tx  =  0  deg,  R  =  infinity,  W  pot 


Figure  45:  Falconer’s  FI  MBET  error  plot. 
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Objects  A  4  B:  Tx  =  45  deg.  R  =  infinity.  W  pol 


Rx  position  (deg) 

Figure  46:  Falconer’s  FI  MBET  error  plot. 

From  these  data,  we  can  approximate  the  angular  size  of  the  valid  regions  for  each 

of  Falconer’s  MBETs  based  on  the  1  dB  and  22.5°  criteria. 


Objects  A  4  B:  Tx  =  135  deg.  R  =  infinity.  W  pol 


Rx  position  (deg) 

Figure  47:  Falconer’s  FI  MBET  error  plot. 

6.5.4  Near-field  MBET  evaluation  (swept  range  bistatic  RCS) 

Finally,  to  evaluate  the  MBET  error  trends  as  a  function  of  range,  we  plot  the  RCS 
error  as  a  function  of  range  for  a  10°  and  45°  fixed  bistatic  angle.  The  ranges  correspond 
to  the  radiating  near  field  defined  in  equation  (80)  for  each  object. 
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fa *♦>“  . 


Figure  48:  Object  A  MBET  vs.  MoM 


Object  A:  Freq:  8  GHz,  Pol:  W.  Tx  Odeg 


Figure  49:  Object  A  range  dependent  error 


Object  B:  Freq:  8  GHz,  Pol:  W.  Tx  0  deg 


Figure  50:  Object  B  MBET  vs.  MoM 


•  4  S.*=\  I  T  V  u  •» 


Object  B:  Freq:  d  GHz,  Pol:  W.  Tx:  0  deg 


Figure  51:  Object  B  range  dependent  error  plot 

6. 5. 5  Conclusions 

These  simulations  have  produced  a  good  idea  of  the  conditions  for  which 
Falconer’s  extended  near-field  and  CW  MBET  are  applicable.  By  varying  the  bistatic 
angle,  range,  and  illuminated  object’s  complexity,  we  can  better  understand  at  what  point 
the  MBETs,  and  consequently,  the  Physical  Optics  and  discrete  scattering  center  models 
fail  to  predict  the  monostatic  cross-section  given  a  bistatic  reference  measurement.  Here 
it  seems  the  FI  MBET  performs  better  for  the  conditions  presented,  although  an 
investigation  into  the  performance  of  the  F2  MBET  under  far-field  conditions  would  be 
in  order.  The  upper  bound  on  magnitude  error  (1  dB)  seems  to  be  more  restrictive  than 
the  phase  bound,  which  would  tend  to  mitigate  the  importance  of  the  phase  correction 
presented  in  F2. 

It  has  been  clearly  shown  that  the  complexity  of  the  scatterer  will  severely  restrict 
the  region  of  validity  for  each  MBET,  where  shadowing  and  multi-path  interactions 
render  the  assumptions  behind  the  theorems  invalid.  The  wide  disparity  between  the 
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agreement  for  Object  A  and  Object  B  (Figure  46,  for  illumination  angle  of  45  )  illustrates 
this  point. 

It  is  not  known  what  effect  the  coarse  quantization  of  the  equivalent  monostatic 
range  and  frequency  values  has  on  the  F2  approximation.  Further  work  will  permit 
simulations  utilizing  more  quantization  levels  in  order  to  more  closely  represent  the 
MBET. 

Further  MoM  simulations  are  necessary  in  order  to  draw  a  clear  conclusion 
regarding  the  maximum  bistatic  angle  and  object  complexity  for  which  the  MBETs  begin 
to  break  down.  It  has  been  clearly  demonstrated  that  a  complex  object  (Object  B)  is  not 
well  predicted  by  any  MBET,  whether  the  approximation  is  based  on  a  scattering  center 
model  or  Physical  Optics  model  of  the  object.  Where  shadowing  and  multi-path  effects 
are  dominant,  the  errors  involved  in  the  prediction  of  a  bistatic  measurement  by  a 
monostatic  measurement  processed  through  a  near-field  or  far-field  MBET  are  great. 
Conversely,  a  simple  object  (Object  A)  which  is  well  predicted  by  the  Physical  Optics 
and  scattering  centers  model  is  likewise  well-predicted  by  the  MBET,  even  at  bistatic 
angles  greater  than  90°. 

Further  work  toward  the  end  of  characterizing  these  MBETs  for  their  applicability 
in  varying  bistatic  configurations  would  include  extensive  monostatic  MoM  simulations, 
with  finely  tuned  range  and  frequency  parameters,  so  that  the  effect  of  range  and 
frequency  adjustments  in  Falconer’s  MBET  can  be  well  understood.  In  addition,  several 
objects  with  varying  degrees  of  complexity  based  on  the  amount  of  shadowing  and  multi- 
path  effects  should  be  evaluated  simultaneously. 
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APPENDIX  B:  MATLAB  scripts 

The  following  is  a  complete  list  of  the  MATLAB  routines  that  were  used  in  the 
processing  of  measured  and  theoretical  data,  as  well  as  in  performing  the  calibrations. 
Where  outside  expertise  was  used  in  creating  the  function,  appropriate  credit  is  given. 

I CAL  D  Compute  4x4  calibration  matrix  (C)  from  the  subtracted 
calibration  object 
%  files 

%  C  =  cal_d (disk, d i s k_d i a m_cm , me  s h ver t , me  s h t i 1 1 ) 

% INPUTS: 

%disk — 2x2xn__freq  scattering  matrix  of  disk 
%disk_diam_cm--diameter  of  disk  in  cm 

%me  s  h  v  e  r  t  -  -  2  x  2  x  n__  f  r  e  q  scattering  matrix  of  vertical  mesh 
%meshtilt — 2x2xn_freq  scattering  matrix  of  tilted  mesh 
%f_sfcart_GHz- -start  frequency  in  GHz 
% f  stop  GHz — stoo  frequency  in  GHz 
% OUT PUT : 

%C---4x4xnfreq  calibration  coefficient  matrix 

function  C  =  cal_d (M_disk, disk_diam_cm, Mjneshvert , M_meshtilt ) 

global  freq_start 
global  freq_stop 
global  n_freq 
global  n_ang 
global  C 

if  n_ang>l 

warning (  1  Calibration  files  contain  more  than  one  angle — this  routine 
will  utilize  only  the  first1) 

M_disk  =  squeeze (M_disk (:,:,:,  1)  )  ; 

M_meshvert  =  squeeze  (M__meshvert  1)  )  ; 

M_meshtilt  =  squeeze (M_meshtilt 1)  )  ; 
else 

M_disk  =  squeeze (M_disk) ; 

M_meshvert  =  squeeze (M_meshvert)  ; 

M_meshtilt  =  squeeze (M_meshtilt) ; 

end 

C  =  zeros (4,4, n_freq) ;  '^initialization  of  4x4  diagonal  distortion 
matrix 

%get  PO  prediction  for  calibration  disk 

S_disk  =  calcdisk (disk_diam_cm, freq__start, f req_stop, n_f req) ; 

%ca lculate  the  co-poiar  calibration  terms 

0(1,1,:)  =  M_disk(l,l,  :)  ./S_disk(l,l,  :); 

C  (4 , 4 ,  : )  =  M_disk (2,2,  : )  . /S_disk (2,2,  : ) / 

‘^calculate  the  absolute  magnitude  of  the  vertical  mesh 

S_meshvert  =  zeros (2 , 2 , n_f req) ; 
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S_meshvert (2,2,:)  =  M_meshvert (2,2,:)  . \  )  ;  * 

% calculate  the  absolute  magnitude  of  the  tilted  mesh 

S_meshtilt  =  zeros (2 , 2 , n_f req) ; 

%with  a  qua  si -mono  static  angle  cf  5  degrees,  the  cross-polar  terms  are: 

theta  =  5*pi/180; 

S_meshtilt (1,2, : )  =  - l/2*S_meshvert (2,2, :) /cos { theta/2 ) ;  ^differs  by 
factor  of  - 1  f  r cm  b i s  s imp o 1 . p r o 
S_meshtilt (2, 1,  : )  =  S_meshtilt ( 1 , 2 , : )  / 


% solve  for  the  remaining  diagonal  terms  of  C 

C  ( 2 , 2 , : )  =  M_meshtilt (1,2, : )  .  /S_meshtilt (1,2,  : )  ; 
C  (3 , 3 , : )  =  M_meshtilt (2,1,  : )  . / S__meshtilt (2,1,  :  )  ; 


%  assuming  the  second-order  crosspc-lar  terms  are  negligible 
%  (the  s i mp .1  e - p o  1  a r i.. metric  assumption)  ,  we  can  leave  the  elements 


%C  (1,2) 
%C  (2,2, 
%  C  ( 3 ,  3 , 


0(1,3),  C  ( 1,  4  )  ,  0(4,1),  C  (4, 2)  ,  and  C(4,3)  as  zero. 
~  Mrnesh  (2,1,  :  )  .  /Smesh  (2,1,  : )  ; 

=■  Mraesh  (1,2,  :  )  .  /Smash  (1,2,:); 


% similarly,  leave  the  crosspolar  terms  C(2,l),  C(2,3),  C(2,4), 
%C (3, 1 ) ,  C (3, 2) ,  and  C(3,4)  as  zero. 


%Christopher  J  Bradley 
%2d  Lt,  USAF 
%AFIT/ENG 

%cbradley@ieee . org 
%  2  9  Septe mb e  r  20  0 0 
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%CAL  D__ MOD  Compute  4x4  calibration  matrix  (C)  from  the  subtracted 
*  ca  1.  i b r a  t  i  on  c b j  e  c t  f  i  1  e s 

cThis  function  is  a  modification  of  the  EMSL  simple  polarimetric 
% calibration  procedure  - 

%C  =  cal_d(N_disk,  P_d.isk,  N_dihedral r  F_dihedral } 

^INPUTS:” 

%disk'-~2x2xn  freq  scattering  matrix  of  disk 

tdihedral — 2x2xn_freq  scattering  matrix  of  dihedral  tilted  22.5  deg 
% OUT PUT : 

%C — 4x4xnf req  calibration  coefficient  matrix 

function  C  =  cal_d__mod (N_disk, P__disk, N_dihedral , P_dihedral ) 

global  freq_start 
global  freq_stop 
global  n_freq 
global  n_ang 
global  C 

if  n_ang>l 

warning (  1  Calibration  files  contain  more  than  one  angle--this  routine 
will  utilize  only  the  first1) 

N_disk  =  squeeze (N_disk( 1)  )  ; 

N_dihedral  =  squeeze (N_dihedral (:,:,:, 1) ) ; 
else 

N_disk  =  squeeze (N_disk) ; 

N_dihedral  =  squeeze (N_dihedral ) ; 

end 

C  =  zeros (4,4, n_f req) ;  %  initialization  of  4x4  diagonal  distortion 

matrix 

%get  PO  prediction  for  calibration  disk 

S__disk  =  P_disk; 

%calculate  the  co-pclar  calibration  terms 

0(1,1,:)  =  N_disk (1,1, : ) . / S_disk (1,1, : ) ; 

C  (4 , 4 ,  : )  =  N_disk (2,2,  : )  . /S_disk (2 , 2 ,  :  )  ; 

^calculate  the  absolute  magnitude  of  the  tilted  dihedral 

S_dihedral  =  zeros (2 , 2 , n_f req) ; 

%with  a  quasi-monostatic  angle  of  5  degrees,  the  cross-polar  terms  are: 

S_dihedral (1,2,:)  =  P_dihedral ( : , 6 ) ;  % vh  pel 

S_ dihedral (2,1,:)  =  P_dihedral ( : , 5) ;  %hv  pel 

Isolve  for  the  remaining  diagonal  terms  of  C: 

C  (2 , 2 , : )  =  N_dihedral (1,2,  : )  .  / S_dihedral (1,2, : ) ; 

C  ( 3 , 3 , : )  =  N_dihedral (2,1, : )  . /S_dihedral (2,1, : ) ; 

^Christopher  J  Bradley 
%2d  Lt,  USAF 
%  A  F I T  /  E  N  G 

%  cb  r  a  d  1  e  y  @  i  e  e  e .  o  r  g 
%7  February  2001 
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%‘CAL  D  Compute  4x4  calibration  matrix  (C)  from 
c  a 1 i b  r a t i o n  ob j e  c  t 

'ft tries  {Type-i  calibration  version) 

%  C  -•  c  a  1  d  ( d  i  s  k ,  di  s  k  _  d.  i  am  cm ! 

c.  T  Isjpf'f'pQ  » 

%disk--2x2xn_freq  scattering  matrix  of  disk 
%  d i s  k  d i am  cm — d i ame  t  e  r  o i  disk  i n  cm 
% f __st ar t_GH z- -s t art  frequency  in  GH z 
%f  stop  GHz- -stop  frequency  in  GHz 
% OUTPUT? 

% C - - 4 x 4 xn f r e q  calibration  coefficient  matrix 


th 


C"  O  Liiw’ 


abtracted 


function  C  =  cal_e  (M__disk,  disk_diam_cm) 


global  freqLstart 
global  freq_stop 
global  n_freq 
global  n_ang 
global  C 

if  n_ang>l 

warning ( 'Calibration  files  contain  more  than  one  angle — this  routine 
will  utilize  only  the  first1) 

M_disk  =  squeeze (M_disk ( 1)  )  ; 
else 

M_disk  =  squeeze (M_disk) ; 

end 

C  =  zeros  (4,4,  n_f req)  ;  -^initialization  of  4x4  diagonal  distortion 

matrix 


%get  PO  prediction  for  calibration  disk 

S__disk  =  calcdisk (disk_diam_cm,  freq^start ,  f  reqjstop,  n__freq)  ; 

%ca.lcula.te  the  co-pc-lar  calibration  terms 

C  ( 1 , 1 ,  : )  =  M_disk (1,1,  : )  • / S_disk (1,1,:); 

C  ( 4 , 4  ,  :  )  =  M_disk (2 ,  2 ,  : )  . / S_disk (2,2,  : )  ; 

%Christopher  J  Bradley 
%2d  Lt,  USAF 
%AFTT/ENG 

%cbradley@ieee . org 
December  2000 
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%CAL_EXEC_C 

% Per form  the  matrix  multiplication  to  calibrate  a  measurement 
%using  the  C- matrix  (4x4  system  distortion  matrix) ,  given 
%the  distortion  matrix  and  raw  measurement. 

%S  =  ca.l_exee_c  (M) 

% INPUT: 

%M  uncalibrated  data,  with  subtraction  and  gating  already  performed 
*Ot»T  ir’UT : 

%S  calibrated  data 


function  S  =  cal_exec_c (M) 
global  n_freq 
global  n_ang 
global  C 

S  =  zeros (size (M) ) ; 

% i n v e  r t  C -ma  t r i x 

for  m=l:n_freq 

C_inv(:,:,m)  =  inv (squeeze  (C  (:,:, m) ))  ; 

end 

q=waitbar (0,  ’Calibrating. .  .  ’ )  ; 
for  m=l:n_freq 
for  n=l.-n_ang 

M_vec  =  squeeze ( [M ( 1 , 1 , m, n) ;M  ( 1 , 2 , m, n) ; M (2 , 1 ,m, n) ;M (2 , 2 , m, n) ] ) ; 
S_vec(:,m,n)  =  C_inv ( : , : ,m) *M_vec; 

end 

waitbar (m/n_f req) ; 

end 

close (q) 

disp( ’Converting  to  scattering  matrix  format. . . 1 ) 

S  =  convert_scat (S_vec) ; 

% Christopher  J  Bradley 
%2d  Lt,  USAF 
%AFIT/ENG 

%  cbradley @ ieee . o  rg 
%  2  9  Sept  ember  2  0  0  0 
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%2 d  Lt  Christopher  J  Bradley 
%AFIT/ENG 

^Calibration  Method  A:  Whitt,  Ulaby,  Pol at in,  &  Liepa,  "A  General 
Polar imetric 

%Radar  Calibration  Technique,  If  IEEE  Transactions  on  Antennas  and 
Propagation, 

%vol.  39,  no.  1,  Jan  1991, 

%This  is  a  general  f ully-polarimetric  monostatic/bistatic  calibration 
technique 

%that  makes  no  assumptions  regarding  the  form  of  the  scattering 
matrices  of  the 

% reference  object  (excepting  that  one  of  the  matrices  is  invertible)  or 
of  the 

%form  of  the  system  distortion  matrices  (R  and  T)  .  It  requires  3 
objects  with 

^precisely  known  scattering  matrices. 


function  [phi , c , R , T] =cala ( PI , P2 , P3 , N1 , N2 , N3 ) 


lObjects  (3) : 

S; 

vertical} 

vertical) 


m e  t  a 1 1 i c  s  phe r e 

long  metallic  cylinder  (45  deg.,  horizontal,  or 
squat  metallic  cylinder  (45  deg.,  horizontal,  or 


% Input s : 

( c  o  1  urrm  v  e  c  t  c  r  s } 

%nT: 

Measured  scattering  matrix  of 

calibration  object  1 

(invertible) 

%P1 : 

T h eoretical  s c a 1 1 e r i n g  ma t r i x 

of  calibration  object 

“5 

(invertible) 

?;N2: 

Measured  scattering  matrix  of 

calibration  object  2 

%P2 : 

Theoretics 1  s c a 1 1  e r i n g  ma t r i x 

of  calibration  object 

z 

%N3: 

Measured  scattering  matrix  of 

calibration  object  3 

%P3: 

Theore t ica 1  s cat t e ring  ma t r ix 

of  calibration  object 

0 

%nrfreq: 

number  of  frequency  points 

^Outputs : 

iphi :  Absolute  phase  correction 

% e :  Magnitude  cor r ect: i on 

%R:  Receiver  distortion  matrix 

Transmitter  distortion  matrix 

^construct  the  matrices  necessary  to  perform  the  calibration: 

Nt=inv  (Nl)  *N2; 

NtB=inv (Nl) *N3; 

Pt=inv(Pl) *P2 ; 

PtB=inv(Pl) *P3; 

[Xt , LtP] =eig (Pt )  ; 

[Yt , Lt] =eig (Nt) ; 

[XtB, LtPB] =eig (PtB) ; 

[YtB, LtB] =eig (NtB) ; 

Nr  =  N2*inv(Nl); 

NrB  a  N3  *inv (Nl )  ; 

Pr  =  P2*inv(Pl); 

PrB  =  P3*inv (PI) ; 

[Xr , LrP]  =  eig(Nr); 
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[Yr,Lr]  =  eig(Pr); 

[XrB, LrPB]  =  eig(NrB); 

[YrB, LrB]  =  eig(PrB); 

%equations  for  relative  phase  of  the  reference  objects  (not  used  here 

dif f_philphi2  =  angle (LtP) -angle (Lt) ; 
dif f_philphi3  =  angle (LtPB) -angle (LtB) ; 
di f  f _phi 2phi 1  =  angle (LrP) -angle (Lr) ; 
dif f_phi3phil  =  angle (LrPB) -angle (LrB) ; 

^identify  each  element  of  Xt  and  Yt  and  assign  them  to  scalar 
variables : 

xll=Xt(l,l) ;xl2=Xt (1,2) ; x2 l=Xt (2,1) ;x22=Xt (2#2) ; 
xllB=XtB (1,1) ;Xl2B=XtB (1,2) ;x21B=XtB (2 , 1 ) ;x22B=XtB ( 2 , 2 ) ; 
yll=Yt (1,1) ;yl2=Yt (1,2) ;y21=Yt (2,1) ; y22=Yt (2 , 2 } ; 
yllB=YtB (1,1) ;yl2B=YtB (1,2) ;y21B=YtB (2,1) ;y22B=YtB (2,2) ; 

^formulae  for  ratios  of  the  elements  of  the  matrix  of  constants,  "C": 
rat_c2cl ( 1)  =  (xll*x21B-x21*xllB) * (y22*yllB~yl2*y21B) /  ... 

(x22  *xllB-xl2  *x21B) / (yll*y21B-y21*yllB) ; 
rat_c2cl (2 )  =  (xll*x22B-x21*xl2B) * (y22*yl2B-yl2*y22B) /  ... 

(x22 *xl2B-xl2 *x22B) / (yll*y22B-y21*yl2B)  ; 
rat_c2BclB(l)  =  (x22*xllB-xl2*x21B) * (y22*yl2B-yl2*y22B) /  ... 

(x22*xl2B-xl2*x22B) / (y22 *yllB-yl2*y21B) ; 
rat_c2BclB (2 )  =  (xll*x21B-x21*xllB) * (yll*y22B-y21*yl2B) /  ... 
(xll*x22B-x21*xl2B) / (yll*y2 lB-y21*yllB) ; 

tone  of  these  ratios  may  be  zero--if  so,  use  the  other  one: 

if  rat_c2cl (1) ==0 

use_c_rat  =  rat_c2cl (2) ; 
else 

use_c__rat  =  rat_c2cl(l); 

end 

if  rat_c2BclB(l)=-0 

use__cB_rat  =  rat_c2BclB  (2 ) 
else 

use_cB_rat  =  rat_c2BclB  (1) ; 

end 

if  ~isf inite (use_c_rat)  |  ~isf inite (use_cB_rat } 

warning ( *  "c"  ratios  are  infinite — system  is  indeterminate*) 

end 

if  rat_c2cl==0  |  rat_c2BclB~=0 

warning (* "c"  ratios  are  zero — system  is  indeterminate*) 

end 


%i f  Pt  and  PtB  have  more  than  one  eigenvalue  in  common, 

1 the  system  is  indeterminate: 

if  [Xt (1,1)  Xt ( 1 , 2 ) ]  ==  (  [XtB(l,l)  XtB (1,2)]  |  [XtB (2,1)  XtB(2,2)] 

)  ... 

&  .  .  . 

[Xt (1,1)  Xt (1,2)]  ==  (  [XtB(l,l)  XtB (1,2)]  I  [XtB (2,1)  XtB(2,2)]  ) 
warning(*Pt  and  PtB  have  common  eigenvalues — system  is 
indeterminate  * ) 
end 

cl  =  det (Yt) / (xll*y22-use_c_rat*xl2*y21) ; 
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c2  =  det (Yt) / (xll*y22/use_c_rat-xl2*y21) ; 

^diagonal  matrix  of  constants  to  solve  for  Tx  distortion  matrix: 

C  =  [cl  0;  0  C2] ; 

check  =  (cl*xll*y22-c2*xl2*y21) /det (Yt) ; 
if  round (100*check) /100  ~=  1 

warning ( 1  first  element  of  matrix  T  should  be  ”1",  but  it  isntT) 

end 

^transmitter  distortion  matrix: 

T  =  Xt*C*inv (Yt ) ; 


%The  development  to  solve  for  the  receiver  distortion  matrix,  "R",  i 
^almost  identical  to  that  of  the  transmitter  distortion  matrix,  "T” , 
%except  the  assignment  of  Nr,  Pr,  Nr 5,  and  PrB  (above)  is  slightly 
%dif f erent , 

%identify  each  element  of  Xr  and  Yr  and  assign  them  tc  scalar 
variables : 

xll=Xr (1,1) ;xl2=Xr (1,2) ;x21=Xr(2,l) ;x22=Xr(2,2) ; 
xllB=XrB (1,1) ;xl2B=XrB (1,2) ;x21B-XrB (2 , 1) ;x22B=XrB (2 , 2 ) ; 
yll=Yr  (1,1)  ;yl2-Yr  (1,2)  ;y21=Yr(2,l)  ,*y22=Yr(2,2)  / 
yllB=YrB  (1,1)  ;yl2B=YrB  (1,2)  ; y2 lB=YrB  (2,1)  ;y22B=YrB  (2,2)  ; 

%  formulae  for  ratios  of  the  elements  of  the  matrix  of  constants,  "D" 

rat_d2dl (1)  =  (xll*x21B~x21*xllB) * (y22*yllB-yl2*y21B) /  ... 

(x22*xllB-xl2*x21B) / (yll*y2 lB-y21*yllB) ; 
rat_d2dl (2)  =  (xll*x22B-x21*xl2B) * (y22*yl2B~yl2*y22B) /  ... 

(x22*xl2B-xl2*x22B) / (yll*y22B-y21*yl2B) ; 
rat_d2BdlB(l)  =  (x22*xllB-xl2*x21B) * (y22*yl2B-yl2*y22B) /  ... 

(x22  *xl2B-xl2  *x22B) / (y22*yllB-yl2*y21B)  ; 
rat_d2BdlB (2)  =  (xll*x21B-x21*xllB) * (yll*y22B-y21*yl2B) /  ... 
(xll*x22B-x21*xl2B) / (yll*y21B-y21*yllB) ; 

%one  of  these  ratios  may  be  zero — if  so,  use  the  other  one: 

if  -isfinite (rat_d2dl (1) )  |  rat_d2dl (1) ==0 

use_d_rat  =  rat_d2dl(2); 
else 

use_d_rat  =  rat_d2dl(l); 

end 

if  -isfinite (rat_d2BdlB ( 1) )  |  rat_d2BdlB ( 1 ) ==0 

use__dB_rat  =  rat_d2BdlB  (2 )  ; 
else 

use_dB_rat  =  rat_d2BdlB (1) ; 

end 

if  ~isf  inite  (use_d_rat)  |  ~isf  inite  (use_dB__rat ) 

warning  (  '  "d"  ratios  are  inf  inite- -system  is  indeterminate1) 

end 

if  use_d_rat==0  |  use_dB__rat==0 

warning (’ "d"  ratios  are  zero — system  is  indeterminate') 

end 

dl  =  det (Yr) / (xll*y22 -use_d_rat*xl2*y21) ; 
d2  =  det (Yr) / (xll*y22/use_d_rat-x!2*y21 ) ; 

% diagonal  matrix  of  constants  to  solve  for  Rx  distortion  matrix: 
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•  ;• 


D  =  [dl  0;  0  d.2]  ; 

check  =  (dl*xll*y22-d2*xl2*y21) /det (Yr) ; 
if  round (100*check) /100  ~=  1 

warning ( 1  first  element  of  matrix  R  should  be  "l",  but  it  is  not1) 

end 

%  r  e  ce  i  v  e  r  d  i  s  t  ort  i  on  in  a  t  r  i  x : 

R  =  Xr*D*inv (Yr) / 

% optimize  precision  for  the  constant  c: 

%best_P  -  input (* Which  theoretical  matrix  to  use  for  "c"  (most  precise 
?  *)  ; 

%ba.ck-substit.ute  to  find  magnitude /phase  correction  factor: 

%Note:  It  is  best  to  use  a  diagonal  element  (1  or  4)  to  find  1 c 1 

ibecause  the  diagonal  elements  are  much  more  stable  with  noise 
%tnan  the  off-diagonal  elements 

%eval ( [ ’ c_complex  =  N { 1 ) ’  num2str (best_  P)  1 . / (R ( 1 ) *P ( 1 ) ’ . . . 

%  num2str (best_P)  ’ *T ( 1 } ) ;  1  ] ) 

c_complex  =  N1 (4) . / (R (4) *P1 (4) *T (4) ) ; 

^magnitude  correction  factor: 

c  =  abs (c_complex) ; 

% pha s e  co r r e c t ion  fa c t or  {re 1 a t i v e } : 

phi  =  angle (c_compl ex) ; 

%Ch ristopher  J  Bradley 
%2d  Ltf  USAF 
%AFIT/ENG 

%cbradley@ieee . org 
%29  September  2000 
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ICALE  Implement at ion  of  Type-1  calibration  technique 
%2d  Lt  Christopher  J  Bradley 
%AFIT/E m 

%Calibration  Method  E:  Jersak.,  "Bistatic,  Fully  Polarimetric  Radar 
C  r  o  ss-Se c t ion 

^Calibration  Techniques  and  Measurement  Error  Analysis,  IT  Ph.D. 
d i s  s  e  r t  a  t i on , 

^University  of  Texas  at  Arlington,  Dec.  1993,  pp.  94-96. 

%This  is  a  very  basic  polarimetric  calibration  technique  requiring  only 
a  single 

%re£erence  object  (sphere) .  It  does  not  account  for  cross -pel  leakage 
of  the 

%antennas,  and  at  bistatic  angles  at  which  there  is  a  null  in  the 
return  from 

% the  sphere,  the  technique  cannot  be  performed.  It  has  the  advantage 
that  it 

%is  completely  insensitive  of  object  orientation. 


% INPUTS : 

%M0:  measurement  of  sphere  (matrix) 

%3G:  predicted  value  of  sphere  (matrix) 


%OUTPUTS : 

%Ch:  modified  distortion  matrix:  pointwise  (hadamard)  product  with 

%  unknown  object  scattering  matrix  reveals  the  calibrated  value 

%C)bjects  (1)  :  metallic  sphere 

function  Ch  =  cale(MO,SO) 
global  n_freq 


Ch 


zeros(2,2);  %init ialization  of  4-element  distortion  matrix 


%  f  i  1 1  distort  i  o  n  ma  t  r  i.  x : 

Ch ( 1 , 1 )  =  SO (1,1) /MO (1,1) ; 
Ch (2,2)  =  SO (2,2) /MO (2,2) 
Ch (1,2)  =  SO (1,2) /MO (1,2) 
Ch ( 2 , 1 )  =  SO (2,1) /MO (2, 1) 


^Christopher  J  Bradley 
%2d.  Lt,  USA.F 
%AFIT/ENG 

%  cbradl ey@ i eee . org 
%29  September  2000 
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%GET___C  Compute  system  distortion  matrix  C  using  the  EMSL  simple 
%polarimetric  calibration  technique.  No  inputs  are  given,  and  the 
k output  (C)  is  returned  to  a  global  variable.  Utilizes  ,,cal_d.mff 
%to  compute  the  C  matrix. 

function  get_C 
global  C 
global  n_freq 
global  freq_start 
global  freq_stop 

f  =  linspace (freq_start , f req_stop, n_f req) ; 

disp(’ Input  raw  calibration  files — ’) 

Cl_s  =  input  (' calibration  disk  file:  h’s’); 

C2_s  =  input ( ’ empty  room  file :  '  ,  ’  s  ' )  ; 

C3_s  =  input ( ’vertical  mesh:  1 , ’s’) ; 

C4_s  =  input ( ’tilted  mesh:  ’,'s’); 

diam  =  input (’ calibration  disk  diameter  (cm):  ’); 

Cl  =  read_emsl (Cl_s) ; 

C2  =  read_emsl (C2_s) ; 

C3  =  read_emsl (C3_s) ; 

C4  =  read_emsl (C4_s) ; 

Cl  =  C1-C2 ; 

C3  =  C3-C2; 

C4  =  C4-C2; 

gl  =  input(’Gate  start  time  (ns):  ’); 

g2  =  input (’Gate  stop  time  (ns):  '); 

t  =  linspace (gl , g2 , n_f req) ; 

Cl_ga  =  time_gate (Cl , f , t ) ; 

C3_ga  =  time_gate (C3 , f , t) ; 

C4_ga  =  time_gate (C4 , f , t ) ; 

C  =  cal_d(Cl_ga,diam,C3_ga/C4__ga)  ; 

s  =  input(’Save  coefficient  matrix  C  as  a  mat  file?  * ,’s*); 
if  s ( 1 ) == ’ y * 

dir  =  input (’ directory :  ’,’s'); 
name  =  input  (  ’  File  name :  ’  ,  ’  s  ’ )  / 

eval ( [name,  ’=C; 1 ] ) 

eval ( [  ’save  dir,  * \ ’, name,  ’  ’,name]) 
size (C) 
else 
end 

^Christopher  J  Bradley 
%2d  Lt,  USAF 
%AFIT/EHG 
%cbradley@ieee . org 
%29  September  2000 
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% CALC DISK 

%generate  theoretical  scattering  matrix  for  a  metallic  circular  disk 
% (Physical  Optics  approximation) ,  in  amplitude,  NOT  dBsm 
%  Inputs: 

%diam  -™  diameter  in  centimeters 
ifreq_ start  =  start  frequency  (GHz) 

%freq_stop  =  stop  frequency  (GHz) 

%nfreq  =  number  of  frequency  samples 

function  S  =  calcdisk (diam, f req_start , f req_stop, nf req) 

% convert  diameter  (m)  to  radius  (cm) 

r  =  diam/2/100; 

% i ni t i a i i z e  s t r uct ure 

S  «  zeros (2 , 2 , nf req) ; 

% p r opa g a  t i on  v e 1 o ci t y 

CO  =  2 . 998E8 ; 

% wave length 

lambda  =  c 0  . /linspace (f req_start*lE9 , f req_stop*lE9 , nf req) ; 

Scross-s e  c  t ion a 1  are a 

A  =  pi* (r) ^2; 

%co -polar  components  are  equal,  cross-polar  are  zero 

S  ( 1 , 1 , : )  =  2*i*sqrt (pi) *A. /lambda; 

S (2 , 2  ,  :)  =  S ( 1 , 1 ,  :)  ; 
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% CALIBRATE  DATA  Calibrate  gated  and  subtracted  data 
%The  routine  is  interactive  within  the  Mat lab  workspace/  and 
%for  the  raw  data  files,  which  must  be  in  the  unprocessed  EMS 
format 

% S  “  calibrate  data (M) 

% INPUT : 

%M  --  input  scattering  matrix:  2  x  2  x  n_freq  x  n__ang 
* OUT PUT : 

%S  —  calibrated  scattering  matrix  of  the  same  dimensions 

function  S  =  calibrate_data (M) 

global  freq_start 
global  freq_stop 
global  n_freq 
global  n_ang 
global  ang 
global  C 

file  =  input (  1  Enter  file  name  of  data:  ','s'); 
read_hdr (f ile) ; 

f  =  1  inspace  (freq_s tart ,  freq_s top,  n__freq)  ; 

% calibrate : 

docal  =  input (  1  Perform  new  calibration?  *  #  *  s 1 )  ; 
if  docal (1) == * y ’ 

which_cal  =  input (  1  Calibrate  with  sphere?  ','s'); 
if  which_cal(l)  ==  ' n ’ 

get_C; 
else 

get_C_sph; 

end 

elseif  docal (1) ==' n 1 

disp('If  not  in  memory,  enter  dir  and  filename — T); 
dir  =  input (  1  directory :  ','s'); 

name  =  input ('file  name:  ','s'); 

if  is empty (name) 

%do  nothing — C  is  a  global  variable, 

else 

eval ( [  1  load  1 , dir,  T \ 1 , name] ) 
eval  {  [  '  C~  ’  ,  name ,  1 ;  '  ]  ) 
end 

end 

read_hdr (file) ; 

S  =  cal_exec_C (M) / 
read_hdr (file) ; 

in_2  =  input ( 'Plot  calibrated  data?  ','s’); 
if  in_2(l)  ==  'y' 

pol  =  input (’View  polarization:  ','s') 
sel_f  =  input ( 'View  frequency  (GHz):  '); 

[sel_f , sel_f_ind]  =  min (abs (f -sel_f ) ) ; 
if  n_ang>l 

sel_a  =  input  (  'View-  angle  (deg):  '); 

[sel_a, sel_a_ind]  =  min  (abs  (ang-sel_a)  )  ; 
view_data_a  =  get__angle_sweep  (S,pol ,  sel_f_ind)  ; 
view_data_f  =  get_f req_sweep (S , pol , sel_a_ind) ,* 
figure (1) 
subplot (2,1,1) 
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G,  pi 


plot (ang, view_data_a) 

title([’RCS  vs.  angle,  '  ,  num2str  (f  (sel_f_ind)  )  ,  '  GHz1]) 
subplot (2,1,2) 
plot { f , view_data_f ) 

title ([’RCS  vs.  frequency,  '  ,  num2str (sel_a) ,  1  deg1]) 
elseif  n_ang==l 

view_data  =  get_f req_sweep (S,pol , 1) ; 
figure (1) 
plot (f , view_data) 
title(’RCS  vs.  frequency1) 
end 

end 

s  =  input ('Save  data  as  a  .mat  file?  ' ,1s'); 
if  s  (1)  ==  '  y* 

dir  =  input (' Directory :  '  ,'s'); 

name  =  input  (  '  File  name :  '  ,  ’  s  ' )  ; 
eval ( [name,  f =S; 1 ] ) 

eval ( [ ' save  1 , dir, ' \ ' , name, ’  1 , name] ) 

else 
end 

^Christopher  J  Bradley 
%2d  Lt,  USAF 
%AFIT/ENG 

% ch r a  d 1 e  y @ i eee.or g 
%29  September  2000 
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% CONVERT  FISC  Convert  . mat  files  generated  from  "fiscread.m" 
tor  T,read_f  ield_2  .m"  into  another  .mat  file  corresponding  to  the 
structure 

%of  the  measurement  scattering  matrices 
% s cat _da t a  =  co n v ert_t  i s c  { a ,  f ,  n__ang ) 

function  scat_data  =  convert_f isc (a, f , n_ang) 

% initialize  vectors: 

f_data  =  round(a (:, 1) *10000) /10000; 
f  =  round(f *10000) /10000; 
n_freq  =  length (f) ; 

scat_data  *  zeros ( [2 , 2 , n_f req, n_ang] ) ; 

%  s  o  r  t  b  y  f  r  e  qu  e  n  c  i  e  s  : 

global  m 

%  fill  scattering  matrix: 
h=waitbar { 0 ,  *  filling  matrix1); 
for  m  =  l:n_freq 

f_ind  =  f ind(f_data==f (m) ) ; 
scat_data  (2, 2 ,m,  : )  =a  (f__ind,  4) 
scat_data (1 , 1  ,m, : ) =a (f_ind, 7 ) 
scat_data (1, 2 ,m, : ) =a (f_ind, 6) 
scat_data (2 , 1 ,  m, : ) =a (f_ind, 5) 
waitbar (m/n_f req) 

end 

close  (h) 


%vv  pol 
%hh  pol 
%vh  pol 
%hv  pol 
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£scat_new  =  convert_scat  (scat ,  n_f req) 

^.convert  structures  between  scattering  matrix  and  scattering 
% vector  format  {toggle} : 

IScattering  matrix  =  [  HH  VH  ] 

%  [  Hv  VV  ] 

% Scattering  vector  ~  [  HH  j 

%  [  VH  ] 

%  [  HV  ] 

%  [  VV  ] 

%  INPUTS: 

%  s  c  a  t :  i  n  p u  t  s  t  r  u  c t  u  r  e  v a  r  i  a  b  1  e 

% OUT PUT : 

%S:  converted  scatterincr  structure 


function  S  =  convert_scat ( scat ) 
global  n_freq 
global  n_ang 

%if  the  number  of  dimensions  of  the  input  structure  is  4, 

%then  we  must  convert  to  the  3-dimensional  (scattering  matrix) 

% format : 

if  n_ang>l 
if  ndims (scat ) ==4 

S=zeros (4 , n_f req, n_ang) / 
for  m=l:n_freq 
for  n=l:n_ang 

S ( : ,m,n) ^squeeze (  [scat ( 1 , 1 , m, n) ; scat (1 , 2/m/n) ; scat (2 , l,m,n) ; scat (2,2 

,m,n) ] ) ; 

end 

end 

% conversely,  if  the  number  of  dimensions  of  the  input 
% structure  is  3,  then  we  convert  to  the  4 -dimensional 
% (scatter! n g  v e  c tor)  f o  rma t : 

elseif  ndims (scat) ==3 

n_ang=length (scat (1,1,:))/ 

S= zeros (2 , 2 , n_f req, n_ang) ; 
for  m=l:n__freq 
for  n=l:n__ang 

S  (1,  l,m,n)  =scat  (l,m,n)  ; 

S (1 , 2 ,m, n) =scat (2 ,m,n)  ; 

S (2 , 1 ,m, n) =scat (3 ,  m,  n) ; 

S (2 , 2 ,m, n) =scat (4 ,m, n) ; 
end 

end 

else 

error (’input  structure  "scat"  has  the  wrong  dimensions') 

end 

end 

if  n_ang==l 

if  ndims (scat) ==3 

S=zeros  (4,n__freq)  ; 
for  m=l:n_freq 

S ( : ,m) ^squeeze ( [scat (1 , l,m) ; scat (1, 2 ,m) ; scat (2 , 1 ,m) ; scat (2 , 2 ,m) ] ) ; 
end 

elseif  ndims (scat ) ==2 
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w  o 


S=zeros (2 , 2 ,n_f req) ; 
for  m=l:n_freq 

S  ( 1 , 1 , m) =scat ( 1 , m) ; 

S  (1 , 2 ,m) =scat (2 ,  m) ; 

S (2, l,m) =scat (3,m) ; 

S (2 , 2 ,m) =scat (4 ,m) ; 
end 
else 

error (? input  structure  "scat"  has  the  wrong  dimensions 1 ) 
end 

end 

% Ch  r i s t  op her  J  B  r a dl e y 
%2 d  Lt t  USAF 
% AFIT/ENG 

bradley@ieee.org 
9  September  2000 
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V; G, AT E  DATA  R e a ci  f  v i e w  3 n cl  /or  gate  unp t o c g s  s 0 cl  o. 3 t. 3 
§  The  routine  is  interactive  from  the  Mat lab  workspace, 
iciata  files  should  be  in  the  raw  format  from  the  EMSL  HP 
output 

%  i  n  t  o  scat  t  e  r  i  n  g  m  a.  fc  r  i.  x  f  o  r rn  a  t . 

%S  =  gat e_dat  a 
% OUT PUT : 

*S_£inal  --  processed  data,  as  defined  by  user  prompts 

function  S_final  =  gate_data 
close  all 

tar_file  =  input ( 'Directory  and  name  of  object  file:  '  , 
ER_file  =  input ( 'Directory  and  name  of  empty  room  file: 
disp (' Loading  data  files...’) 

ER  -  read_emsl (ER_f ile) ; 

tar  =  read_emsl (tar_f ile) ; 

global  object 

global  mnemonic 

global  freq_start 

global  freq_stop 

global  f 

global  n_freq 

global  n_ang 

global  ang 

tarjname  = [mnemonic,  ' , object] ; 

S  =  tar  -  ER; 
df  =  f  ( 2 )  -  f  { 1 )  ; 

p  =  input (’Plot  frequency  domain?  (y  or  n) :  '  ,'s'); 

if  p(l)  ==  * y ' 

pol  =  input (' Polarization?  ’,'s'); 
if  n_ang>l 

scat _j?lot_2D(S,pol , f ,  ’ f req ’ ) 
label (  *  f req ' , 0 , pol ) 
else 

plot  (f ,get_freq_sweep(S,pol) ) 
label {  '  freq’ , 0,pol) 

end 

elseif  p ( 1 )  ==  'n’ 
end 

g  =  input (’Apply  range  gate  (y  or  n) ?  ’,'s'); 

if  g ( 1) ==  ? y ' 

pol  =  input ('View  polarization:  ’,’s'); 
if  length (pol ) ==2 

wl  «  input (  'Time  window~~start  (ns):  '); 
w2  =  input  (  ’Time  window- -stop  (ns)  :  '); 

t__range  =  1 inspace (wl , w2 , n_f req) ; 
if  wl>w2 

error ( 'improper  windowing  parameters:  wl>w2 1 ) 

elseif  w2>round( (1/df )  ) 

warning ( 'window  is  beyond  unambiguous  range') 
end 

if  n_ang==l 
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time_view  =  tf convert (S , f , t_range ,  ’time’); 
time_view  =  get_f req_sweep ( time_view, pol # 1) ; 
figure (1) 

plot (t__range, time_view) 
label (  ’ time  1 , 0,pol) 
elseif  n_ang>=l 

time_view  «  tf convert (S , f , t_range ,  ’time1); 
figure (1) 

scat_plot_2D (time_view,pol , t_range,  ’time1 ) ; 
label { ’ time ' , 0 , pol ) 

end 

end 

gl  =  input (  'Select  gate  cutoff — start:  ’); 
g2  =  input  (  1  Stop  :  1 )  ; 

if  gl>g2 

error ( 1  improper  windowing  parameters:  wl>w2 1 ) 
elseif  g2>(l/df) 

warning ( 'gate  is  beyond  unambiguous  range’) 
end 

S_f inal  =  time_gate (S, f , 1 inspace (gl , g2 , n_freq) ) ; 
v  =  input ('View  data?  ’,’s’); 
if  v(l)«'y' 
if  n_ang>l 

sel_a  =  input (  ’angle  (deg):  ’); 
sel_f  =  input (  ’frequency  (GHz):  ’); 
pol  =  input ( ’View  polarization:  ’.’s’); 

tq,  sel_f_ind]  =  min (abs (f -sel_f )  )  ; 

[q, sel_a_ind]  =  min (abs (ang-sel_a) ) ; 
view_ga_a  =  get__angle_sweep  (S_f  inal ,  pol ,  sel_f_ind)  ; 
view_ga_f  =  get_f req_sweep (S_f inal , pol , sel_a_ind) ; 
figure (1) 
subplot (2,1,1) 
plot  (ang,  view_jga_a) 
label ( ’angle’ , f (sel_f_ind) ,pol) 
subplot (2,1,2) 
plot  ( f ,  view__ga_f ) 
label ( ' f req ’ , ang ( sel_a_ind) , pol ) 
elseif  n_ang==l 

pol  =  input ( ’View  polarization:  ’ , ? s ’ ) ; 

view_ga  =  get_f  req_sweep  (S__f  inal ,  pol ,  1)  ; 

figure (1) 

plot (f , view_ga) 

label ( ’ freq1 , 0,pol) 

end 

else 

end 

elseif  g (l) == ’ n ’ 

S_final  =  S; 

end 

s  =  input (’Save  data  as  a  .mat  file?  T,’s’); 
if  s (1) == ’ y 1 

dir  =  input  (  ’  Directory:  '  ,  ’  s  ’ )  / 
name  =  input (  ’ File  name :  ' , ’ s  ’  )  ; 

eval ( [name,  ’ ~S_final; 1 ] ) 
eval(['save  ’, dir , ’ \ ’, name, ’  ' , name] ) 

else 
end 
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U  CM  rfj  O  CM 


ri st op her  J  Bradley 
Lt  ,  USAF 
I? /ENG 

r adl e y @ i e e  e . o r g 
September  2000 
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%GET  ANGLE  SWEEP  Read  the  scattering  matrix  1 S 1 
^format),  select  a  single  frequency,  and  convert 
*S__dB  =  get_angle_sweep (S,  pol,  freq_index) 

% IN PUTS: 

%S  —  input:  scattering  matrix  (complex  power) 
%pol  --  desired  polarization 
% f req_in.de x  —  desired  frequency  index 
% OUT PUT : 

%3_  dS  —  vector  of  RCS  (dBsin)  for  a  single  polar, 
frequency 

function  S_dB  =  get_angle_sweep (S, pol , freq_index) 


(in  complex  power 
to  dBsm. 


zation,  single 


if  pol== ' hh ' 
r  =  1; 
c  =  1; 

elseif  pol==,vhf 
r  =  1; 

C  =  2; 

elseif  pol == ’ hv ’ 
r  =  2; 

C  =  1; 

elseif  pol== 1 vv 1 
r  =  2; 

C  =  2  ; 
else 

error ( T Unrecognized  variable  "pol"  T) 

end 


S_dB  =  20*logl0 (abs (squeeze (S (r, c, f req_index, :)))); 

% Christopher  J  Bradley 
■r2d.  JUt ,  USAF 
a AF IT /ENG 

%cbradley@ieee . org 
%29  September  2000 
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%GET___ANGLE__SWEEP  Read  the  scattering  matrix  ?  S’  (i: 
% format),  select  a  single  angular  position,  and  con 
%S__dB  =  get__f  req_sweep  ( S ,  pol ,  ang_index } 

%  I  INPUTS: 

%S  —  input  scattering  matrix  (complex  power) 

%pol  --  desired  polarization 

%ang_index  desired  angular  position  index 

function  S_dB  =  get_f req_s weep (S, pol , ang_index) 


complex  power 
rt  to  dBsm. 


if  pol=='hh’ 
r  =  1; 
c  =  1; 

elseif  pol==fvh? 

r  -  1;  • 

C  =  2; 

elseif  pol== 1 hv T 
r  =  2; 

C  =  1; 

elseif  pol == ’ vv ' 
r  =  2; 

C  =  2  ; 
else 

error ( 1  Unrecognized  variable  "pol"  ') 

end 


if  length (size (S) ) <=3 

S_dB  =  20*logl0 (abs (squeeze (S (r, c,  :))))  ; 
else 

S_dB  =  20*logl0 (abs (squeeze (S (r, c, : , ang_index) ) ) ) ; 

end 


% Christopher  J  Bradley 

•3Zd  -Lit  / 

%AFIT/ENG 

%cbradley@ieee . erg 
%  2  9  S ep t ember  2  0 0  0 
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%  I DLCONVERT 

%  input  and  convert  a  1  *£r e*  (w/o  frequency  and  time  values} 
if ile  from  IDL  and  read  into  MAT LAB  also  convert  the  values 
%  i  n t  o  c crap  1  e x  sea 1 1  e r i  ng  ma t  r i  x  f  o  rm a  t : 

%  [  HH  VH  1 

?;scat_mat  =  [  ] 

i  ’  [  HV  VV  ] 

i  The  scattering  matrix  is  (in  general}  4-dimensional — 

% s cat  jna t ( k , 1 , m ,n): 

%  k,l:  polarization  dimensions  (traditional  2x2  scattering  matrix) 

%  m:  frequency  dimension 

%  n:  angle  dimension 

iThis  function  is  different  from  "read  emsl.m"  in  that  it  uses  T.freT 
lor  ’ . tim1  files  which  do  not  contain  a  header. 

%  INPUTS: 

?;fiie_  name:  string  containing  the  directory  and  file  name 

%numfreq:  number  of  frequencies  per  sweep 


%The  output  tensor 
%  f o rm .  T o  conv e  r t 


’ s  ’  is  i n  c omp 1 e x  p owe r  ( a b s o 1 u t e , 
to  dB  -->  20*  1  og  1 0  ( abs  { s cat_ma t }  ) 


not  dB} 


function  s  =  idl convert (file_name , numfreq) 


%dats  files  always  contain  4  polarizations: 

numpol  =  4; 

fid  =  fopen(f ile_name,  1 r?,’bf); 

[a,  count]  =  f read (fid,  ’float’); 
numangle  =  count/2/numf req/numpol ; 
stmp  =  zeros (2, 2, 2/numfreq/ numangle) ; 
stmp ( : )  =  a; 

S  =  Stmp (1,  ) + j  *Stmp (2 ; 

s  =  squeeze (s); 


%Christopher  -j  Bradley 
%2d  Lt,  USAF 
%AFIT/£NG 

%cbradley@ieee . org 
% 2  9  Sept ember  2000 
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%MCHIRPZ  Chirp-z  transform 
%g  =  mchirpz (x, f , t, sgn) 

% INPUTS: 

%x:  input  vector 

%f:  vector  of  frequency  values  (or  initial  domain) 

%t:  vector  of  time  values  (or  tr as form  domain) 

%sgn :  1  to  convert  to  time,  -1  to  convert  to  frequency 
sJoaquim  Fortuny,  27  Aug  2000 

function  g=mchirpz (x, f , t , sgn) 


X=X ( : ) . ' ; 

f=f  (:) ; 

t  =  t  (  : )  ; 
n=length (x) ; 
m= length ( t) ; 
tpi=2*pi ; 

dw=tpi* (f (2) -f (1) ) ; 
dt=t (2 ) -t(l) ; 
wO=tpi*f (1) ; 
tO-t  (1)  ; 
l=n+m-l ; 

l=pow2 (floor (log2 (1) ) +1)  ; 

xa=zeros (1,1) ; 
xb=xa ; 
xc=xa; 

xa(l:n)=x. *exp ( j  *sgn* (dw*t0* (0 :n-l) +dw*dt* (0 :n-l) . A2/2) ) ; 
xa-f f t (xa) ; 

xb (1 :m) =exp ( - j  *sgn*dw*dt* (0 :m-l) .  A2/2) ; 
xb(l-n+2:l)=exp(-j*sgn*dt*dw* (1- (l-n+l:l-l) ) . A2/2) ; 
xb=f f t (xb) ; 
xc=if f t (xa . *xb) ; 
if  sgn==l 

factor=l/length(f )  ; 
else 

f actor=l ; 

end 

g-xc (1 :m)  . *exp ( j  *sgn* ( (t0+ (0 :m-l) *dt ) *w0+dw*dt* (0 :m-l)  . A2/2) )  ; 
g=g*factor; 
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%READ_EMSL  Read  the  raw,  unprocessed  data  from  the  EMSL  HP  int 
^-scattering  matrix  understandable  to  MAT  LAB .  Uses  global  varia 
^initialized  in  f,read_hdr .m" 

% S  =  read  em s 1 { f i 1 e ) 

% INPUT: 

%file  —  directory  and  file  name  of  the  EMSL  unprocessed  data 
% OUTPUT : 

%S  --  data  in  2  x  2  x  n  freq  x  n  ang  scattering  matrix  format 

function  S  =  read_emsl {f ile) 
global  n_freq 
global  n_ang 
global  n_hdr 

%read  the  header  of  the  data  file  and  give  values  to  the 
%  g 1 oba 1  var i abl e  s  a bo v e : 

read_hdr (f ile) ; 

%  u s  e  f  u 1  cons  t  a  n t  s : 

bytes  =  9;  % number  of  bytes  per  line  of  data 

bad_l i ne  s  =  3  ;  %  n u mb e  r  o  f  1 i ne s  c f  b a d  da  t  a  pe  r  po s i t i on 

Breaci  x Du  data: 

f close ( 1  all 1 ) ; 

fid  =  fopen(file,  'r’/s’l; 

f=f read (fid,  ’double1); 

%select  binary  data: 

%The  data  begins  3  lines  after  the  end  of  the  header 

hdr_length  =  bytes*  (n__hdr)  ; 
f_data_l  =  f (hdr_length+l:end) ; 

%The  raw  data  has  3  lines  of  garbage  at  the  beginning 
% o f  each  frequency  sweep  which  have  to  be  removed: 

n  =  []  ; 

for  m-1:  (bytes*n_freq+bad_lines*bytes)  :  length (f_data_l)  ; 
n  =  [n  m:m+bad_lines*bytes- 1]  ; 

end 

f_data_l  (n)  =  []  ; 

%find  and  remove  frequency  indices: 

freq  =  []  ; 
f__data  =  f__data_l; 

f_data (1 : bytes : length ( f_data) )  =  [] ; 

%{2  channels  (real,  imaginary),  4  polarizations,  n  freq  frequencie 
% convert  the  vector  into  scattering  matrix  format: 

stmp  =  zeros (2 , 2 , 2 , n_freq,n_ang) ; 
stmp ( : )  =  f_data; 

S  =  squeeze (stmp (1 +j *stmp (2 ; 

% C  h  r i s  t  op  h e  r  J  B  r a d 1 e y 
%2d  Lt,  USAF 
% A FIT/ENG 

%  cb  r  a  d  1  e  y  @  i  e  e  e .  o  r  g 
%  2 9  S  ep t  embe r  2  0 0  0 
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%  R.E  aD  i!  x.bLiJ^2 

%field__mat  -  read_field_2  (filename) 

%This  function  takes  the  . field  file  generated  by  FISC 
%and  converts  it  to  a  MAT LAB  matrix  that  is  output  to 
%the  workspace.  The  output  matrix  is  of  the  same  form 
las  the  original  file,  i.e.  7  columns: 

% frequency  :  cbs-EL  :  obs~AZ  :  VV  :  HV  :  VH  :  Hh 

%2 d  Lt  Christopher  J  Bradley 

%AFIT/ENG 

%29  Nov  2000 

function  field_mat  -  read_fi el d_2 ( filename) 

% e I erne n t  1 :  i n c i de n t  elevati o n 
%  e 1 erne n t  2 :  i n c i de n  t  a  z i mu t  h 
%element  3:  observation  elevation 
^element  4:  observation  azimuth 
%  e  1  erne n t  5  :  f  r  equency  (GHz ) 

% element  6:  vv-real 
%  e 1 emen t  7 :  vv - i ma g i na  r y 
^element  8:  hv-real 
%  e  1  e  m  e  n  t  9 :  h  v  -  i  m  a  g  i  n  a  r  y 
% element  10:  vh-real 
^element  11 :  vh-imaginary 
%  e 1 emen t  12:  hh -real 
^element  13:  hh-imaginary 


inc_el  =  []  ; 
inc_az  =  []  ; 
obs_el  =  []  ; 
obs_az  =  []  ; 
f  =  []  ; 


w_r  =  []  ;  W_i  =  []  ; 

hv_r  =  []  ;  hv_i  =  []  ; 

vh_r  =  []  ;  vh_i  =  []  ; 

hh_r  =  []  ;  hh_i  =  []  ; 

w  =  []  ; 
hv  =  []  ; 
vh  =  []  ; 
hh  =  []  ; 

tic 

data  =  textread (filename, 

1 whitespace 1 , 1  ' , 1  headerlines  1 ,5) ; 

toe 

tdata  =  textread  ( T  short_cyl___10tll .  field1 ,  ?  %s  * ,  'whitesp 
* ,  *  headeriines ' , 5 ) ; 

tic 

N  =  length (data) ; 
q_char  =  []  ; 
x_data  =  []  ; 
for  n=l : N 

x  =  data (n) ; 

q_char  =  cat(l,x{:}); 

if  length (q_char) >1 

k_open  =  findstr  (eschar,  1  (  1 )  ; 
k_clos  =  f indstr (g_char ,  '  )T); 
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k__comm  =  f  indstr  (eschar ,  T ,  1  )  ; 
if  -isempty (k_open) 

x_data(n)  =  str2num (q_char (k_open+l :k_comm-l) ) 
+  j  *str2num  (q_char  (k_comm+l  :k__clos-l)  )  ; 
else  x_data(n)  =  str2num (q_char) ; 
end 

elseif  length  (q_char  ====!) 

x_data(n)  =  NaN; 
end 

end 

%Njdats  =  length  (x  jiata)  ; 

n=l  ; 

while  n<length(x__data) 
if  isnan(x_data (n) ) 
x_data(n)  =  []  ; 
end 
n=n+l; 

end 

N_data  =  length (x_dat a) 

for  n  =  l:9:N_data 

inc_el  =  [inc_el  x_data(n)]; 
inc_az  =  [inc_az  x_data (n+1)  ]  ; 
obs_el  =  [obs_el  x_data (n+2 )  ]  ; 
obs__az  =  [obs_az  x_data (n+3 )  ]  ; 
f  =  [f  x_data(n+4)  ]  ; 
w  =  [w  x_data (n+5)  ]  ; 
hv  =  [hv  x_data (n+6)  ]  ; 
vh  =  [vh  x_data (n+7)  ] / 
hh  =  [hh  x_data (n+8)  ]  ; 

end 


field_mat  = 
f ield_mat ( : 
f ield_mat { : 
f ield_mat ( : 
f ield_mat ( : 
f ield_mat ( : 
f ield_mat ( : 
field  mat ( : 


zeros (N_data/9, 7) ; 

,1)  =  f '  ; 

,2)  «  obs_el ' ; 

,3)  =  obs_az ' / 

,  4 )  =  w  ’  ; 

,  5 )  =  h v 1  ; 

,  6 )  =  vh '  ; 

/  7 )  =  hh' ; 


$  save  as  .mat  file 

leva!  ( [filename  (1 :  end-7  5  ,  ’___field  =  f  ield_mat;  ’  ]  ) 

%eval { [ *  save  i : \igars2000\f isc\£ iscmat \ 1 ,  . . . 

%  filename  ( 1 :  end-7 )  ,  '___field  1 ,  filename  ( 1 :  end- 7 } ,  1 

?;eval  {  [ 1  save  i  :  \igars2000\fisc\f  iscmat  \  * ,  filename  ( 1 :  end 
%  1  field  1 , filename { 1 : end-7 ) , 1  field ' ] } 


field 

7), 
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%READ_HDR 
^and  store 
Us  present 
file . 


read  and  interpret  information  from  EMSL  file 
the  pert inent  i n formation  in  global  variables 
ly  only  capable  of  handling  up  to  two  tasks  i 


header 

The  routine 


n  the 


lli/iO  -i-5  UC 


%  r  e  a  d__hd  r  ( f  i  1  e ) 

% INPUT: 

%file  --  directory  path  and  name  of  unprocessed  data  file  from  the 


%EMSL  HP. 

%The  function  returns  the  pertinent  variables  as  global 
%variables,  n a m e 1 y : 

^object:  object  name  as  identified  in  header 

Amnemonic:  object  mnemonic  name  as  identified  in  header 


%freq_start : 
%  f  req  jstop : 
%n__  freq: 

-If : 

%avg: 

%tx__sled: 

%rx  sled: 


start  frequency,  in  GHz 
stop  frequency,  in  GHz 
number  of  frequency  samples 
vector  of  frequencies 
number  of  averaging  points 
sled  used  for  transmitter 
sled  used  for  receiver 


axis  that:  is  moving 


%A__p°sit  ion : 
%B  position: 
% Deposition : 
& Deposition: 
%n_ang : 

%ang: 

%n_hdr : 

*hdr  : 


initial  position  of  Sled  A 

initial  position  of  Sled  B 

initial  position  of  rotation  axis 

initial  position  of  linear  translation  axis 

number  of  angular  samples 

vector  of  angular  movements/separation  of  relevant  axis 
number  of  lines  contained  in  the  header 
complete  text  of  header 


function  read__hdr  (f  ile) 
global  object 
global  mnemonic 
global  freq_start 
global  freq_stop 
global  n_freq 
global  f 
global  avg 
global  tx_sled 
global  rx_sled 
global  movement 
global  Ajposition 
global  B_j?osition 
global  Deposition 
global  L_position 
global  n_ang 
global  ang 
global  n_hdr 
global  hdr 


%read  header  again  to  get  information  about  the  EMSL  data  file: 

f close { 1  all f )  ; 
f id=fopen (f ile,  1 rt ' ) ; 
hdr=fscanf (fid,  f %c 200*72 ) ; 
hdr  size=size (hdr) ; 


%get  size  of  header  (line  4): 

n  head  1  =  72*3+18; 
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n_head_2  =  4*72; 

n__hdr  =  str2num  (hdr  (n_head_l :  n_head_2 )  )  ; 

% ge t  ob j  e  c t  a nd  mnemon i c ( 1 i n e  s  1  a nd  1 1 ) : 

mnemonic  =  hdr(7:72); 
object  =  hdr (72*10+10:72*11) ; 

%get  initial  axis  positions  (lines  28-31) : 

Apposition  =  str2num(hdr (27*72+1:28*72) ) ; 
B_position  =  str2num (hdr (28*72+1 : 29*72) ) ; 
D_position  =  str2num (hdr (29*72+1 : 30*72 )) ; 
L__position  =  str2num(hdr  (30*72  +  1:31*72)  )  ; 

%  get  LB/UHF  pa  rarnet  ers  : 
s _ f  =  1 Microwave  LB/UHF1 ; 

a  =  f indstr (s_f ,hdr) ; 

LB  =  hdr (a(2) :a(2)+24*72) ; 

f req_start  =  str2num (LB (72+1+12 : 2*72) ) ; 
freq_stop  -  str2num(LB (2*72+1+12 : 3*72) ) ; 
n_freq  =  str2num(LB (3*72+1+12 :4*72) ) ; 
f  =  linspace (freq_start , f  recast  op,  n_freq) ; 
avg  =  str2num (LB (6*72+1+12 : 7*72 )) ; 
tx_sled  =  str2num (LB (17*72+1+12 : 18*72) ) ; 
if  tx_sled  ==  1 

tXpSled  =  1 A 1  ; 

elseif  tXpSled  ==  2 

tXpSled  =  1  B 1  ; 

else 

warning (  1  This  is  a  monostatic  measurement') 

end 

rx_sled  =  str2num (LB (21*72+1+12 : 22*72 )) ; 
if  rXpSled  ==  1 
rXpSled  =  'A1; 

elseif  rx_sled  ==  2 
rx_sled  =  '  B 1  ; 

else 

warning ( 1  This  is  a  monostatic  measurement1) 

end 


% which  antenna,  if  any,  is  moving: 
s8  «  ’Sled  A  movement1; 
s9  =  'Sled  B  movement1; 
slO  =  'Axis  rotation1; 

a  =  f indstr (s8 , hdr)  ; 
b  =  f indstr (s9,hdr)  ; 
c  =  f indstr (slO, hdr) ; 

if  n_hdr==82 

movement  =  ' none  1 ; 
n_ang  =  1 ; 
ant__start  =  'N/A'; 
antpStop  =  'N/A'; 

%get:  antenna  movement  parameters: 
elseif  -isempty(b) 
movement  =  '  B  '  ; 

SLED  =  hdr (b(2) :b(2) +4*72); 
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el seif  -isempty(a) 
movement  =  1 A 1  ; 

SLED  =  hdr (a ( 2 ) : a (2) +4*72) ; 
elseif  -isempty(c) 
movement  =  1  D '  ; 

SLED  =  hdr (c(2)  : c  ( 2 ) +4  * 72 )  ; 

end 

if  n_hdr>82  %if  the  measurement  is  more  than  a  single  frequency  sweep, 

mov__start  =  str2num  (SLED  ( 72  +  1  +  13  :  2*72)  )  ; 
mov_stop  =  str2num (SLED (2*72  +  1  +  13  :  3*72 )) ; 
n_ang  =  str2num (SLED (3*72+1+13 : 4*72 )) ; 

end 

% interpret  ante nna  movement  pa  r amet ers : 

if  tx_s  1  ed= = rx_s  1  ed  &  n_ang~=l  %if  the  measurement  is  monostatic, 

ang  =  linspace (mov_s tart ,mov_stop, n_ang) ; 
i f  movement = = 1 A ' 

initial  =  A_position; 

ang  =  linspace (mov_start, mov_stop,n_ang) ; 
elseif  movements  1  DT 

initial  =  D_j>osition; 

ang  =  linspace (mov_s tart , mov_s top, n_ang) ; 
elseif  movement == r none 1 

ang  =  abs (A_position  -  Exposition); 

end 

else  %if  the  measurement  is  bistatic, 

if  movement == 'A' 

initial  =  Bjposition; 

ang  =  linspace (abs ( initial -mov_start ), abs (initial - 
mov^stop) , n_ang) ; 

elseif  movement== 1 B 1 

initial  =  Apposition; 

ang  =  linspace (abs ( initial -mov^start) , abs ( initial - 
mov_stop) , n_ang) ; 

elseif  movement== * D 1 

initial  =  Deposition; 

ang  =  linspace (mov_s tart ,mov_s top, n^ang) ; 
elseif  movement== T none  1 

ang  =  abs (Apposition  -  B_position) ; 

end 

end 

^Christopher  J  Bradley 
%2d  Lt,  USAF 
% AFX T/ ENG 
%  c  b  r  a  d 1 e  v @ i e  e  e . org 
%29  September  2000 
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%SCATJPLOT_ 2D  Plot  normalized  2-D  bistatic  RCS  with  a  dynamic  range 
i of  80  dBsm. 

Ascat__plot_2D  (M,  pol, x,  domain} 

% INPUTS: 

%M  --  input;  scattering  matrix  (2  x  2  x  n_£req  x  n_ang) 

%pol  desired  polarization  to  plot 

%x  --  vector  of  values  of  x-axis 

%  domain  --  string  with  which  to  label  x-axis 

function  scat_plot_2D (M, pol , x, domain) 

global  n_freq 

global  n_ang 

global  mnemonic 

global  object 

if  pol—'hh1 
r  =  1; 

C  =  1; 

elseif  pol==,vh’ 
r  =  1; 

C  =  2; 

elseif  pol== ’ hv 1 
r  =  2  ; 

C  =  1; 

elseif  pol== ' w 1 
r  =  2; 

C  =  2  ; 
else 

error ( T Unrecognized  variable  "pol"  ') 

end 

figure (1) 

theta  =  1 inspace (4 , 184 , n_ang) ; 

fig  =  abs (squeeze (M(r, c, :,:)))  ; 

fig_norm  =  20*logl0 (fig/max (abs (M (:)))) ; 

imagesc (x, theta, rot90 (fig_norm) , [-80  0] ) ; 

xlabel (domain) 

ylabel ( ’bistatic  angle ’ ) 

title ( [object (1 : 50) , pol,  1  pel’,1,  1 , domain]); 

^Christopher  J  Bradley 
%2 d  Lt,  USAF 
?.AFIT/ENG 

% cbradl ey@ i eee . erg 
%2 9  September  2000 
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between 


%TFCONVERT  convert  bistatic  data 
% S_t f  =  t f con ve rt  ( S ,  f ,  t ,  ) 

% INPUTS: 


time  and  frequency  domain 


*S  —  input  scattering  matrix 

%f  —  vector  of  frequency  values  (in  GHz) 

%t  —  vector  of  time  values  (in  ns) --same  length  as  "fM 
%t  or  f  --  to  con  vert  to  time,  t___or__f  -  1  time  ’ 

%  to  convert  to  frequency,  t__or_f  =  ' freq* 

%OUTPUTS: 

%  S__t  f  -  -  t  r  an  s  form-domain  s  ca  1 1  er  ing  mat  r  lx  ( iden  t  i  ca  1  dime  ns  i  ons 
las  input  structure  "S" } 


function  S_tf  =  tf convert (S , f , t , t_or_f) 
global  n__ang 


if  t_or__f  ==  ’time1 
sgn  =  1; 
vecl  =  f; 
vec2  =  t; 

elseif  t_or_f  ==  1 freqT 

sgn  =  - 1 ; 
vecl  =  t; 
vec2  =  f ; 
else 

ERROR ( 1  Unknown  input  argument  for  "t_or_f " :  Enter  -1  or  +1 ’ ) 

end 

q  =  waitbar(0,  ’Computing  Chirp-z  transform...’); 
for  m=l:n_ang 
for  n=l : 2 

for  p=l : 2 

S_tf (p, n, : ,m)  =  mchirpz (S (p, n, : ,m) , vecl, vec2 , sgn) ; 
end 

end 

waitbar (m/n_ang)  ; 

end 

close (q) 


^Christopher  J  Bradley 
%  2d  Lt,  USAF 
%AFIT/ENG 

icbradley@ieee . org 
%29  September  2000 
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%TIME__GATE  apply  time  gate  to  data  using  "mchirpz .m" 

%S  -----  time_gate  (M,  f,  t) 

% INPUTS : 

%M  --  input  data 

%£  --  vector  of  frequency 

%t  —  vector  of  time  {desired  gate) 

% OUT PUTS : 

%M_f re  --  gated  frequency-domain  data 
%M  tim  --  gated  time-domain  data  (optional) 

function  [M_fre,  M_tim]  =  time_gate (M,  f , t ) 

global  n_ang 

global  gate_start 

global  gate_stop 

gate_start  =  t(l); 

gate__stop  =  t(end); 

n  =  length (f) ; 

% create  Hamming  window: 

H  =  hamming (n) ’ ; 
big_H  =  zeros (2 , 2 , n, n_ang) ; 
for  p  «  1:2 
for  n  =  1:2 

for  m  =  l:n__ang 

big_H(p,n,  :  ,m)  =  H; 
end 

end 

end 

%  w i ndo w  i n  f r eq  u en c y : 

M_win  =  M.*big_H; 

% t r a n s  f o r m  t o  t ime-d oma i n : 

q  =  waitbar(0,  ’Applying  gate...1); 
i=l; 

for  p=l : 2 
for  n=l : 2 

for  m=l:n_ang 

M_tim (p, n, : ,m)  =  mchirpz (M_win (p, n, m) ,  f , t , 1) ; 
waitbar  ( i/4/n__ang/2  ,  q) 
i=i+l ; 

end 

end 

end 

^transform  back  to  frequency-domain : 

for  p=l : 2 
for  n=l : 2 

for  m=l:n_ang 

M_win_ga (p, n, : ,m)  »  mchirpz (M_tim(p/n/ : ,m) , t, f, -1) ; 
waitbar ( i/4/n_ang/2 , q) 
i = i + 1  ; 

end 

end 

end 

close (q) 
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%  un do  H  ammi ng  w i  ndo  w : 

factor= (t (1) -t (end) )*(f(2)-f(l)) 
M_fre  =  f actor *M_win_ga. /big_H; 

%  C  h  r  i  s  fc  o  phe  r  J  Bra  d  1  e  y 
%2d  Lt ,  USAF 
%  AFT.  T /ENG 
ncbradleyQieee ♦ org 
%  2  9  S  ep  t  ember  2  0  0  0 


APPENDIX  C:  The  EMSL 

http://www.tdp.sai.irc.it/TDP/test  and  evaluation  facilities.htm 


Radar  measurements 


The  European  Microwave  Signature  Laboratory  (EMSL) 


1 .  Main  Door 

2.  Service  door 

3.  Object 

4.  Object 
Positioner 

5.  Object 
Positioner  rail 

6.  Radar  Sleds 

7.  Sled  rail 

8.  Fixed 
Antennae 

9.  Free  standing 
structure 

10  .  Microwave 
absorber  cover 


Table  1:  Positioner  Mechanical  Specifications 


Positioner 


Range  Repeatability  Max  Load 


Radars 

Object  (linear) 

Object  (rotation) 

Table  2:  Measurement  Systems  Specifications 


-110/110°  +/-  0.005°  50  Kg 

-2. 5/2. 5  m  +/-  0.5  mm  6000  Kg 

0/360°  +/-  0.05°  6000  Kg 


Measurement  configuration  mono/bistatic/multistatic 

Frequency  Range  0.5  -  26.  5  GHz 
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Polarisation 
Measurement  mode 
Detection 
Sensitivity 
Dynamic  Range 
Transmitted  power 


Full  polarimetric  (HH,  VH,  HV,  W) 
Continuous  Wave  (Step  or  Ramp) 
Coherent  (phase  locked) 

-60  dBsm 
100  dB 
<  1  W 


The  European  Microwave  Signature  Laboratory  (EMSL)  of  the  Space  Applications 

Institute  at  the  Joint  Research  Centre  at  Ispra  (Italy)  is  designed  for  experiments  in  the 
field  of  radar  signature  research.  It  is  a  unique  European  facility  which  complements  in¬ 
field,  air-  and  space-  borne  radar  measurements  by  providing  stable  and  reproducible 
environmental  conditions  and  operational  modes  for  well  controlled  experiments  of 
microwave  scattering. 

The  EMSL  consists  of  a  large  anechoic  chamber,  providing  the  environment  and  the 
facilities  to  perform  radar  measurements.  The  measurements  essentially  can  provide  the 
spatial  scattering  pattern  of  the  object  under  test  in  a  three-dimensional  characterization. 
The  radar  scattering  matrix  is  measured,  which  represents  the  Microwave  Signature  of 
the  object,  the  most  comprehensive  information  from  active  microwave  measurements. 
Other  widely  used  parameters  such  as  the  Radar  Cross  Section  (RCS)  can  be  directly 
derived  from  the  scattering  matrix.  The  variety  of  measurement  options  can  be  used  in 
two  basic  ways:  scatterometric  and  imaging  modes. 

The  facility,  in  addition  to  the  monostatic  configuration  of  usual  radar  systems, 
offers  polarimetric  measurement  features  in  bistatic  and  multistatic  configurations.  This 
gives  additional  dimensions  in  the  measured  data  for  characterizing  an  object,  validating 
a  model,  or  testing  of  new  sensors. 


143 


The  overall  structure  is  made  of  a  hemispherical  and  a  cylindrical  part,  both  with  a 
diameter  of  20  m.  A  circular  rail  is  mounted  between  the  two  parts  on  which  two  sleds 
carrying  antennas  and  other  sensors  are  moved  independently.  Further,  37  receiving 
antennas  are  fixed  on  the  inner  surface  of  the  hemisphere  and  are  integrated  in  the 
microwave  measurement  system. 

Objects  of  varied  sizes  and  shapes  can  be  located  on  the  object  positioner.  The 
positioner  moves  on  a  rail  to  allow  the  precise  linear  and  rotational  positioning  of  the 
object  under  test  inside  the  chamber  for  the  microwave  measurements.  The  spatial 
resolution  depends  on  the  measurement  mode  and  parameters;  it  ranges  from  few  meters 
to  less  than  1  cm. 

The  electro-mechanical  devices  and  the  microwave  measurement  system  are 
remotely  controlled.  The  control  program  allows  user-defined  sequences  of  operations  to 
create  measurement  scenarios,  including  microwave  measurements  and  data  storage, 
sensor  positioning  and  object  movements.  Complex,  extended  measurement  sequences 
may  be  executed  automatically  by  this  control  system. 

All  collected  data  (microwave  data,  system  parameters,  environmental  conditions, 
optical  images)  are  transferred  to  the  EMSL  Information  Management  System.  This  data 
archiving  and  handling  system  was  especially  developed  to  cope  with  the  complex 
information  management  for  all  data  generated  by  the  EMSL  activities  and  the  related 
research. 

The  EMSL  was  first  applied  to  the  investigation  of  the  use  of  radar  to  counter  the 
APL  threat  some  three  years  ago  following  a  request  from  the  European  Parliament. 
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